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Об'єкт дослідження - математичні моделі геологічного середовища та класи єди- 
ності розв'язку обернених задач теорії потенціалу. 

Мета роботи - розробка теорії і методів розв'язування лінійних і нелінійних обер- 
нених задач для тіл з відомими класами єдиності і для тіл, близьких до цих класів. 

Метод дослідження - формально-теоретичний аналіз існуючих підходів щодо пос- 
тановки і методів розв'язання граничних обернених задач теорії потенціалу. 

В роботі представлено нові теоретичні результати дослідження єдиності та стій- 
кості локальних розв'язків конкретних інтегро-диференціальних рівнянь, до яких зво- 
дяться обернені задачі гравіметрії, що пов'язані з визначенням форми чи густини ано- 
малієстворюючих тіл. Розроблено основи ,лінеаризованої" постановки оберненої стру- 
ктурної задачі логарифмічного потенціалу в класі неперервних, разом з першими похі- 
дними, функцій для границі поділу однорідних пластів різної густини за умов задання 
вхідних даних на істотно обмежених множинах. Побудовано чисельні алгоритми з мож- 
ливістю керувати швидкістю збіжності відповідних ітераційних процесів і проведено 
порівняльне тестування для конкретних математичних моделей середовища. 

Розв'язано обернену задачу логарифмічного потенціалу для обмеженого тіла, 
близького до круга, за відомої постійної щільності без істотних обмежень на розмірність 
функціональних просторів. Запропоновано певні методи непрямого визначення коефі- 
цієнтів характеристичного полінома матриці, що засновані на ортогоналізації деяких 
послідовностей векторів. Охарактеризовано області застосування даних методів при 
відшуканні розв'язків систем лінійних рівнянь, до яких найчастіше зводяться обернені 
задачі геофізики. 


ВЛАСНІ ВЕКТОРИ, ВЛАСНІ ЗНАЧЕННЯ, ГРАНИЧНІ ЛІНІЙНІ ТА НЕЛІНІЙНІ, 
ПРЯМІ І ОБЕРНЕНІ, КОРЕКТНІ І УМОВНО КОРЕКТНІ ЗАДАЧІ ТЕОРІЇ ПОТЕНЦІАЛУ, 
КОНТАКТНА ЗАДАЧА, СИСТЕМИ ЛІНІЙНИХ РІВНЯНЬ, ТЕОРЕМИ ІСНУВАННЯ, 
ЄДИНОСТІ ТА СТІЙКОСТІ РОЗВ'ЯЗКІВ, ЧИСЕЛЬНЕ МОДЕЛЮВАННЯ. 
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ВСТУП 


Дослідження, результати яких викладаються в даному п'ятирічному звіті, продовжу- 
ють розпочате кілька десятиліть тому назад вивчення обернених задач теорії потенціалу 
100-105) та містять розробки ряду нових проблем теорії інтерпретації гравітаційних 1 магні- 
тних аномалій |18-21, 78-91, 97-99, 106-109, 111-114|. Можна з впевненістю сказати, що 
створено основи теорії граничної коректності певних лінійних і нелінійних обернених задач 
теорії потенціалу, чільне місце в якій займає вивчення проблем єдиності розв'язків. Важли- 
вість вивчення цих проблем не піддається сумніву математиками, бо тільки коректні в кла- 
сичному сенсі чи в сенсі коректності за А. М. Тихоновим задачі фізично детерміновані. Про- 
те в геофізичній літературі, поряд з широким застосуванням теорії регуляризації, де, до речі, 
проблеми єдиності одні із центральних, не приділялось належної уваги цій проблематиці. Бі- 
льше того, наприкінці 70-х років минулого століття, В. М. Страхов висунув концепцію є- 
еквівалентності в задачах інтерпретації гравітаційних 1 магнітних аномалій. Розвиваючи цю 
концепцію, О. І. Кобрунов довів її до такого рівня, у відповідності з яким розв'язки задач з 
цілком неперервними операторами слід шукати не в компактних класах єдиності, а в більш 
потужних екстремальних класах, що містять (або ж 1 ні) класи єдиності як власні підкласи. 
Безпідставність цього твердження особливо чітко проглядає з його логічного висновку, у ві- 
дповідності з яким ,дбажані" розв'язки обернених задач здобуваються в нуль-множині опера- 
тора задачі. Інакше кажучи, відповідні цим розв'язкам красиві геологічні побудови знахо- 
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дяться за значеннями нульового поля, тобто, за його повної відсутності. Тому вивчення про- 
блем єдиності залишається актуальним в загальній теорії регуляризації обернених задач гео- 
фізики аж до сучасності. 

В даній роботі наведено результати дослідження єдиності локальних розв'язків конк- 
ретних інтегро-диференціальних рівнянь, до яких зводяться ті чи інші обернені задачі, 
пов'язані з визначенням аномалієстворюючих тіл в рамках теорії малих збурень. Також роз- 
роблено основи теорії ,лінеаризованої" оберненої задачі логарифмічного потенціалу для ко- 
нтактної поверхні в локальній області за умов задання вхідних даних на істотно обмежених 
множинах без обмежень розмірності просторів даних та розв'язків в рамках загальної теорії 
наближеного вирішення умовно коректних задач геофізики; продовжено дослідження прик- 
ладних аспектів задач відновлення границі зірчастих тіл, близьких до заданих; запропонова- 
но метод відшукання розв'язків систем лінійних рівнянь, в якому використовуються процеси 
доортогоналізації деяких сімей векторів, що виникають при побудові характеристичних по- 
ліномів матриць. 

За результатами тематичних досліджень відбувся захист двох кандидатських дисерта- 
цій: Якимчиком А. І. ,Гранична задача відновлення потенціалу за значеннями модуля його 
градієнта" (2001 р.) 1 Дубовенком Ю. І. ,,Визначення контактної границі за значеннями похі- 
дних логарифмічного потенціалу на істотно обмежених множинах" (2005 р.). 


1. ДОСЛІДЖЕННЯ ЗАДАЧІ ВИЗНАЧЕННЯ КОНТАКТУ ЗА ПОЛЕМ, 
ЗАДАНИМ НА КОРОТКОМУ ПРОФІЛІ 


Вступ. Теорія розв'язання плоских обернених задач потенціалу для контактної повер- 
хні з інтегральними операторами Пуассона сформована 1 вивчається здавна в численній гео- 
фізичній літературі, як то (39, 44, 51, 57, 1191, яка періодично оновлюється. Вирішення кон- 
тактних задач грунтується на апроксимації, запропонованій ще Андреєвим |7|. В (103) отри- 
мано набір лінійних інтегральних рівнянь, якими на класі Ми (ай, п) меж С (х) поділу одно- 
рідних середовищ з високою точністю замінено нелінійне інтегральне рівняння оберненої 
задачі потенціалу для контактної поверхні і побудовано збіжні ітераційні процеси для їх 
розв'язання. Згідно |103|, нелінійна обернена задача логарифмічного потенціалу для контак- 
тної поверхні зведена до розв'язання нелінійного інтегрального рівняння типу Урисона 

ПК А ЖАО, 0 (б) 
2-2 Па ; 45 з и(х), (х)о ОТ (11) 
Зп, (8-х) «СЧЕ) 2 
де, на відміну від рівняння Б. В. Нумерова |44| дано замість середньої глибини й до контакту 
значення 0 (х) цього контакту у фіксованій точці х (можливість такої підстановки вказав 





О. А. Шванк (110), а обгрунтував А. В. Чорний |96, 1051). Розв'язок цього рівняння знайдено 
на компактному структурному класі Страхова!, на якому розв'язок нелінійної контактної за- 
дачі (1.1) існує, єдиний і стійкий |103|, але досягається лише при заданні значень поля на не- 
обмежених інтервалах, а в польових умовах елементи поля вимірюють на істотно коротких 
профілях. В (21| вивчалась контактна задача за умов задання вхідних даних на істотно об- 
межених множинах як лінеаризований варіант задачі |76|, який дозволяє відновити контакт 
за значеннями поля на істотно обмежених інтервалах через заміну інтегрального рівняння з 
повільноспадним ядром Пуассона рівнянням зі швидкоспадним ядром. Способи її коректно- 
го розв'язання з урахуванням результатів (21, 97, 103| наведено в |18|. Але досі не з'ясовано 
питання, з якими особливостями задачі пов'язано розмаїття її описів: чи то з неозорою мно- 
жиною моделей задачі; чи то з наявністю теоретичних проблем її розв'язності, що важко під- 
на практиці. Відсутність у літературі детального аналізу постановок задачі і методів її 
розв'язування, з якого випливала б відповідь на поставлене питання, спонукало до спроби 


1 Перетині класу функцій густини, обмежених і локально інтегрованих майже всюди в необмеженій замкнутій області, що 
містить тяжіючі маси, і класу неперервних за Гельдером (обмежених) функцій контакту. 
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з'ясувати ситуацію, що склалася. Поза спробою вичерпного огляду всієї літератури з про- 
блем коректної розв'язності 1 чисельної реалізації розв'язків задачі для контактної поверхні 1 
ефективності їх використання для тлумачення гравімагнітних аномалій в статті |97| зроблено 
якісний аналіз основних праць щодо стану вивченості математичної моделі задачі, який по- 
тім продубльовано в |21|. Наводимо звідти основні міркування. 


1.1. ОСНОВНІ МОДЕЛІ ЗАДАЧІ. Задачу про відновлення межі поділу двох однорідних 
середовищ уперше сформулював і наближено розв'язав Б. В. Нумеров в 1930 р. у такій пос- 
тановці |44). Відомо, що під плоскою поверхнею Землі ОД: у «0 є одна контактна границя 


др. ібес (х), -о0« хаос, що розділяє однорідні шари легких "верхніх " порід з густиною 
с, з соп5і від важчих "нижніх" порід з густиною с» - соп5і . Окрім різниці є є с, - б» гус- 
. . . . ж . 
тин порід задано значення сили тяжіння (/, (2,0), -со«х «со, і віддаль й - Дх') між пове- 
рхнями 0, і дО, у певній точці х , яку відраховують за нормаллю до дО, (рис. 1.1. а). За 
цими вихідними даними потрібно визначити "поверхню" ОР,, припустивши, що невідома 
функція Є «С(х), яка описує межу поділу, є однозначною й обмеженою. Для визначення ко- 
нтакту в цих умовах Б. В. Нумеров вивів нелінійне інтегральне рівняння 
п 2 2 
1 о ВСІХ Й (20) 
ре | Я З РЧР ДО и(х)- -е 
2 (5-х) я с 


що залежить від двох параметрів - стрибка є густини на контакті 0 (х) і глибини Й залягання 


12) 


контакту в деякій довільно обраній точці х', яка не збігається, взагалі кажучи, з точкою х 
спостереження поля; 7 - гравітаційна стала. Наближений розв'язок цього рівняння Б. В. Ну- 
меров отримав у вигляді 

сн (х) РЕ Ані) 
за додаткової умови, що рельєф контактної границі не надто контрастний, в крайньому випа- 
дку такий, що коливання функції ДЛ ,0 (х) менше від /. За цієї умови однозначну функцію 


Йон (х) визначено як обмежений розв'язок лінійного інтегрального рівняння першого роду 


175 Й М 
Пе руайнєкои М 


яке залежить від двох параметрів є 1 й, тобто Ані) а пуд ЕЯ о; п). 





М. Р. Малкін |39| задачу визначення контактної поверхні С), розділу двох однорід- 
них косо намагнічених шарів за заданим магнітним потенціалом і стрибком ./ з у хм ,) век- 


тора намагніченості порід звів до розв'язання нелінійного інтегрального рівняння, яке за вер- 
тикальної намагніченості порід збігається з рівнянням (1.2). Наближений його розв'язок 
ак (х)ел зара (х) одержано методом Нумерова і тому залежить від трьох параметрів: 
двох природних - складових ./, і ./, стрибка вектора намагніченості порід і одного "привне- 
сеного? | - | товщини Й виділеного однорідного магнітного шару, тобто 
А ЗАЄГ ав 

Г. Рейнбой П1138), порівнюючи лінеаризоване рівняння (1.3) для контактної поверхні з 


інтегралом Пуассона для подання гармонічної функції м(Х, у) в півплощині за значенням її на 
рівніу - 





Мо ИН 
і сего 


виявив, що лінеаризована задача для контактної поверхні еквівалентна задачі визначення 
гармонічної функції у(х, у) за значеннями у(х,0), тобто що А/(х)- м(х, п). Тим самим він 
фактично показав, що точний розв'язок інтегрального рівняння (1.3) повинен бути слідом 
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гармонічної функції (хоча цей факт ще довго залишався неусвідомленим) 1, окрім того, відк- 
рив перспективний напрямок в теорії інтерпретації аномалій, пов'язаний з аналітичним про- 
довженням полів у півпростір зосередження їх джерел (не тільки необмежених типу шару, 
але Й обмежених) | 37, 59, 70, 1191. 

О. О. Заморєв |22, 23| дослідив у комплексній множині 2 -х-кіу задачу визначення 
контактної поверхні ОД, у вигляді скінченної ундуляції у(х) є й - п), уУ(а)- у(Б) є Й, на скі- 
нченному інтервалі |а,Б| горизонтальної прямої у(х) «Йесопії » 0, хе Їа, ь| (рис. 1.1. б). 
Проблему відновлення форми збурюючого тіла (виповненого масами однорідної густини с і 
обмеженого, на відміну від з нумерівського ї шару, прямою  у-й і кривою 
У(х) я ячнт(х), х є|а, В|, па)- т) - 0), тобто задачу визначення параметрів а,Р, однозна- 
чної на інтервалі |а,Р| обмеженої функції лох) за заданими на осі у - 0 значеннями (відпові- 
дним чином прокаліброваними) похідної логарифмічного потенціалу (/ б» у) Заморев звів до 
нелінійного рівняння, яке можна записати в двох еквівалентних один одному зображеннях: 


і Є-хчій о ічп(Є) 1 М 
не оо Б оннеотоеруз гарно (14) 


дочка се)- 90) 1 92», у) :дЧ(х,») 


л/с да 2 Ох бу 


Обидва ці подання наведено у статті |22|, до того ж друге, яке детально вивчено, оде- 
ржано з першого інтегруванням частинами за умови та) б ть) - 0. В результаті аналізу дру- 











гого подання напруженості и(2) поля Заморєв знайшов його точне обернення. На жаль, ефек- 
тивність цього обернення при розв'язуванні прикладних задач не перевірено ще й досі. Звер- 
немо тепер увагу на перше зображення напруженості поля. Неважко виявити, що дійсна її 
складова 





оо Че в, 
По 

2. (5-х) "(5) 
за зовнішнім виглядом нагадує рівняння Б. В. Нумерова (1.2) без доданку Й в правій частині, 
що відтворює за вибраної калібровки поля притягання однорідного шару товщиною /?0. Бі- 
льше того, величина інтегралу в лівій частині рівняння (1.5) не змінюється на заморєвських 
тілах у разі заміни скінченних границь а і Б на нескінченні -оо і со, тому ліві частини рівнянь 
(1.2) 1 (1.5) за зовнішнім виглядом повністю збігаються, хоча і зображають в математичній 
символіці різні за суттю предметні моделі: перше - притягання необмеженим шаром між 
границями у-й і у-с (х), - 00 « х 200; друге - притягання скінченним тілом, обмеженим 


ає зи(х), У(х)е пт), х а, | (1.5) 


кривими у- й і уз Ачн(х), х Є а, |, та) (Б) - 0. Можливо, зовнішня подібність лівих 
частин рівнянь Нумерова й Заморєва відігравала свою роль у подальших прикладних дослі- 
дженнях, серед яких осібно вирізняється праця |110). 

У цій праці, написаній в 1941 р. і виданій в 1947 р., О. О. Шванк помітив, що коли в 
рівнянні Б. В. Нумерова замість параметра Й підставити значення контакту С (х) в точці х 
спостереження поля, то з рівняння (1.2) виходить нове нелінійне інтегральне рівняння 


ДНО СО НН 
пара б 


яке за не дуже великих коливань шуканої функції 0 (х) можна розв'язувати методом послі- 





довних наближень. На одному окремому прикладі він показав, що обчислені "вручну" послі- 
довні наближення 12, (х)); контакту 2(х) з першого 4, (х)- ц(х) і до 6-го включно мають те- 
нденцію до збіжності. Проте зауваження О. О. Шванка, засноване на виводі Б. В. Нумерова, 
залишилось, на жаль, без уваги, можливо, через те, що і в самого автора пропозиції не скла- 
лося чіткої думки про єдино можливий функціональний зв'язок поміж контактом і полем. 
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Дійсно, в 5 2 розд. 8 (стор. 361-374) його книги (110) для поля від контакту розгляда- 
ються, як рівноправні, два вирази, один з яких, а саме, 


тезу оббіє со ЗАоДІЗИ ал) 


являє собою лінеаризацію зв'язку 





що АН 
- Га (є з) Бл 7 Рой (18) 
ях, (з-х) х6"(з) 
а другий - зв'язок вигляду (1.6). Так чи інакше, але в наступних роботах О. К. Маловичка 
140), В. М. Новоселицького (|43| 1 багатьох інших авторів (див. Бібліографію в |971) за рівнян- 
ня для контактної поверхні брали вираз (1.3) або його лінеаризацію (1.7). Кожен з цих вира- 
зів залежить від двох параметрів - сі п - й описує, як це випливає з інтегрального рівняння 
(1.5) Заморєва при а-»-о і Б -» со, вплив необмеженого в латеральному напрямку і обме- 
женого за товщиною геологічного тіла на зразок синклінорію, антикліналі й синкліналі якого 
виповнені масами однорідної густини протилежних знаків (рис. 1.1. в). Тіло це відділено від 
рівня спостережень шаром товщини /, дія якого виразом (1.8) не враховується. Разом з тим, 
у цій постановці було отримано цікаві результати, які бажано використати в теорії структур- 
ної задачі, що грунтується на рівнянні (1.6). До числа таких розробок в першу чергу відно- 
сяться відкриті М. Г. Сербуленком |52| 1 Ю. В. Антоновим |б| способи визначення двох ї бі- 
льше контактних границь, що розділяють однорідні шари. Слід згадати також дослідження 
О. К. Маловичком |40| питань уточнення розв'язків лінеаризованого рівняння та єдиності 
його розв'язків, якими встановлено, що одному й тому ж полю можуть відповідати як об- 
межені однорідні випуклі тіла, так і різноманітні необмежені однорідні тіла типу шару. За- 
слуговують на увагу перші результати чисельного розв'язку нелінійного інтегрального рів- 
няння (1.8) у межах теорії регуляризації (14, 16, 24, 64, 651, вивчення В. І. Старостенком |53| 
можливостей стійкого чисельного визначення декількох контактних границь, що описуються 
нелінійними рівняннями виду (1.8), і спроби О. І. Кобрунова та Р. П. Денисюка |31-32| від- 
шукати за допомогою одного чисельного алгоритму границі розподілу неоднорідних компле- 
ксів. 
Аналогічну модель задачі для контактної поверхні в комплексній площині 2 -х-ніу 


запропоновано в 70 - х р. р. минулого століття В. М. Страховим |60). Саме для шаруватого 
середовища у смузі, відокремленій від нуля і обмеженій прямими у-Пп 1 узи, 
йо «й" «0, що складаються з двох однорідних шарів порід з густинами б, і б», які розділяє 
жорданова крива дД:препі(х), - со « х «со зі стрибком с з с, - є, густини на контакті і з 


асимптотою й, й «П«й"', однією й тією ж ліворуч та праворуч, яку крива ОЇ) може пере- 
тинати не більше, ніж скінченне число разів (рис. 1.1. г), він вивів нове нелінійне інтеграль- 
не рівняння 
5 - 5 - ЗіЙ и"(2) г2є р" 
Г---бе (г) (19) 


ня МЕРА: 
СЕ и (2), зер 
де Р)" - частина площини "вище" кривої 2), ЩД)' - "нижче" кривої ОД), 5-счій, 
5 «6 -її; ца) - в калібровці (1.4) задано за результатами спостережень на "земній поверх- 


ні? у-0. Для цього рівняння Страхов довів теорему єдиності його розв'язку 1 редукував 
задачу визначення контакту за полем до задачі побудови допоміжної функції, що конформно 
(1 однолистно) відображає нижню півплощину певної комплексної площини на область 1) 
під шуканою кривою. Цей підхід до розв'язування задачі передбачалося побудувати за ана- 
логією з розвинутим В. К. Івановим методом відшукання однорідно розподілених джерел на 
скінченних однозв'язних областях |26|. Незабаром після вибраної редукції було знайдено 
цілі класи розв'язків рівняння (1.9), Їх подавали або як сімейства розв'язків у скінченному ви- 
гляді, коли границю шуканої області визначають скінченним числом параметрів зі скінчен- 
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ної системи (нелінійних) рівнянь, або у вигляді двопараметричних сімейств розв'язків "6 ма- 
лому", коли границю визначають у вигляді (невеликих) відхилень контура шуканої від кон- 
тура заданої області. Ці класи отримали, відповідно, О. В. Цирульський |71, 72|1 В. Г. Чере- 
дніченко |74|. До того ж В. Г. Чередніченко зняв обмеження на розв'язки рівняння (1.9), до- 
пустивши, що крива 0) необов'язково має бути графіком і може перетинати асимптоту Й 
нескінченну кількість разів. Детальніше про це можна довідатись в бібліографічному комен- 
тарі дисертаційної роботи автора цього розділу, автореферат якої вказаний в переліку поси- 
лань як |211|. 


1.2. МОДЕЛЬ І ПОСТАНОВКА КОНТАКТНОЇ ЗАДАЧІ. Для опису фундаментальної залеж- 
ності гравітаційного поля від шаруватого середовища мались на початок 1980 - х років дві 
принципово різні математичні моделі: модель Нумерова та модель Маловичка - Страхова 
(рис. 1.1. а, г). У моделі Нумерова як прикладну модель середовища розглянуто необмежений 
шар між контурами ОР, та ОР, його притягання враховано у вигляді суми притягання ша- 
ром постійної товщини Й і відхиленнями від цього шару між прямою у- й ікривою у-( (х) 
. У моделі Маловичка - Страхова предметну модель шаруватого середовища ототожнено з 
флуктуаціями у - лох) - п кривої уз ло) відносно своєї асимптоти ЛП, що обмежують "дода- 
тні" та "від'ємні" маси густиною с; поле цих мас аж ніяк не могло збігатися з полем моделі 
Нумерова, оскільки не було враховано вплив шару товщиною Л. Проте кожне з інтегральних 
рівнянь (1.2) 1 (1.9), що відповідало тій чи іншій моделі, виявилось однаково залежним від 
двох параметрів: надлишкової густини с 1 сталої товщини (асимптоти) Пп шару. Відповідно і 
розв'язки рівнянь також стали залежними від цих параметрів. У літературі відзначений па- 
раметричний зв'язок подано навіть як визначальний для класифікації обернених задач граві- 
метрії. З виходом у світ праць |60, 71-74, в яких закладено основи теорії оберненої задачі 
для предметної моделі Маловичка - Страхова (рис. 1.1. г), устоялося положення, що обернені 
задачі для обмежених тіл принципово відрізняються від обернених задач для необмежених 
тіл. Їх відмінність полягає саме в тому, що розв'язки задач для обмежених тіл залежать тіль- 
ки від одного параметра с, тоді як розв'язки для необмежених тіл залежать як від параметра 
с , так і від параметра ЛП ("середньої глибини" за Нумеровим або асимптоти за Страховим). І 
довгий час ніхто не піддавав сумніву переконливість цього твердження, хоч воно явно не уз- 
годжувалося з аналізом рівнянь (1.9) та (1.6). Уже побіжний погляд на них переконує, що рі- 
вняння (1.9), а отже, й його розв'язок має залежати як від параметра с, так і від асимптоти Й 
шуканого контура 0), тоді як рівняння (1.6) для того ж контура (як і численних інших обме- 
жених) від параметра / не залежить. Рівняння (1.6), якщо воно вірне, давало зразок опису за- 
дачі для контактної поверхні, для якого одне з положень теорії, розвинутої в згаданих робо- 
тах, потребувало уточнення. З невідомих причин цього не сталося. Або автор пропозиції (1.6) 
чи його послідовники не були впевнені в адекватності опису моделі задачі, або авторам 
праць |60, 71 - 74| не була відома праця (110). Так чи інакше, але рівняння (1.6) було проїг- 
норовано під час побудови теорії обернених задач для шаруватих середовищ 1 питання про 
вирішальну відмінність задач для обмежених та необмежених тіл залишилося відкритим. 
Окрім того, в створеній теорії не розглянуто узагальнення відкритого Г. Рейнбоєм факту, що 
розв'язки лінеаризованого інтегрального рівняння для контактної поверхні слід шукати на 
обмежених множинах гармонічних функцій. Математичне обгрунтування згаданих та інших 
супутніх проблем стало предметом подальших досліджень |76, 96, 103, 104). 

Початок цих досліджень було закладено аналізом предметних моделей Нумерова і 
Маловичка - Страхова. В результаті адекватною для двовимірною евклідового простору 


в в в), в - (- оо, го) виявилась предметна модель Нумерова, яку подано у вигляді 
шару 

рейху) -о«хао; 04у«х)), 
обмеженого контурами ОД): у-010Ю,: у-с (х), хек), зануреного в простір КО?) і випо- 
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вненого масами постійної густини с «ос, - о, (рис. 1.2). Відповідну предметній математич- 
ну модель контактної поверхні знаходять за допомогою регуляризації відповідного означено- 
го інтегралу, що розбігається на границі С) - дО), (/ СД, збурювальної області Р. Регуляри- 
зацію здійснено за правилом |76, 103| 








с Є) 





Рис. 1.2. Предметна модель: контактна поверхня порід з різною густиною. 


Ця операція теоретично обтрунтовує адекватність нелінійного інтегрального рівняння 
(1.6) геологічній моделі шаруватого середовища. 


У зв'язку з вивченням шаруватих структур введено |21, 103| клас Нумерова Ми СЯ п) ш 

2-5 А )х сів) як декартовий добуток множини 5, (р) обмежених і (майже всюди) 
локально інтегрованих функцій, що описують густину сасі(х, у) (х, у)єр мас в області 
рек? та сукупності Сз (кі?) Коб; О«а «1 функцій ух), хе К(? непере- 
рвних (за Гельдером) з їхніми похідними до К - го порядку включно з показниками ай, що 
описують контакти різнорідних шарів зі смуги П- (С, у): -о« ха; Й (х) «Л ) ,дел- 
обмежена стала. У працях |76, 96, 103, 104) розв'язки рівняння (1.9) розглянуто на класі 
мине, п) зі сталою густиною (відзначеною символом 1). Показано (103|, що на класі об- 


межених східчастих контурів обернене відображення поля в параметри середовища є розри- 
вним. Оскільки нелінійний оператор 


І Ср а 
Ах) - бен 1 





прямої задачі для контактної поверхні обмежений, неперервний і компактний лише на мно- 
жині Ми, п) функцій, неперервних (з показником а) на ЕФ разом зі своїми першими 
похідними, то і задача розв'язання рівняння (1.6) поставлена некоректно на множині 
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ми (І, п) неперервних функцій (це випливає ще Й з того, що кожну східчасту функцію 
можна подати як границю певної послідовності неперервних функцій). Доведено теорему 
єдиності розв'язку рівняння (1.6) на Мие (І, п). 

Доведення здійснено узагальненим на випадок необмежених зірчастих тіл методом, 
створеним зусиллями П. С. Новікова |42|, О. О. Заморєва |22), Л. М. Сретенського (541, 
І. Т. Тодорова, Д. П. Зідарова |671, М. М. Лаврентьєва |341 1 О. І. Прилєпка |48| для обмеже- 
них зірчастих тіл. Суть методу полягає в побудові й оцінці певного функціоналу, значення 
якого слугують ідентифікатором дійсності висловленої пропозиції. Функціонал будують з 
використанням фундаментальної інтегро - диференціальної тотожності Остроградського - 
Гріна, він відображає рівність зовнішніх потенціалів (або їх похідних) і властивості густини 
тіла, що не породжує зовнішнього потенціалу. Функціонал визначають на певній гармонічній 
функції, яку, у свою чергу, будують за визначеними даними граничної задачі Діріхле в 
об'єднанні областей, що зумовлюють згадані потенціали. 

Вибором методу доведення теорем єдиності розв'язку визначають і сам структурний 
клас Ми, п), до якого цей розв'язок належить. Проблему існування розв'язку задачі та 
питання стійкості розв'язків розкривають майже "дарма" завдяки спеціальному поданню 
лівої частини рівняння (1.6). Задавши її у вигляді двох складових, одна з яких (невідома фун- 
кція), як виявилося, домінує за нормою над іншою (інтегральним членом, залежним від неві- 
домої функції), автори отримали можливість побудувати таку послідовність компактних 
операторів, яка збігається до необмеженого оберненого до Й Га ; х) оператора, і за його допо- 
могою, у свою чергу, визначити точний розв'язок рівняння з множини допустимих. Таким 
чином було встановлено, що задачу відновлення контуру необмеженого шару за даними спо- 
стережень потенціалу поставлено коректно за Тихоновим на нумерівському класі 
Ми (1, п), тобто, що на класі Ми (ма, п) розв'язок задачі існує, єдиний і стійкий. До того 
ж, цей розв'язок залежить (як і в задачах для обмежених тіл) тільки від одного параметра сб і 
зовсім не залежить від параметра Й (середньої глибини до контакту за Нумеровим або його 
асимптоти за Страховим). 

На основі створеної теорії розроблено два чисельні регулярні способи розв'язування 
оберненої задачі за наявності перешкод |76|. Побудова одного з них заснована на тому, що 
значення поля, за якими визначають фрагмент контуру, задають на практично необмеженому 
інтервалі. В цій ситуації апроксимацію оператора прямої задачі Й Га 23 можна виконувати з 
високою точністю. Побудову другого здійснюють за умови, що значення поля задано на ко- 
роткому інтервалі і апроксимацію оператора прямої задачі можна здійснити з невисокою то- 
чністю. У першому випадку, за умови високоточної апроксимації оператора АС - 23) для чи- 
сельного пошуку каркасу наближеного розв'язку рівняння (1.6) конструюють регуляризую- 
чий алгоритм за способом М. М. Лаврентьєва |34|; у другому - для знаходження каркасу 
розв'язку використовують алгоритм, побудований в межах теорії регуляризації за Тихоновим 
165). У цьому алгоритмі вперше обтрунтовано вибір саме такого стабілізатора, який відпо- 
відає характерним властивостям розв'язків й забезпечує пошук розв'язків на множинах єди- 
ності, 1 рекомендовано певну нестандартну процедуру для точнішої порівняно з традиційни- 
ми апроксимації рівнянь Ейлера, що не потребує перемноження матриць високої розмірності. 
Ретельно вивчені питання скінченновимірної апроксимації оператора Й Га . х) і послідовних 


наближень оберненого оператора 4" (2 ; х). Особливу увага приділено залежності чисельних 


розв'язків від багатьох параметрів, які супроводжують дискретизацію задачі, 1 вказано спо- 
соби вибору таких значень цих параметрів, які сприяють одержанню оптимальних чисельних 
розв'язків. Доведено збіжність каркасів наближених розв'язків до точного розв'язку рів- 
няння (1.6). Детально з теорією оберненої задачі для контактної поверхні читач може при 
бажанні ознайомитись в дисертації Н. Н. Чорної |761. 

Переходячи безпосередньо до задачі з істотно обмеженим інтервалом, перш за все слід 
сказати, що поняття "істотно обмеженого інтервалу" запозичене з роботи |77| з цим уточ- 
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ненням: для визначення в рівномірній дискретній мережі точок х, з постійним кроком Ах 
деякого фрагменту с(х)х є Їсьсь | контактної поверхні сх) за виміряними на інтервалі ее 





еерец відомою точністю (Сх) - и; (2), «6 значеннями поля и(х) в деякому ітераційному 


процесі (|103| на основі рівняння (1.1) слід обчислити таке значення межі інтегрування а 
прямого оператора (1.1), за якого похибка його апроксимації не перевищить величини, крат- 


ної похибці поля ІєСолу| |єОса)| « 5,145, де є(х,а) - точність скінченновимірного проек- 





тування, «(х,п) - точність чисельного інтегрування (додатковий член формули чисельного 
інтегрування). Якщо поле и хх) задано на інтервалі ро так, що виконується вкладання 
ге -а-Дх /2,с,-а- Ййх /2| од то інтервал є майже необмеженим, в іншому разі - іс- 
тотно обмеженим. 

У першому випадку за співвідношенням (1.1) непогано відновлюється контакт, у ін- 
шому - надійно відновити контакт з допомогою ітерацій з ядром Пуассона неможливо і не- 
обхідно застосовувати інші, більш швидкозбіжні схеми обчислень?. 

Моделлю геологічної будови для вищезгаданої задачі є два зіркові відносно нескін- 
ченно віддаленої точки однорідні циліндричні шари різної щільності є, що залягають один 
над одним Й обмежені смугою 

Паіїкхб»-о«х «ой, «С(х)ЄВ,), (1.10) 


де й зіпі Сх), й" звиро(х) - глобальні екстремуми контакту СХ), а реіобі| і 


Й, з 8ир ен - постійні, описують множину контактів сх) зі смуги (1.10) при 
о«п «зх)є й, хе К (рис. 13). 





Рис. 1.3. Модель геологічної будови для визначення контакту способом Нумерова. 


Якщо в класі Страхова Й, п) |76) виділити на підставі співвідношень 
оС) « СО), о) ох є В, (11) 
оби «Сх «141, О« «Ак, (1.12) 
що обмежують "коливання" 710а ) контакту, деякий компактний клас Нумерова миє, п) 


то, як показано в |21, 103|, для контактних границь, описаних функціями з цього класу, голо- 
вну роль у гравітаційному притяганні підземних мас відіграє перша складова в лівій частині 
рівняння (1.1) , а друга - лише роль деякої малої поправки до неї. Наслідком цього є можли- 


вість подати на класі Ми, п) нелінійне рівняння (1.1) з точністю 0 29 щ Їмх)- пРр лі- 


2 Одну з таких схем для рівняння (1.1) якраз і наведено в роботі |771. 
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нійним інтегральним рівнянням Фредгольма 1-го роду з ядром Пуассона 


15 й | яр Ю 
нео сор ух) мс) с м(х)- Ах є КМ), алу) 


розв'язок сх), хе В) якого на класі Нумерова існує і єдиний при заданні правої частини на 





довгому профілі. Наша задача полягає у відновленні контакту при обмеженнях (1.10-1.12) за 
полем, заданим на "необмежених" со» і істотно обмежених сьо | множинах. Для відно- 


влення контакту при заданні поля на істотно обмежених множинах рекомендується еквіва- 
лентне рівнянню (1.13) лінеаризоване рівняння зі швидкоспадним ядром Шварца |211: 








РИ Це яа- 9) сема -и(), (114) 
УРОИРІЕ, (зе - З а - : умами 


розв'язок якого існує (як наслідок лінеаризації) 1 єдиний. Це не єдиний шлях визначення на- 
ближень 0 (ж, п): за належності шуканого розв'язку до компакту 0! (х)є ми, п) лінійне 
наближення контакту 0 (2, п) є сумою розв'язків 0 (ж, п) У "(х, п) 5 (ол) лінійних інтегра- 
льних рівнянь 


А Їзен зно 46 дія! 2л)в ць дове о ц(х) (1.15) 


Існує ще один варіант знаходження ненні С (2, п) через розв'язання рівняння 


Іов 
Сп) 25" (хоп) ейронза 


де 5 "(а п) визначено з першого рівняння системи (1.15). Таке наближення відрізняється від 





(845, (116) 
х) 


точного розв'язку 0 (х) рівняння (1.1) не більше за квадрат відхилення 0 (х)е Ох) - В, 


п епай,хеє КО), однак значно гірше за якістю від наближень (1.15) через наявність інтег- 


ралу з повільноспадним ядром Пуассона. Розв'язки рівнянь (1.13-1.16) залежать від апріорі 
невідомої середньої глибини до контакту П (що, до речі, є наслідком модельних зображень 
складних середовищ, а не принциповою характеристикою задачі). Її значення з геологічних 
міркувань не завжди відомі надійно, тому важливо розробити спосіб знаходження контакту, 
у якому було б усунуто подібну параметричну залежність. Цього можна досягнути шляхом 
введення нових узагальнених співвідношень для контакту 


п уявно СЕС 0 ат 


Йо чать, але 


Фо с3 
ще УРИХР ВИГРАЄ В  АЬ 2 Р 
п 2 (2к6(х))/ (5 - х) 

Подібна заміна справедлива й для рівнянь (1.15) 1 (1.16). Практична значимість цих виразів 
полягає в можливості отримувати значення глибин контакту в кожній фіксованій точці рів- 
номірної мережі в ході ітераційного поліпшення моделі, що надійніше, ніж підстановка апрі- 
орних значень ЛП. Співвідношення (1.13-14, 1.17-18) цілком працездатні 1 у випадку магнітоа- 
ктивних контактів, коли в їхніх правих частинах стоятимуть похідні напруженості магнітно- 
го потенціалу; контакт відновлюють із системи спряжених рівнянь типу (1.15) як суму зна- 
чень 
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и(ху-п) ш 57 (», в) 57 (», р); и(х, 1) ш 57 (», п)- 57(», п). 

Ця обставина відкриває перспективи для комплексної інтерпретації гравімагнітних 
полів і вимагає ретельнішого вивчення. 

Крім того, в дисертації (автореферат якої подано як |21|), наведено також відповідні 
формули для постановки контактної задачі в просторовому випадку, що теж є результатом 
досліджень за останні роки. Наведемо тут лише вигляд просторової моделі контакту 
(рис. 1.4) та основні вирази відповідної задачі. 








- є? 









іо чн! 
дріх,о0 якої 
чн 
Х» І п! ек 
ба 2, -''ПТЬн ОТ. 
Гл 

ою чно! то Ії, о чн! бовче 
чн! ЛЕН 

чно! Гл я і і чо 


Рис.1.4. Просторовий вигляд моделі контакту. 
Отже, з інтегрального рівняння контактної поверхні 0 (2, п) у постановці Нумерова 
ел) 


азатреу (о у, 1). 
Ме-с(му) 
після нескладних РАВ маємо для обох постановок Нумерова 1 Заморєва наступні вирази 


слі) зда У орг), ел) НС), 
я х) ї. зано о 
рай унуо ЛЕК зп) 
ЗЕ я (по у) | у) «Р 

ца, у, б.) є (5, п) - (2, у), 


справедливі з точністю до членів другого порядку малості за умови 





2лусб! (х, у) у)у-Уус і П 


зз00--00 








Я асат Ух, у, 2), 





««1, де 


НеЕН (2 п) чи Нас (х, у) , відповідно. Перше рівняння лінійне, а друге переходить в ліній- 
не, якщо звести його до формулювання Б.В. Нумерова й спростити вигляд. Подальше пере- 


ведення отриманих вище виразів на , просторові рейки" і чисельна регуляризація відповідних 
дискретних співвідношень становитиме основу подальших досліджень. 


1.3. АЛЬТЕРНАТИВНА ПОСТАНОВКА КОНТАКТНОЇ ЗАДАЧІ. У працях |78-81, 105| струк- 
турну задачу гравіметрії трактовано як задачу визначення за даними спостережень контуру 
необмеженого зіркового тіла, близького до шару постійної товщини. Основним рівнянням 
задачі в цій постановці є співвідношення (1.2), одержане Б. В. Нумеровим понад 70 років то- 
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му. Для відновлення границі ОЇ): у -с! (х), хек? збурювального тіла відомої сталої гус- 
тини с за значеннями похідних логарифмічного потенціалу, заданими на множині єдиності, 
або, що те ж саме, для розв'язування нелінійного інтегрального рівняння (1.2) запропоновано 
П105) такий спосіб. Розв'язок рівняння (1.2) відшукують у вигляді границі послідовності 
розв'язків певної послідовності лінійних інтегральних рівнянь першого роду, що апроксимує 
вихідне нелінійне рівняння. Встановлено, що оператори послідовності лінійних рівнянь є 
обмеженими, неперервними і компактними в просторах функцій, які найчастіше використо- 
вуються при розв'язуванні обернених задач потенціалу. Виконано спектральний аналіз опе- 
раторів, на основі якого вказано ефективні критерії належності вихідних даних 1 розв'язків 
відповідним їм області значень оператора і області визначення нормальних розв'язків. Дове- 
дено локальні теореми існування та стійкості розв'язків задачі визначення необмежених зір- 
частих областей, близьких до заданих. Показано, що задачу поставлено коректно тільки на 
компактних множинах банахових просторів вихідних даних та розв'язків. Внаслідок умовної 
її коректності для отримання ефективного розв'язку слід застосовувати спеціальні регулярні 
методи. В основі підходу закладено теорію згладжувального функціоналу над просторами 
даних і розв'язків. Детально вивчено його складові: функціонали нев'язки або функціонали 
типу нев'язки та стабілізатори |78|. Показано, що функціонали нев'язки чи типу нев'язки лі- 
нійного інтегрального рівняння першого роду з невід'ємним компактним оператором є стро- 
го випуклими, двічі неперервно диференційованими функціями. Множини точок мінімумів 
функціоналів над гільбертовим простором даних є множинами другої категорії по відношен- 
ню до множини всіх точних (нормальних) розв'язків інтегрального рівняння (множини єди- 
ності). Тому не кожна послідовність, яка мінімізує функціонали, збігається до множини єди- 
ності, через що задача пошуку мінімуму функціоналів над гільбертовим простором даних є 
некоректною. Недостатніми є Й додаткові вимоги, які широко використовують, стосовно мі- 
німальності норми розв'язку при мінімізації функціоналів нев'язки. Для орієнтування пошу- 
ку нормального розв'язку рівняння за допомогою мінімізації функціоналів нев'язки й типу 
нев'язки запропоновано спеціальні стабілізатори, визначені на множині єдиності. Їх конс- 
труюють за допомогою певних диференціальних операторів, власні функції яких тісно 
пов'язані з власними функціями операторів прямої відповідності в задачах визначення зірча- 
стих областей. Спеціальні стабілізатори -- це невід'ємні, сильно випуклі, двічі неперервно 
диференційовані функціонали, що не мають локальних мінімумів, окрім глобального. Впер- 
ше встановлено існування такого регуляризуючого оператора для лінійного інтегрального 
рівняння першого роду, який серед відомих його модифікацій має найкращу обумовленість. 
Кожній функції з гільбертового простору даних оператор ставить у відповідність функції Гі- 
льбертового простору розв'язків, що мають узагальнену похідну, яка інтегрується з квадра- 
том. У свою чергу, кожна з функцій, яку обчислюють за допомогою цього оператора, є екст- 
ремаллю згладжувального функціоналу, який є зваженою сумою функціоналу типу нев'язки 
з одним із запропонованих стабілізаторів, що забезпечують однозначну розв'язуваність вихі- 
дних задач. Глобальні регуляризуючі оператори визначення зірчастих областей, близьких до 
шару постійної товщини, будують у вигляді послідовності локальних регуляризуючих опера- 
торів для відповідних послідовностей лінійних інтегральних рівнянь першого роду, в правих 
частинах яких використовують згладжені вихідні дані. У свою чергу, локальні регуляризуючі 
оператори для лінійних інтегральних рівнянь конструюють у вигляді параметризованих пос- 
лідовностей, що мінімізують згладжувальні функціонали |79). 


1.4. УТОЧНЕННЯ КЛАСИЧНИХ НАБЛИЖЕНЬ КОНТАКТУ. В статті |99)| наведено стійкі 
способи уточнення апроксимації Андреєва |7| та інших, що, як там доведено, з неї виплива- 
Ж і ь 
ють. Нехай и(х) - вертикальна складова напруженості поля (х), а розв'язок 
С (ж, п)є Ми (І, п) є її аналітичним продовженням в бік тяжіючих мас на середню глибину 
й сх) шуканого контакту: 


22. 








що ї. (2-х) 4075) й р)- по п ей ІЛ 
созаодоо Й порно (вл) ма) ер Іще)-пиє. (119) 
Тлумачачи рівняння (1.1) через інтеграл Пуассона, маємо и, (х, п) С(х), м. (х, п) -С(), та 
17 П 
соопузид ос Ге те, ше, п), п ДАЄ (1.20) 


Апроксимацію ве п) розв'язку 0! (», п)є Ми мВ п) з точністю до квадрату величини 
п'є 2), де ЛЕ п' - ПО, визначаємо значенням каліброваного поля на рівнях у-0 і 
у -П (57): на першому рівні їх задають, а на іншому обчислюють як аналітично продовже- 
ну в бік тяжіючих мас гармонічну функцію за заданими граничними значеннями. 

За виразом (1.20) не можна строго визначити функцію их, П) через аналітичне продо- 
вження значення и(х) на рівень у-/»0, оскільки відповідна йому пряма у- / розтинає 
тяжіючі маси надвоє, одна з частин якої при С(х)« й виявляється нижче рівня у - ї, тому 


функція и(х, у) при у « п - негармонічна. Формально продовживши через інтеграл Пуассона 


о 


П 
исх,О) ча орто 


- 


(еле (121) 


"контакт" и(х) до рівня у - й, тлумачимо отримане значення и(х) як граничне значення по- 
ля мас, сконденсованих на рівні у - п. Значення сх, п) знаходимо як границю послідовності 
таз п) функцій при розв'язанні рівняння (1.21), генерованих запропонованим процесом 
121) послідовних наближень Лаврентьєва - Андреєва 


само - 2 є сема до біо 129 


Такий алгоритм обгрунтовує доведена в |99| теорема. 
Теорема Андреєва. Якщо розв'язок нелінійного інтегрального рівняння для контакту 





В! (х) є Ми (І, п), то його лінеаризоване наближення обчислюється з виразу 


п МС р у ЕМ азз 


РЕ й і шоб 
й від істинного ухиляється не більш, ніж на п "Ф (2) при П-»оо. 





З останнього виразу при п - 0,1,2 очевидні, як окремі випадки, відомі формули набли- 
жень Нумерова (44, Х'юза (119) і Страхова, відповідно: 








у -их 1 і й и 5 
СУ хв) ц(х) е Бопераеса «м, (1.24) 

9 х,п)еиіх 11 ей а и , . 
Ск п) ц(х) я пн (дає (1.25) 





оба 1Ї п 2 М зи " і 
ення Цеуодусуу рої" ин 


де ли(Є ) б и(25,0) - (20). Останні дві формули як 1 (1.23) при п 2 3 непридатні для обчислень 


--о0 


через наростання похибок з номером ітерації при обчисленні в околі точки 2 - х; уникають 
цього, підсумовуючи ядра інтегралів з різними знаками і зростаючими біномними коефіцієн- 
тами. Подібне не властиво для (1.24); перетворення формул з індексами п » 0 (підсумування 
елементарних дробів в ядрах перетворень) за схемою (1.23) при п-1,2,3 визначає додатні 
ядра перетворень, а при п - 4 - ні, тому в практиці обчислень наближень контакту 0 (ж, п) 
(1.24) є нульовим наближенням, а обмежитись варто однією з наступних формул (181: 
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о 2 4 
ебеьдію І збі знізиЧ(8 - х)? 2-х фа ри 
нн гориуснвча стук сту се з 
12 00 реобно н124пЧ2 - х)? 2902-х) (2-х) р 
сп) чих) ли(5 аа. 
Дучнір кому зх Дак ну Пан 
Значне уточнення зазначених наближень досягається через поправки до відомих спо- 
собів наближень Малкіна і Сенька, суть яких стисло окреслено в доведених в |99) теоремах, а 
детальніше виписано наприкінці статті |181). 





Теорема Малкіна. Корекція кожного з лінеаризованих наближень (1.23) С (а, п) на 
величину поправки 





АССоп)- Ем ее Ім(е, в) - м ие, (128) 


підвищує точність наближень С! ГР п)-с (ж, п) АС) до По (Є); в практиці обчис- 
лень достатньо подати поправку АС(х, П) у вигляді АС (х,п)- «суда І) а де цЄ,) - де- 


яке «середнє» значення поля на поверхні спостережень. 
Теорема Сенька. Уточнення кожного з наближень (4 (ж, п) поправкою, яка є 
розв'язком лінійного засо рівняння 1-го роду 


С А) 
дя (о п) бо 











підвищує точність наближень (29 (ж, п)- що и ПИ Кк -0,1,2,3.. до величини 
АС С уз тах|АС(Х). 


А. К. Маловичко |40) запропонував природне узагальнення ітерацій Б. В. Нумерова й 
А.К. Сенька для уточнення контакту, яке на підставі ітераційного процесу (1.22) нам вдалося 
удосконалити. Узагальнюючи поправки Андреєва 1 Сенька, запропоновано альтернативний 
спосіб, суть його |в обчисленні послідовних | наближень | й, (х)- 6, (х)-- ДАЄ, (х), 


ЛС, (х)- шт АС) за ітераційним процесом, названим в |21| процесом Нумерова - Мало- 











вичка: 
найк - к «01,23... о 
вер Кен з 
дні о де 
оо дя неї ее) 
о и КР тезою, пабу (1.29) 


Теорема | Маловичка. Якщо  С(х)є ми, П), то послідовні | наближення 
бабьп) є 0 (о п) А, (ж, п), ПО «АЙ, генеровані ітераційним процесом (1.29), збігають- 
ся при п-о»оо до граничної функції с (х), незалежної від параметра Й, зі швидкістю геомет- 
ричної прогресії. 

Доведення трунтується на тому, що процес (1.29) за індексом т є процесом Лавренть- 
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єва - Андреєва та збіжності 2, (ж, п) з відповідною мажорантою прогресії. 


Кожен спосіб по суті зводиться до розв'язання окремого лінійного інтегрального рів- 
няння Фредгольма 1-го роду з ядром Пуассона за допомогою чисельних методів |3|. Резуль- 
тати чисельної перевірки деяких з цих способів та порівняння їх ефективності щодо запро- 
понованих вище засобів розв'язання контактної задачі на короткому профілі запропоновано 
в останньому розділі дисертації (21) та наведено реферативно в останньому пункті даного 
розділу. 


1.5. МЕТОДИ РОЗВ'ЯЗАННЯ КОНТАКТНОЇ ЗАДАЧІ. Розв'язність рівняння (1.13) дослі- 
джував Б. А. Андрєєв 1 для його розв'язання запропонував спосіб послідовних наближень |71, 
а М. М. Лаврентьєв |34| обгрунтував його застосовність для розв'язку операторних рівнянь 
першого родуз, а згодом для рівнянь (1.13) 1 (1.17) в працях (103, 21) запропоновано ітера- 
ційні процеси Лаврентьєва- Андрєєва 


сонет (мен) ЄРОО-м 0 (130) 
лай нЕ-х) 


кана о еЕМенимеЙонх 0 аз 


Розв'язність рівнянь (1.14, 1.13) вперше досліджена А. В. Чорним |1031, а для обчислення їх 
розв'язків в |21| введено нові ітераційні процеси, названі процесами Лаврентьєва - Чорного 


ЧУТЬ; Дао) МО ВУЗИ РИ ТР УР РЕА ТЕР ТР О; 











ДП 2 





орд) СРО -| С; б | СЕМ чно), сх) дп о 0,2, (133) 
о с3 0) 
МОГ РУ ПАРТІ 1) до (3 и(5)45 
С ана 
Дослідження збіжності цих процесів елементарне (індукцією за збіжності мажорую- 
чих часткових сум) внаслідок еквівалентного подання розв'язків через наступні наближення: 


)- "Їемеї) У (1-е 91)! сова - ха, (1.34) 





га ян га У (сво) ") ско - х Мао, п є 0, (135) 


о кед 

Твердження про збіжність цих набіннані водночас, згідно |411, є конструктивними теорема- 
ми існування розв'язків рівнянь (1.13) 1 (1.17); аналогічні співвідношення, чисельно складні- 
ші, виведено для узагальнених рівнянь (1.14) 1 (1.18). Вивчення швидкості збіжності різних 
ітерацій? для рівнянь (1.34) 1 (1.35) виявило, що послідовні наближення Лаврентьєва - Андрє- 
єва С ЩЕ); і Лаврентьєва - Чорного С о рівномірно збігаються до граничної функції 
і (х) є Ми!-о) (1, П), яка задовільняє кожному з двох еквівалентних рівнянь (1.13) 1 (1.17) зі 
швидкістю геометричної прогресії, оціненої , відповідно, величинами 
ОК Р РЕ СЛАР ТО РЕ 

Ма-- 2) ач 2) 
а модулі неперервності обернених операторів цих рівнянь оцінено як 











З Ці обставини відбито у назві запропонованих в |21| ітераційних процесів. 
4 Відсутність таких даних заважала оцінці їх чисельної ефективності і можливості порівняння різних алгоритмів з метою 
вибору оптимального для розв'язуваної задачі. 
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рдчаное ау нео 
де М(п) « 2пах Сх); Мп) «є 2 пах (у 21); Мп) є 2 пах |мх) чн м 2), х є Ва м(х 21) 


і моди) - аналітично продовжені значення »(х) і м) як гармонічних функцій з рівня у - 0 














на рівень у - 21. Схожий результат справедливий і для узагальнених рівнянь контакту. 
Порівняння наведених модулів неперервності вказує обернену пропорційність швид- 
кості збіжності ітерацій глибині ЛП і майже однаковий її порядок для обох процесів за невели- 
кої (усього вдвічі) переваги у швидкості збіжності процесу (1.32) над процесом (1.30). В обох 
процесах швидкість спадає з наростанням глибини Й й індексу ітерації п, але практично іте- 
раційний процес Лаврентьєва - Чорного має перевагу над "андрєєвським": завдяки швидкому 
згасанню по горизонталі ядер перетворень (1.32) для відновлення наближень 5 (п) з визна- 
ченою точністю можна обмежитись значеннями поля и(х) на істотно коротших профілях, 


ніж за використання будь-яких схем з ядром Пуассона. Однак, останніми за певних умов мо- 
жна користуватись, істотно прискоривши їх збіжність. 


1.6. ПОСТАНОВКА І РОЗВ'ЯЗНІСТЬ ЗАДАЧІ В КОМПЛЕКСНІЙ ОБЛАСТІ. Вище ми довели 
коректність постановки контактної задачі, яка описується інтегральним рівнянням (1.1), що 


залежить єдино лише від густини аномалієтвірного шару, на заданому класі Мис П) 
віє (р)х С (ав) в дійсній площині. Згадані здобутки, як показано в праці |98|, цілком 
можна поширити і на комплексну область. Дійсно, ще в роботі |103| така задача (1.6) була 
редукована до розв'язання на класі Ми, п) нелінійного інтегрального рівняння типу 
Урисона, яке в комплексній площині 2 - х - іу має такий вигляд |104): 





Уа есмнои 1 зоб 15 
д5 з 5 1.36 
2 о оо і це) і у 
де и) пе С(2)/ л/с - калібрована напруга С(2) поля логарифмічного потенціалу притягання 
С (2) за виразом (1.3), 7 - гравітаційна стала; с - густина активного шару, 
С(2)- Ще) Р 90(ху) Сх.) 5о0тХхяк їСС), 5з-бз ЩЕ) - відповідно фіксована 
да 2 дх бу 


точка контура ОР і точка, що його пробігає, а 5, «ях іс! (х), 5ївб-іс (2) - спряжені з ними 
точки. Перетворимо це рівняння до зручного для порівняння з поданням (1.9) вигляду. Оскі- 
льки, як доведено в |97|, 





(1.37) 
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а зіл, є р, 
іл, іє В), 


др, 
де Р)" - область над кривою дР),,а Д) - під дІ),, то, очевидно, 


р 


2лі бр, 5-2 





Якщо ж серед множини обмежених знайдеться крива ОД, що має асимптоту Й ліво- 
руч й праворуч, то, взявши у співвідношенні (1.36) 5 - 5, - Ци іс! (х) -2ій, отримаємо 
хі-»о0 








(оче бе Аз зцє)ч й. 


аз з ца) або ж | Ї з 

др» лі др» 
Порівнюючи останній вираз з рівнянням Страхова (1.9), дійдемо висновку, що це рівняння є 
частинним випадком від наведеного (1.36), в якому до того ж не враховується (як і у попе- 
редньому рівнянні (1.8)) притягання шаром постійної товщини ЛП. Окрім того, подання (1.9) 


від рівняння (1.36) суттєво відрізняється не тільки тим, що воно описує значно вужчу мно- 


5-2 
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жину можливих контурів збурювальних шарів, а й, головне, тим, що воно залежить від па- 
раметра й, якому в теорії |60, 71-74) надано визначальну роль. 
Розв'язування рівняння (1.36) розпочнемо з обтрунтування співвідношення (1.37). 


Лемма 1.1. Якщо контур ОРД: у-с (х) х є В) належить до класу неперервно дифе- 
Рі гер, 
-1і сер 





ренційованих (за Гельдером) функцій Ми! ай, п), то -- а 
Ж да 


Доведення проведемо для випадку, коли хє Д)". Оскільки функція виду / (з) 


з (5 - 2) є аналітичною в незамкненій області Р" |2) з виколотою точкою г -х-кіу, то на 


основі відомої теореми Коші |35| інтеграл, що нас цікавить, можна подати в такому вигляді 
(рис. 1.5): 


1 аз а аз аз Я з аз аз аз аз аз 
пнів пере рерао. р | 











«пт де ДР, і дО; -- частини контуру Фр 
між точками (х,, С(х, )), (0,2(0)) і 
(х,С(х,); С, - розтин між точками 
2 « (0,0) і (0,20)); С, і С, - коло ра- 
 діуса є -|5- | Й концентричне до 
нього півколо радіуса К-|5 -2| 
центром у точці 2 «(0,0); С, і С, 
дуги секторів, кути розхилу яких від- 
- повідно дорівнюють: 

















Кк 
"а Ф, з агсів о) о агоїв ---о- т 
Х,, 
Ф, ВР С О ЛИНЕ ОРОРНН ЗЛ -»0. (1.38) 
Хе 
ха | зо 








Криволінійні інтеграли за дугами ОД), і-12, очевидно, існують тільки тоді, коли 





контур с (х) належить до класу Ми, п). У о випадку б 


1 4 1 НО ЕВ С ак 
| | зе ана. не ія. 0 а) 


дрі др 








цілком визначають фіксованими числами х,, і х,,аїї границю при х, -2-00, х, -»00 з ура- 
хуванням відношень (1.38) знаходять так: 

2 о Й Р 
1 ) іє шо я змМ4Шф Ї Кіе?аф і і і 


-- т 
2-50 . ріф . оід 
пі о ЗВО зма Се К-е 





Випадок, коли 2 є Р)" , розглядають аналогічно. 
На основі доведеної леми 1.1 нелінійне інтегральне рівняння (1.36) подамо у вигляді 
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ЯН потав ПЕР фам) (1.39) 


зручнішому для подальшого дослідження. Далі в статті |98| наводиться наступна теорема. 

Теорема 1.1. Якщо контур ОР: у-б (х), хе В належить до класу ми о, П), 
тобто задовольняє вимогам Ма (2Ї М ЇМ 2) сят, т, М «є со, то оператор прямої за- 
дачі 


Ас) аву но ГП 


для 





: 1 
є обмеженим, неперервним і компактним оператором на св! ) і 
Доведення обмеженості оператора випливає з очевидної оцінки 


УЗА РИСА Р РА Р А С 
паро Я По55-2 


ор 


«3М. 














Для доведення його неперервності розглянемо різницю 
1 ж ж -Ж Й І У кер, 
АС я пуг)- А(32)- 2 Да зо Пав, 
лізь 5 -ї длі оо 
де 5 «Єчідзіпезвчіт, др: ус СС) т), хе В; лох) -- достатньо мала варіація фун- 
кції 0 (х) у тому сенсі, що функція лох) для довільного числа є » 0 задовольняє вимогам 








Іо, з пахрб|«є; 0 лсо)-п(оо)-0. 
. 5 2 
Оскільки варіація лох) мала, то Е. г і І М і росі Я г.) ков 
5-2 з-21чц(я-) 8-7 5-2 Ц5-2 


. Я . РЕЯ І а) по 
і тому Аа ту 2)- Аг) є іх) 2лі р з 2лі Пежо 





з Ї ПФ .., Звідки неважко 


вивести рівність АС А 1) б АС - 2) « Зітоо) « 7є/2, яка засвідчує неперервність оператора 


АС г) в ск) через те, що число є » 0 довільне. 

Компактність (повна неперервність) оператора АС ,2) означає, що він переводить 
довільну обмежену множину, зокрема клас Ми (1, п), функцій в компактну (за метрикою 
с(к?)) множину функцій це) АС ,2)є У(). Раніше ми переконалися, що оператор 
Й Га :2) обмежений на довільній обмеженій множині функцій, тому для 0 є Миа) (1, п) ма- 
тимемо ім(2)| Й «у (2 е- Цією нерівністю встановлено рівномірну обмеженість множини 
УА). Доведемо її одностайну неперервність. Нехай 2 що | «є для довільного числа є»0). 


Тоді різниці 





Мед моде) оно ДП 5 - 


з урахуванням розгортання 


1 е 1 1 1 зокре п ронавй б 1 Рая А 
5-2, 5-2, 1-(2, 72,)/(5-2,) 5-7» зом доня а 


можна надати вигляду 
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(2, - 2; з 5-5 з 
лі 2(8 - 2/5 -2,)" 


Проте функція 0 (х) неперервно диференційована на КО) і, до того ж, 








ху о ха| «2, з 25|, че- 


рез це норма різниці набирає вигляду: 


Ім на о «де 





а Ї зо а ««(дт 2). 
2|2лі 25(5 - 1) 

Одержана нерівність означає, що множина У(А) функцій ц(2) є одностайно непере- 
рвною. За відомою теоремою Арцела |291, множина (4) є компактною, а оператор АС ,2) 
- цілком неперервним (компактним) у метриці простору неперервних функцій. 

Визначальна роль доведеної теореми 1.1 для розв'язування рівняння (1.39) полягає в 
тому, що вона орієнтує напрям пошуку розв'язків за їх характеризацією. На конкретних еле- 
ментах того чи іншого класу функцій можна показати, що неперервна залежність розв'язку 
інтегрального рівняння (1.39) від його правої частини відсутня не тільки в просторах обме- 


жених (східчастих) функцій, а й у просторах неперервних С (кі?) або ж інтегрованих (абсо- 
лютно чи з квадратом) І9(к?) Кк - 1,2 функцій. Для того, щоб у цьому переконатися, до- 
сить розглянути значення оператора АС ,2) на певній послідовності неперервних функцій, 
що збігаються до тієї східчастої, яку розглянуто у праці (103|. Оператор Й Га : 1), як встанов- 
лено теоремою І, неперервний принаймні на класі Ми (а, п), і щонайменше тільки на цьо- 
му класі потрібно відшукувати розв'язок рівняння (1.39) для вхідних даних з просторів 
св?) та Іов 2) К - 1,2, анаступні теореми обгрунтовують цю процедуру. 

Теорема 1.2 (єдиності). Розв'язок нелінійного інтегрального рівняння (1.39), якщо він 
існує, єдиний на множині Ми (ад (1, п) неперервно диференційованих функцій. 

Дійсно, однорідне нелінійне рівняння 


лето За ГР оо (140) 


ге ЗД реж 





на множині Ми, п) має тільки тривіальний розв'язок. У цьому нас переконує суперечли- 








ва нерівність |, оо - З ен а « - |з Я Ге. « Е |з З» виведена на основі 
С 


леми 1.1. Отже, однорідне рівняння (1.40) має тільки тривіальний розв'язок, а рівняння 


ор ор 


(1.39), якщо й має розв'язок, то він однозначно визначається на класі Ми аа, п). З огляду 


на те, що довільну функцію 0 (х) З МНОЖИНИ миЧЗ1, п) можна подати у вигляді 


х 
С (х) шо (0) . Ї С 1977 , ця множина є рівномірно обмеженою і водночас неперервною, тобто є 
0 
компактною множиною функцій. 
Теорема 1.3 (існування). Процес послідовних наближень 


і (со) 
зані) и 5ні) о У Ї 5 5 





0 - (0 
су ща), по 02, 8-5 са) 0 (141) 
арії) У па; 
збігається на компактній множині Ми маг (1, п) до (нормального) розв'язку нелінійного інте- 
грального рівняння (1.39) зі швидкістю геометричної прогресії. 
Для доведення цієї теореми достатньо показати збіжність ітераційного процесу (1.41) 1 


визначити його швидкість. Твердження стосовно того, що гранична функція 0 (х) послідов- 


ності У ох) буде розв'язком нелінійного інтегрального рівняння (1.39), випливатиме з тео- 
реми єдиності розв'язку. Задля з'ясування збіжності послідовних наближень розглянемо різ- 
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ницю 
МОВО 


М РА ; , І ( п 
(зн) хо ?) (зо ще нови ах? | ето з ах) 








Для спрощення викладок позначимо 5(ані) - 50 -Й кон -С о п), п -- 0,1,2,..., 1, скорис- 
тавшись тим, що прирости пі? за метрикою простору неперервних функцій св ) набагато 


менші самих наближень 4 (п) контуру 5, матимемо дз щ дз) з іаті?), пвеО12.. 


І І КР РН А ; 
попе 5) за Лак (59) - 2) 51) -2 5-2 180- и 


У цих позначеннях попередня рівність виглядає так: 


























і (а) І (а-і) І (ау очі) 
пан о | рано | рю 
Злі злраоеа 2лі зр -І 2лі ОНИ 
п-1 п-1 
ва 1 Ї З, а (п) 1 Ї 4 дз) я 
2лі ні 5) -2 2лі нс 870 -з 
2 
1 по іо" іо" іі і 
- 1 - Лазо па я 
Злі 57 СЛИВ Р А і 
(п-1 п-1 п-1 2 п-1 
Ж 1 Ї 4 і іт! р іт! | і (ах? іт) І Ї а ах (по?) 
Злі аа) 5 (1-1) Й 5 (1-1) шир я (п-ї) Й Злі дат з (пі) -7 


Звідси з точністю до величин другого порядку щодо приростів пФ отримуємо спів- 




















відношення 7)" бе Ї По п по р з якого встановлюємо нерів- 
ТІ орні) | 5 чн 5337 -7 
ність 
г М 2 (п-1) | с-(п-і) д фо 
р, є Но, Дн - б щи у чем 
Ця нерівність слугує основою для наступного ланцюжка 
З Р КА Р З Р Р 

















який, у свою чергу, засвідчує, що оператор (2; 2) задачі для контактної поверхні є стиска- 
ючим. Наближення 0 (с), п 0,1,2,..., збігаються до розв'язку 0 (х) рівняння (1.39) зі швид- 


кістю геометричної прогресії зі знаменником 2". 





Зауваження. Оператори Ви, С "ан ог Ї ; су 2 дз? І) послідовності 
лі -7 


Ви, С (а) яка, очевидно, апроксимує при 509) 5 519) б ц) обернений (нелінійний) опера- 
тор А41(С52) рівняння для контактної поверхні, є обмеженими, неперервними і компактними 
операторами на компакті У (4). Виявити ці властивості можна за допомогою засобів, які було 
застосовано під час досліджень оператора Й Га ; 1) прямої задачі. До того ж, кожен з операто- 
рів Ви, С (2) п 0,1,2,..., проектує на нумерівську множину Ми (ад, п). Переконуємось в 
цьому, переходячи у зображенні оператора до його параметричного подання, а саме 


(п) Яа" 
"й Кз ен З 
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Наявність у зображенні (інтегрованої) похідної від наближення 0 є) контуру 


4292) 
4є 


2), Є є В), є ознакою функції, що належить до множини Ми"(1, П). 

Отже, в перебігу доведення теореми існування розв'язку знайдено таку послідовність 
Ви, С (а) компактних операторів Ви, (а) 2) п -0,1,2,..., яка обмежену множину УА) 
функцій и(2) однозначно проектує на компакт Ми еУ(1, п) і прямує за метрикою простору 
неперервних функцій ск) на певному елементі и) множини (4) до нелінійного обер- 
неного оператора 4" (2 ; 1) рівняння для контактної поверхні. 

Теорема 1.4 (стійкості). Нехай функції с, (х), хе В), і-12, що описують границі 
розподілу однорідних шаруватих середовищ, належать до класу Ми, п). Якщо за дос- 
татньо малого числа є» 0 поля и, (2), г є В), які породжуються відповідно контурами 


2.(), підкоряються вимозі и (2)- и, (2) 





| с 26, тобто їх вважають близькими, то самі гра- 
ниці також незначно відрізняються одна від одної в сенсі виконання умови 


ІС (х)- 6» (2 «де. 


Для доведення теореми зо різницю 











4) Р ПИ З РИ 
2 (зів ни Й Рой ш(е)чни, (2). 
І «д-аід і о ПР А р 
Оскільки ре ІТ хол з Ар 5 а М т Я Бора дз? 








цу авее рн ой очна 


дз? С, Ге аз) аз) 1 
| 5) за Ї 50) оре -510-6СЇ 


грі) о 2л грі) тя 


то | (х)- є, (Їх « Зм (2)- и (2, «25, 


Вищезгадані теореми окреслюють умови коректної розв'язності оберненої задачі ло- 
гарифмічного потенціалу для контактної поверхні в комплексній площині. Задачу поставле- 





но коректно на компактній множині ми, п) неперервних разом з першими похідними 
функцій, бо на класі Ми аа, п) розв'язок задачі існує, єдиний 1 стійкий. Особливість прове- 
деного аналізу - редукція задачі до розв'язання нелінійного інтегрального рівняння (1.39) 
типу рівняння Урисона. Ліву частину отриманого рівняння, залежного від невідомої функції, 
подано у вигляді двох складових, одна з яких - невідома функція - домінує (за нормою) над 
іншою, що є інтегральним оператором над невідомою функцією. Вигода такого подання (на 
зразок лінійного інтегрального рівняння другого роду) чітко вбачається в напрямі дослі- 
дження рівняння. По-перше, подання виявилося доступним для вивчення умов єдиності, іс- 
нування та стійкості розв'язків винятково елементарними засобами. По-друге, виглядом рів- 
няння визначився і спосіб його точного розв'язання. А з огляду на те, що точний розв'язок 
відшукують за допомогою певного процесу послідовних наближень, сконструйованого, у 
свою чергу на основі вихідного подання рівняння, це подання стає вирішальним і в практич- 
ній реалізації розв'язку задачі в умовах, коли праву частину рівняння задано з перешкодами 
на дискретній множині точок. 


1.7. ЧИСЕЛЬНЕ МОДЕЛЮВАННЯ ЗАДАЧІ. Основні результати з чисельної перевірки ви- 
щезгаданих алгоритмів здійснені в дисертації, де окреслено загальну стратегію вирішення 
некоректних задач, здійснено алгебраїзацію пропонованих алгоритмів і наведено, у рамках 
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регуляризації за Тихоновим, методику чисельного розв'язання лінійного інтегрального рів- 
няння Фредгольма першого роду, до вирішення якого зведено всі досліджувані процеси, по- 
дано результати чисельного тестування на комп'ютері та порівняльного аналізу тестових 
прикладів, отриманих за різними алгоритмами. Описано програму обчислень на комп'ютері 
наближених розв'язків поставленої задачі на рівномірній мережі вузлів. Згідно |21|, характе- 
ристики прямих операторів АС - х) задачі 


АСе і к(оєУ(ЄЗАЄ чиїх), сс, (1.42) 


дозволяють, в разі точного задання його та правої частини, без проблем будувати таку послі- 
довність компактних операторів, яка збігається на певному елементі класу Мида) (1, р) до 
точного розв'язку. Через нестійкість оберненого оператора внесено корективи в цю страте- 
гію. Скінченновимірну апроксимацію системи (1.42) здійснено методом проекцій, що дозво- 
ляє зменшити розмірність системи, оскільки вхідна матриця симетрична 1 додатно визначена. 
В підсумку для операторів рівнянь (1.30) 1 (1.32) відповідно маємо з точністю похибок апро- 
ксимації: 











Р А Р . г Й 
АВ, х)є Р, А"(В,х) Г онузар (2 тає 
Ап, х) є Р, АВ, х)- | ев те Це,тУе, (143) 
ь пр и з 1 ТЕ -А 
де похибки є (ха, . - Ге «У аг, |є (за, ба з іі ад. 


З цих оцінок знаходять межі інтегрування а, 1 за формулою Гауса 7-го ступеня точ- 
ності (12 вузлів) здійснюють чисельне інтегрування. По завершенню дискретизації оператора 


Ап г ча в) отримуємо відносно шуканого вектора розв'язку 2 (х, )", такі системи лінійних 


па! 


алгебраїчних рівнянь: 


2п-1 Я Ах а т СС») Й т аб 
ех) 1222 печу: Є, 2 ) 


2п-2 Ах о 1 т 


2 о й Да л(є, - А (о Ум) (144) 


В ході регуляризації за Тихоновим ця система редукується до вирішення рівняння Ей- 





лера 
п-1 пост 
аУ є(ху я уУ(х ;К Же. зб ук, (х ба у, )- У», (ан, ун, (х,) (1.45) 
ка уві ізі 
за умов застосування квадратур Гауса і Гах Ма Ума ун к( А ) де залишковий член 
ізі 


(р,)- і ах - У» Л(о) має смисл похибки наближеного інтегрування і дорівнює 


2 


адлчі (пі)? 
2 дл) (дл) 





1 
ій пп), -1«п7 «1. Для переходу від інтеграла Ї / («мах до інтег- 
-1 


КЗ 


Б п 
рала з довільними межами використано відоме співвідношення Ї дОдах « УА (ох, - В) 
Я ке! 


для обчислень кратних інтегралів з пропущеними вузлами 
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Й (хе кає з ау У ав (В нах б кру я ГОВ зок ру ГВ ах бу 


-ї ізі 








-1(8- ах, - ук де азбов вет, р дн для по2 - мірного куба, 


де, К(р)- 19/35. В ході Уо на інтервалі - 11 отримано КОД, ) 211-109, 


5 
Розв'язуючи цю систему для ряду значень є2-Сб., С-У 10, знаходимо шуканий 
ую у ряду ук 
ї«0 


озв'язок 02, (х, ), |. Оптимізація розв'язків із знаходженням ш/,, для модельних прикладів 
а Кк /Ккеї І 


ор: 
не проводилась і тут ще є простір для вдосконалення. Модельне поле, яке фігурує в правій 
частині всіх алгоритмів, отримане розв'язанням прямої задачі, тобто обчислення полям, ге- 


нерованого модельним контактом типу Іванова: 


с(х)- ма) і ко) (146) 
й, іє | сокс5| 


Значення отриманого поля, збурені 5Уо-ним шумом, служать вхідними даними для вирішення 
оберненої задачі. Під час тестування вивчено чисельну збіжність другого рівняння системи 
(1.32) (рис. 1.6) та поведінку їх інтегральних ядер. Останні нерівнозначні за точністю і швид- 
кодією Працездатність алгоритмів (1.30 - 1.32) проілюстровано на модельних прикладах з 
параметрами: довжина профілю / є ЇсоезТЕ І- ід- 2 км; крок Ах «0,1 км, перепад густини 
с-03г/см'; коливання контакту - 1624 (х)« -1,4 км; апріорна глибина контакту 
й - 1,2.км; амплітудні величини 


азеї; па15,Ю, «2,А-10,фаейекеазь-і, Ах «04133, АС « 0,066, б «107, 


-(-1.А-ніДх).. Чисельний алгоритм такий: 





17. введення вхідних даних /,п,К., а, Р,Ф,П, ,К , значень поля є, та розрахунок додаткових па- 
раметрів чисельної схеми т,т і У ра о 


22. вирішення прямої задачі - визначення вно поля, генерованого модельним ко- 


ха та (ач (Є, зх.) 


нтактом , за формулою и(х,)е 216 оно у; - коефіцієн- 
ізі 


ти квадратури Гауса, його екстраполяція за формулою 


ба ок - Оп 
(хе дхУ фо У; ау "| ЕР т(є.) 5 


ізі узі 





т 


т(є )- У С1)/ у (віп МАЕ віп (2) 612, - іп З Є, сов(2) 22); 
уні 
32. зашумлення отриманого поля й 02) псевдовипадковою похибкою з рівномірним розпо- 


ділом и, (х,)- У, Ар д; /т), Д;г- су ом г/59 - 0.5) для використання вектора 
із1 


ізі 


го 2тчі - . 
ї(х, 1 як початкового наближення оберненої задачі 


ізі 
«Ума УС зем 1)-н1, т - шкфуп-н1. 
ізі ізі 


472. розрахунок спряженої правої частини Кр через розклад значень и(х | у ряд Фур'єї 
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наступний синтез 5 (х,)є ув о "АП я зптодьмо оз 
ї«і узі їз! Езі 
ко Злі . Олі 
сі У 1-со8 -8іп : 
руч Жтяі Жт і 


59. формування моделі початкового наближення: задання апріорної інформації Л,є,б та 
звуження області еквівалентності введенням обмежень 0 ЖЕ. 1) та са на шукані 
розв'язки 





У 5) - У 3"(а)є Уа уві 


ша! ізі 
672. розрахунок границь інтегрування а,,, ініціалізація параметрів ітерацій с, йег, І, т, п, К, 


а, Б, Ка, 5 та обчислення значень інтегральних ядер К, (х, і 1) за квадратурами Гауса 


к.(х,п)е тра Ро р еттмі?й ей із 5 й 0.5)Ду 
ізі узі 











кн ен 
е'-е"! ях, -т,). 
Кк, Со н)з ЗМ поту ерня б ль ше зо 28, 


2. обчислення наближень 0 х 3 із відповідного ітераційного процесу: 


5"(х,)є У (і. АХ КС Д8 (юн 82 Ди, екс 


Кепдчі 9-1 
570. 23 . «Ук х С, 5 (х,)- 5, Де, КД 
Кепочі чі 
«(хе У5"()ч 5"(х,) 
ізпочі 
87. регуляризація отриманого розв'язку: обчислення параметра регуляризації а за принци- 
пом нев'язки та визначення регуляризованого розв'язку з розрахунком середньоквадратичної 





похибки с; з | Ус оп) и -п -1) ; при о - 5| , 5 Ємаш Вертаємось на крок 79, інакше 


йзпо-і 


фіксуємо в файлі результати обчислень 5) бе У; бо (х, У, г збо йек. 


а 

і ? С,» 
Тестування здійснено так: для заданих точних розв'язків 0, (2) іядер К, оо) розра- 

ховують аналітично або чисельно праві частини и, (х) , потім вносять в них випадкові похиб- 


ки (нормальний розподіл), а точні ядра замінюють наближеними К, (ж, 2) тієї ж гладкості з 





с» 
точністю є « Ї ЇК З (х)- К («Р ах .З цими даними вирішено два модельні приклади: 


с 


1. Модельне рівняння Ї к(оЄУ (2 з ц(х), с кс, точний розв'язок 


2 
Ця ання точне ядро к(хл)-е, при цьому Їсо»є, Ї-|-1.414), 
2 


п 29, Лх «0.1, Їсузеь Ї - І- 1,1, порядок регуляризації - нульовий (гладкі функції контакту, 
правої частини, їх майже однакова ефективна ширина) 


шах 
2. Все те саме, але точне ядро К (х, п)- 4/ле"" , що означає внесену похибку оператора 
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є я. 0.12 (тобто є,., я 12.59), За такої похибки складно знайти шукану функцію 2(х). 






ь х У х)«П- хЇ Мат 
З - ко ето 
4/ле"" 








Рис. 1.6. Теоретичне поле за модельним контактом (а) та тестові приклади (6). 

Ітерації закінчують по досягненні заданої середньоквадратичної нев'язки чи по вико- 
нанні заданої кількості кроків ітерації. Похибки отриманих розв'язків не перевищують зада- 
них неточностей поля. В дисертації (21| здійснено порівняння розроблених в роботі 1 відомих 
класичних алгоритмів Чорної, Нумерова, Маловичка, яке виявило більшу ефективність алго- 
ритмів Лаврентьєва - Андреєва, Лаврентьєва - Чорного порівняно з класичними схемами, зу- 
мовлену кращим ступенем обумовленості систем алгебричних рівнянь. Чисельні експериме- 
нти виявили перевагу алгоритмів відновлення контакту на короткому профілі над класични- 
ми способами за точністю та обсягом обчислень (рис. 1.7 -1.10). Деталі експерименту можна 
віднайти в дисертації. 
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Рис. 1.7. Відновлення модельного контакту за алгоритмом (1.32). 
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Рис. 1.8. Результати розв'язання прямої (а) і оберненої (6) задачі для алгоритмів Нумерова, 
Маловичка ії Н. Чорної. 
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Рис. 1.9. Модель контакту (а) і результат його обчислень (б) за алгоритмом (1.30). 
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Рис. 1.10. Відновлення контакту за узагальненим алгоритмом Чорного (1.33): 
а) - заточними значеннями поля; б) - за 3Уо-збуреними значеннями поля. 


Висновки. Дослідження представляють перш за все науковий інтерес, але алгоритми 
можна застосувати у складі автоматизованих систем для моделювання 1 експрес-аналізу гра- 
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віметричних даних на маловивчених територіях нескладної будови. Точність обчислень кон- 
такту залежить не лише від похибок поля 0, але і від довжини інтервалу його задання; за о0д- 
накової точності обчислень у випадку задання значень у(х) на короткому інтервалі для від- 
новлення їх за допомогою операторів зі швидко спадаючими ядрами потрібний менший об- 
сяг інформації, ніж у випадку задання поля на нескінченних профілях. Відтак, отримано те- 
оретичне узагальнення і нове вирішення оберненої задачі потенціалу для контактної границі 
в лінеаризованій постановці. Головні наукові і практичні досягнення такі: 

19, Описано математичну модель задачі відновлення контакту на класі Нумерова за 
вертикальними похідними логарифмічного потенціалу у вигляді розв'язків лінеаризованих 
інтегральних рівнянь Фредгольма | - го роду. В їх узагальнених модифікаціях усунута зале- 
жність розв'язків від апріорі невідомого параметра Й - середньої глибини до контакту. До- 
ведено існування 1 єдиність розв'язків на компактному класі єдиності мибе1, р). Модель 
задачі поширюється на магнітоактивні контакти. 

29. Доведено, що відомі алгоритми, які базуються на виразах з інтегральним ядром 
Пуассона, не забезпечують прийнятної точності відновлення контакту за полем, заданим на 
коротких профілях, які властиві практиці спостережень. Для випадку задання значень поля 
на істотно короткому профілі для визначення контактної границі г(х) запропоновано кілька 
алгоритмів, що базуються на інтегральних рівняннях із швидкозбіжним ядром типу Шварца. 
На основі рівнянь з ядрами Пуассона і Шварца введено два різних типи задач для відновлен- 
ня контакту (з параметром Й та без нього). Вивчено структуру класу єдиності зазначених 
рівнянь. Доведено умовну коректність контактної задачі. 

32. Для розв'язання рівнянь кожного типу задач запропоновано збіжні ітераційні про- 
цеси Лаврентьєва - Андреєва (1.30) 1 Лаврентьєва - Чорного (1.32) та їх узагальнені аналоги 
(1.31) 1 (1.33). Практичною перевагою ітерацій (1.31), (1.33) є швидша збіжність та можли- 
вість задання значень контакту на короткому профілі. Збіжність усіх ітерацій обгрунтовано 
за новим способом, що базується на еквівалентному зображенні підінтегральних виразів в 
ітераційних процесах через узагальнені функції. 

492, В ході дослідження вищезгаданих ітерацій поліпшено ітераційні способи уточнен- 
ня нумерівського наближення контакту, пропоновані свого часу Андреєвим, Малкіним 1 Се- 
ньком, введено новий спосіб Нумерова - Маловичка і обгрунтовано стійкість уточнень. Їх 
можна використати для аналітичних апроксимацій поля в автоматизованих системах інтерп- 
ретації. 

59. Для чисельного вирішення задачі розроблено регуляризуючі алгоритми на базі 
концепції регуляризації Тихонова. Модернізовано знаходження оптимального параметра ре- 
гуляризації а/,, мінімізовано похибки округлень. Створено пакет програм 1 перевірено на те- 


стових задачах. Порівняння чисельних розв'язків довело перевагу алгоритму (1.32) над алго- 
ритмом (1.30). Досягнута точність розв'язку є, я 107 при похибках оператора є, 10" і спо- 


стережень 0 2107" прийнятна для застосувань. 
Однак, ця теорія ще чекає своєї практичної перевірки 1 не позбавлена деяких недолі- 
ків, які дають перспективу подальшого вдосконалення, а саме, на нашу думку: 

1. поширення математичного апарату на тлумачення даних в площинному варіанті; 

2. залежність поведінки контакту (точності рахунку) від різних густин і глибин до конта- 
кту - вивчення впливу низької контрастності вхідних даних на роботу методу та ін- 
ших впливів: великі шуми, багатошаровість, еквівалентність, оцінка достовірності, 
статистична обробка; 

3. розробка діалогово-графічного інтерфейсу для промислового застосування програм та 
поширення задачі на просторовий випадок в рамках класу Нумерова єдиності 
розв'язків і заданих обмежень на функції поля й контакту (в дійсній 1 комплексній об- 
ластях); 

4. вивчення практичної залежності мінімальної довжини профілю від глибин контакту -- 
нині це лише змодельовано як вплив похибок на поведінку розв'язку в оцінці загаль- 
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ної похибки; 

5. вивчення залежності розв'язків задачі від складної форми контактної границі та мак- 
симальної глибинності методу за різного рівня вхідних шумів на реальних прикладах 
поля; 

6. перевірка швидкісних характеристик і схеми екстраполяції на розроблених алгорит- 
мах, чисельне порівняння прискорених і не прискорених алгоритмів розв'язання обе- 
рненої задачі для контактної поверхні та тлумачення результатів; 

7. вивчення можливості відновлення кількох контактних поверхонь одночасно; 

8. відображення контакту в інших полях, перспективи комплексування та застосування 
на складніших моделях середовища і полях більшої розмірності, чисельні проблеми 
методу; 

9. вивчення можливості застосування розробленого методу інтерпретації даних структу- 
рної гравіметрії для визначення положення покришок нафтогазових колекторів у ком- 
плексі з даними сейсморозвідки нафтогазоносної території, та для супроводу остан- 
нього методу при дорозвідці і моніторингу діючих родовищ вуглеводнів; 

10. вивчення можливостей методу відновлення контактів для завдань інженерної геології 
та моніторингу рівня техногенного навантаження в умовах міських агломерацій. 


2. ПРО ВІДНОВЛЕННЯ ГРАНИЦІ ОБМЕЖЕНОГО ЗІРКОВОГО ОДНОРІДНОГО ТІЛА, 
БЛИЗЬКОГО ДО КРУГА 


Вступ |89). Задача визначення форми тіла за його потенціалом поставлена давно в 
зв'язку з дослідженнями фігури Землі, хоча математично вона була чітко сформульована 1 
наближено розв'язана лише на початку 30-х років минулого сторіччя Л. Ліхтенштейном (|38| 
в припущенні, що потенціал шуканого тіла близький у сенсі деякої функціональної метрики 
до зовнішнього потенціалу заданого тіла. Тому згодом вона іменувалася задачею для тіла, 
близького до заданого, а її розв'язок - - як локальний чи "в малому". Незабаром, наприкінці 
30-х років, П. С. Новіков (42) довів теорему єдиності глобального розв'язку задачі, а Л. М. 
Сретенський (53| - - локальну теорему існування в припущенні, що заданим тілом є однорід- 
на куля. І тільки в 1943 р. А. М. Тихонов |62| встановив на підставі теоретико-множинних 
концепцій умови стійкого розв'язку задачі як безпосередній наслідок теореми єдиності. Ци- 
ми роботами був відкритий новий клас так званих обернених задач потенціалу, що визначив 
найбільш перспективний напрямок теорії обчислень. Правда, поки це було усвідомлено, 
пройшло довгих 20 років. Саме тоді М. М. Лаврент'єв |34| сформулював нові умови корект- 
ності за А. М. Тихоновим задач математичної фізики з центральною вимогою однозначного 
розв'язування задачі на певній множині функцій. Це відкриття докорінно змінило існуючі 
уявлення про те, що будь-яка задача, яка не задовільняє вимогам коректності за Ж. Адама- 
ром, не відповідає ніякій реальній фізичній постановці. Дослідженнями в різних галузях 
знань за останні 30-40 років встановлено, що множина процесів і явищ природи, математич- 
ний опис яких дається в рамках коректності за А. М. Тихоновим, по потужності, якщо не бі- 
льша, то, принаймні, не менша в порівнянні з множиною процесів 1 явищ, що описуються в 
рамках класичної коректності за Ж. Адамаром. Суть справи тут зовсім не в адекватності ма- 
тематичного опису того чи того реального процесу чи явища, а в способі опису їхніх власти- 
востей: або за безпосередньо вимірюваними їхніми характеристиками, або за опосередкова- 
ними їхніми проявами. Опис за опосередкованими характеристиками, як правило, приводить 
до некоректних задач. У свою чергу некоректні задачі найчастіше зводяться до інтегральних 
рівнянь із цілком неперервними (компактними) операторами. 

Для деяких класів лінійних і нелінійних компактних операторів і певних множин у ві- 
дповідних функціональних просторах М. М. Лаврент'єву вдалося побудувати функції корек- 
тності та за їх допомогою отримати ефективні розв'язки сформульованих задач. Наступний 
крок у розвитку теорії розв'язування некоректних задач зробили В. К. Іванов | 251 А. М. Ти- 
хонов (631. В. К. Іванов, прагнучи усунути ускладнення, пов'язані з відсутністю розв'язку 
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(лінійного) інтегрального рівняння першого роду з неточно заданою правою частиною, за- 
пропонував шукати його квазірозв'язок, тобто елемент, який мінімізує функціонал нев'язки 
рівняння на опуклому компакті, що містить точний розв'язок. А А. М. Тихонов для обчис- 
лення наближеного розв'язку некоректної задачі побудував регуляризуючий оператор за до- 
помогою мінімізації певного функціонала, який задається на множині можливих розв'язків. 
Роботи М. М. Лаврент'єва |34, 36, В. К. Іванова |24, 25| 1 А. М. Тихонова |63, 65| дали мо- 
гутній поштовх бурхливому розвитку теорії і методів розв'язування некоректних задач, а та- 
кож їх застосуванню у різних галузях знань, у тому числі й у геофізиці. Чільне місце в зага- 
льному потоці публікацій зайняли роботи М. А. Алексідзе (1-5|, А. В. Бакушинського |10), П. 
І. Балка (111), В. Б. Гласко ЦП 31, І. І. Данилюка та М. І. Михалюка |17|, В. М. Ісакова |27, 28, 
Д. В. Капанадзе |30), А. С. Леонова |66), В. П. Маслова |41|, С. М. Оганесяна (45, 46), А. І. 
Прилепка |47, 49), В. Г. Романова |50), В. І. Старостенка (551, В. Н. Страхова |38, 611, В 
Чередніченка |73-75|, Н. Н. Чорної |76, 77, А. В. Чорного |92-96), А. В. Цирульського |71, 
72| та ін. І проте всі ці надзвичайно важливі досягнення знаменують собою тільки початок 
великого наукового пошуку, котрий повинен завершитися створенням загальної теорії. До- 
тепер залишаються невирішеними питання коректного розв'язання за А. М. Тихоновим неко- 
ректних задач. Якщо М. М. Лаврент'єв (34) пропонує не торкатися зовсім проблеми існуван- 
ня (глобальних) розв'язків обернених задач потенціалу, а вважати їх заздалегідь відомими 
для певного класу даних, то В. П. Маслов |41| орієнтує на те, що одержання чисельного 
розв'язку некоректної задачі з заданою правою частиною методом регуляризації само по собі 
засвідчувало б існування розв'язку, а розбіжність регуляризованого процесу означала б не 
існування цього розв'язку. Іншими словами, за В. П. Масловим існування чи відсутність гло- 
бального розв'язку некоректної задачі пов'язується з реалізацією певного регулярного алго- 
ритму. Однак, останній критерій інколи викликає сумніви, тому що реалізація методу регу- 
ляризації не завжди буває адекватною. Так, у геофізичних обчисленнях метод регуляризації 
використовувався зазвичай формально, без врахування особливостей операторів рівнянь 1 
характеристик їхніх правих частин. Він застосовувався в найпростішому варіанті, яким його 
пропагували класики, у вигляді побудови мінімізуючої послідовності функціоналів нев'язки з 
додатковою вимогою мінімуму норми розв'язку і не завжди приводив до позитивного ре- 
зультату. Тому в загальній теорії розв'язування некоректних задач питання побудови адеква- 
тних регуляризаторів стають одними з першорядних. Як і раніше залишаються значні труд- 
нощі в питаннях єдиності та існування локальних розв'язків обернених задач потенціалу при 
відновленні різних класів тіл навіть за постійної їхньої щільності в зв'язку з винятковою 
складністю відповідних інтегро-диференціальних рівнянь, до яких редукуються ці задачі. 
Деякі з перерахованих проблем визначили коло тих питань, у рамках яких нижче розгляда- 
ється обернена задача логарифмічного потенціалу для обмеженого тіла, близького до круга, 
за відомої постійної щільності без істотних обмежень на розмірність арифметичних просто- 
рів, що апроксимують вихідні дані і розв'язки. 


2. Математична модель задачі. Нехай на нескінченній прямій 5 « (-с0,- оо), уз- 
довж якої направимо вісь (Ох, декартової прямокутної системи координат Ох,х,, заданий 
логарифмічний потенціал притягання 


Об; РА зер» 


2 
2 2 
хер| ем ду зуду,, ххх), 

ізі 
тіла (9, с) постійної щільності с, де 7 -- гравітаційна стала. За відомими значеннями по- 
тенціалу і заданої щільності потрібно визначити границю 2 9 області 9? , що притягає, при 
додатковій умові, коли задані на 2) значення потенціалу (/(х; Ф, с) мало відрізняються (у 
сенсі визначеної метрики) від значень потенціалу (/, (х; З ,, С) круга 9,(2, р) радіуса 2 з 


39 


центром у точці 2 - (8 |, 2 .), заповненого масою тієї ж постійної щільності с, з якою розпо- 
ділена маса по області 9 . Вісь Ох, направимо в бік збурюючого тіла так, щоб вона прохо- 
дила через центр круга, тобто через точку 4 - (0,п), р « й. Тоді при малому відхиленні зов- 
нішніх потенціалів М(х; ?, с) 1 (х; 9 ,, с) областей ? і 9,, шукане тіло ? буде мало 
відрізнятися як за своїми розмірами, так і за своїм розташуванням від апроксимуючого круга 
Ф о. Тому область З буде зірчастою відносно центру круга 5, 

Наведену задачу зручно розглядати в полярній системі координат з полюсом у центрі 
С е(0, й) круга 9, 1 з відліком кутів проти часової стрілки від полярного радіуса, що збіга- 
ється з віссю Ох, . Позначаючи через о(с/), 04 с 227, однозначну функцію, що описує кон- 
тур зіркового тіла ? , а через С (а) з р(а)- р -- відхилення контуру тіла З? від круга 9, 
одержимо вираз для потенціалу (/(х; 9, с) з врахуванням того, що область Ф близька до 
круга 7, 








и(х)а ен У ду зи, (Ф, г)- би(Ф, г); (2.1) 
ао)«С(х;9Ф,с)/2/о, би(ф, к)- Їзіа. Да), РФ, гаа, 
(0) 
чо» 1 зд 
ик(Фф,г)з Га Па зд в лр'п-, 
о 7 5 - 2зг сов(д - ФУ)" г 





Фа, 2,ріфу)- сову 7 й п К"(у, р, г)- 


не зо еХа газ), ба) 
4 К (ши, р, г) Кк (м, р, г) | 
нраваут іа б'ЄЙе гово), а) 
4 К (ші, р, г) Кк (м, р, г) | 
2 - 
о й 2 аксіє "Ца)зіп у - 
і ву, різ З аХо -поову) 
В (уз рог) 
В наведеній формулі (2.1) доданок и, (Ф,г) є потенціалом круга 9 ,, а би(ф,г) -- по- 


га) 2рс(а)н са) З2рс(а)-с"(а) 
2 




















тенціалом "надлишкового" тіла ДФ - 9 19, стосовно цього круга. Справедлива 
Теорема 2.1. Якщо область 9 , що розглядається, мало відхиляється від круга 9, в 
сенсі виконання умови 
пах | (а) | «піп Во - фу, р, г), 
а а 





де | "ау р(а)-р, тіК(а-, р, г) - тіл? -2ргсов(й-ф)ч «Н01-у), 


узрі" «1, то підінтегральна функція 
р(а) 
1 


Ф(а, Є ,руф,т)- | то РУ 
Р Ме -2зг сова - ФУ) ЧІ 








що визначає варіацію потенціалу 





би) з Ї Прой 


зо се У| 


уз і Ф(а, (ад), рі р, где, 


хр ег5пФ, х, -Ппаегсо5Фф, О«ФХ27, 
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розгортається в абсолютно збіжний ряд 

уба і 2 ЗОдБИ 1 

го, є А оо бр'"  Ка-ф,р,т) 
Доведення. Подамо підінтегральну функцію Ф(а, С,рФ, й що обумовлює варіацію 








себто 


би(ф,т) в зображенні (2.1), в нових змінних (/8, т), котрі зі старими змінними (и/, у), 


уз ради, 





у | « Ї, зв'язані наступними співвідношеннями 


шт РЕ АС 
К(у, р, г) К(у, р, г) др 
У результаті одержимо 


В? 2 р? 2 р? 
ФВ, т; р, Юертв-", совра і Я отв - ва рейа рвеор з 
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Кк Кк І 
- ог 2в|наз 2тгсо8ф-- 238 - рами В- ее До ка р Д 


Ін тсо8 8! 
де К - К(м, р, г). А так як (105| 


ЩІ тов Вк) - ра ту" б сор, Кано 





т8іп зб 
агсі 1 и К т|«1, 2.1 а) 
ет Зі й | | 


в" 
со8 кВ а ( Рон (м, р, г) б 
(Е-1! др 
те попереднє зображення Ф(), г; р, К) можна переписати у вигляді 


2" б о дія 1 б з 252 
Фа, С, р Фуг)Е-рстК-32-|пКоер--пК|- 5 2- цу п К |, 
(а, Є» руфьг)- рат во, Р нота дові 
яке збігається з тим, яке потрібно було довести, якщо за згодою прийняти для п - 2 вираз 
2 2: 2-2 
є |рАнроюВ -5 о ря бокі. 
2 ор 4 р""др Ор" 


На підставі теореми 2.1 нелінійне стосовно функції 0 - р(0), 0 227, інтегральне 


рівняння (2.1) перетворюється в наступне нелінійне інтегральне рівняння для функції 
сас(а), Очажоя, 





п К(у, р, г), 
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з віп ду агсія 


"с (а)зт ши | 
Ву, р.гуя с(а(р- сову) 

причому функція /(и, 2; р, г) відновлюється за функцією Ф (г., С; р; Ф, г) з вихідного вира- 
зу (2.1). 

Але якщо замість відхилення С((/) контуру тіла З від кола 29, підставити в рів- 
няння (2.1) наближення 2, «2. (а), п- 0, 1,2,..., цього відхилення, то наведене перетворен- 
ня рівняння (2.1) можна розглядати як лінійне інтегральне рівняння першого роду для визна- 
чення наближень. Таким чином, нелінійне рівняння (2.2) для функції 2 (и) перетворюється на 
лінійне рівняння (2.3) для функції 2, - 2, (2) , якщо покласти 


2 
Ї ка, фу р, г)Є, (адаа є ур б, рег), 
0 
поа1,2,,ю0; С24(2)-0; (2.3) 
1 
УФ, оо різ Ге гучяріпг, 


у, Ф) З МФ бут) рот); 
Беручи до уваги все сказане, уточнимо постановку задачі визначення контуру однорі- 
дного фінітного джерела за формою, близького до круга. Потрібно знайти однозначну функ- 
цію с(а), 0 ж 227, що характеризує відхилення границі однорідного фінітного джерела 


2 


(9, с) від кола 29, ху -(х, - Я) е р' і визначає контур 29: 2(а) - р-С(а), за зна- 


ченнями логарифмічного потенціалу області 9 на нескінченній прямій поза областю за 
умови, що потенціал мало відрізняється (у визначеному сенсі) від потенціалу круга за відо- 
мої щільності с розподілу по області 9 гравітуючих мас. У даній постановці відшукання 
функції с(а) здійснюється у вигляді границі послідовності (2, (//)) , що генерується проце- 
сом наближень (2.3). Надалі співвідношення для відшукання наближень 2, (с) будемо коро- 
тко називати інтегральним рівнянням для визначення контуру тіла, близького до круга. 

2.2. Властивості оператора прямої відповідності. Оператором А: Ж -» хх прямої ві- 


дповідності простору розв'язків )ж у простір даних хх будемо називати інтегральний опера- 
тор АС - у, заданий рівністю 


Їна, рФус(ааа, -м(Ф), С є, уєхх, (2.4) 


де для стислості позначено К(а., Ф) - К(а, Ф5; р, г). З очевидної його лінійності Й елементар- 
них обчислень легко встановлюється, що оператор А є обмеженим, неперервним і компакт- 
ним оператором із С(0, 2л) у С(0, 2л), із С(0, 2) у 8(0, 27), із 8/9 (0, 2л) у 8(0, 27) 
, Кк 1,2. 105) 

Функція К(а, Ф) -- ядро оператора (2.4) -- за природних обмежень на параметри р, 
г 1 функцію (а) « г, зв'язаних з розмірами шуканого тіла, по сукупності «кутових» змінних 
неперервна (при фіксованих 2; і г) в основному квадраті 0 а, Ф.«27 і симетрична відно- 
сно цих змінних у тому сенсі, що К(о., Ф) - К(Ф, а). 


- ж У . 
Тому спряжений оператор А , що задається рівністю 
21 


Г Ка, фу УФ 4р з та), 


0 
- - . - ж . 
збігається в гільбертовому просторі за своєю дією з оператором 4, тобто 4 - А 1, отже, 
2 
оператор 4 у ЗУ7(0, 2) буде самоспряженим. А це в свою чергу означає, що вся сукуп- 
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ність (Л,) власних значень 42, самоспряженого (компактного) оператора А -- власна час- 
тина його спектру с(А) -- складається з дійсних чисел, до того ж множина |Л,) власних 
чисел і відповідна їй множина (х,) власних функцій х,, Ах, « Лх,, не більше ніж зліченна 


в З? (0, 2л). Для відшукання власних функцій і власних значень самоспряженого компакт- 


ного оператора задачі (2.3) злегка перетворюємо його ядро |105| та отримаємо наступний ро- 
зклад ядра 
т уї 
к(а,Ф,У)З- -ПРЯ 2. (З ка со8 КФ -к58іп Ка 58їп КФ) 
Кеї 
в абсолютно збіжний ряд за умови у - ог. «1. Звідси випливає, що послідовності 


п п 


у у у 
-О2лійк, лу, Пр лайно, СКУ У 
1 1 п п 
1 1 1. 1 1. 
сто ро 080, - , ЗІ,  - СО8ЛИ, ІПН 5 
мало «Я 4 4П Я 


утворюють повну систему власних чисел і відповідних їм власних функцій неперервного си- 
метричного ядра К(а, ФУ), Оса, Ф2ОЗл, у«1, самоспряженого компактного оператора 
А рівняння (2.4). 

А це в свою чергу означає, що спектр с(Д) оператора А дискретний і цілком визна- 
чається зліченним числом значень 


п 


с(Ау з -дліп г, 0, 
п 


Далі зауважимо, що області визначення операторів 4 і А" суть Ф(А) з 8 (4 Ф'В(А) і 
Ф(Аузв(А)ФФ(А), де 'Ф(А), ЯР(А ) -- області значень цих операторів, а 
9 (Ау «(2 АС) «0, С є 8У(0,2л)) та З (А) «(01 А"(2)-0, С « З8Р(0,2л)) -- їх нуль- 
множини. Причому в гільбертовому просторі 8 (0,2л): 8 (Ау- З (А), Ф(Ау «я (А). 
Що ж стосується лінійного стосовно наближення 2, (с) інтегрального рівняння задачі (2.3), 


то воно є рівнянням Фредгольма першого роду. 
З наведеного випливає, що самоспряжений компактний оператор А задачі (2.3) вста- 


новлює взаємно однозначну відповідність тільки між підпросторами ЯР (А") і ЯР (А) просто- 
ру 87 (0, 27). Тому оператор А має обернений тільки на "Р (А). У зв'язку з цим необхідно 
прояснити два теоретико-множинних аспекти. Перше, що необхідно тримати в полі зору при 
розгляді розщеплень для 9 (4) й 9 (4 ) угільбертовому просторі, - - це відсутність ефекти- 
вних критеріїв приналежності правої частини рівняння задачі області значень його операто- 
ра. Справді, якщо г є 2 (4 ), то для всіх // є 897 (0, 27) будемо мати 


0 з (472, 2) « (2, АС) «(2 у). 
Звідси видно, що нуль-множина 2 (А ) ортогональна області значень "Р (А) оператора А 
задачі (2.3). У зв'язку з цим не можна вказати ефективний критерій приналежності функції 
у(а) з ЗР7(0, 27) -- правої частини рівняння 44 зу -- до множини "Я? (А), тому що в да- 
ному випадку він має вигляд нескінченного числа співвідношень (типу скалярного добутку). 
Залишається єдина можливість «розпізнавання» елементів у(//) множини  ЯР(А) з 
ЗР (0, 27) за допомогою різних процедур «згладжування», що зводяться власне до побудо- 
ви середніх функцій за Стєкловим. Друге, що заслуговує на увагу в зв'язку в розщепленнями 
ФА) 1 Ф(А"), так це те, що для розпізнавання області визначення "Р (4) нормальних 
розв'язків рівняння задачі (2.3), необхідно ввести додаткові критерії, тісно зв'язані з внутрі- 
шньою його структурою. Дійсно, оскільки для будь-якого 2 є 8 (А) й усіх ує 8 (0, Ол) 
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має місце рівність 
ж 
0 -(А2,у) «(2, А у), 
то нуль-множина 2 (А) ортогональна області значень "Р (4 ) спряженого оператора А" за- 
М у ж 
дачі (2.3). А так як ортонормована система власних функцій оператора А (- А ) виявляється 
- - - ж - - - 

повною не тільки в підпросторі ЯР (4 ), але і у всьому просторі 83 (0, 2л), то, як і в попере- 
дньому випадку, за ознакою ортогональності підпросторів 5 (А) і ЯР(А") не можливо виді- 


лити однозначно елементи структури ЯР (4/7). Для розпізнавання підпросторів у розщепленні 
доводиться залучати додаткові характеристичні ознаки. Так за відмітні особливості множини 
Я ( А") можна обрати рівномірну в середньому обмеженість його підмножин і їх одностайну 
в середньому безперервність (чи еквівалентні їм критерії). Для матеріалізації цих ознак мно- 
жини ЯР ( А") визначимо для довільного елемента у є 8 (0, 27) й відповідного йому елеме- 


нта (а) є ЯР (7) за досить малого приросту |Ло/ « є аргументу а різницю 


аз ла) - а) з Ма Ла, фу у) - Ка» фу у (руд - 


- лУ у"(-а, зіппа - 8, созпа)Ла -- ода) з ллас"(а)я о(даї ), (2.5) 
па 


де а, ф, -- коефіцієнти Фур'є функції у(Ф). Переконаємося, перш за все, що 


С(а) є ЗР (0, 27). Безпосередні обчислення приводять до рівності 
2 и 
|с"а)| зя"У "(ар ЯР). 
пе! 


А оскільки за нерівністю Коші-Буняковського 








то 





хЇУ| 


М 
і, отже, С"(а) є 8 (0, 27). 

На підставі отриманого результату і безпосередньо з подання (2.5) видно, що ознаку 
одностайної неперервності множини функцій 2 є ЯР (А") можна замінити еквівалентною їй 
ознакою рівномірної обмеженості перших похідних функцій цієї множини. Що ж стосується 
ознаки рівномірної обмеженості, то для функцій множини ЯР (А") вона виглядає, як неважко 
переконатися, таким чином 








р 
Іе |єяуї ковку нар «ліні; ар. зу (ол). 
У зв'язку зцим множину ЯР 5 ( А"), яку згодом зба називати областю визначення норма- 
льних розв язків задачі (2.3), можна описати у вигляді 


ФА) а ау 29 є РО, 27) (ач 21) -с"(а); || «с, кої, 2, 0 де 


Са) а Са), С з шана), с) с ладен (іп А)? 
Г2 





сен - рі). 








2.3. Проблеми існування і стійкості розв'язків |89). Ми вже переконалися, що опе- 
ратором А прямої дії між підпросторами "Р (А) і ЯР (А) гільбертового простору 3 (0, 27.) 


. . . - є 
встановлюється взаємно однозначна відповідність, за якої для кожного елемента 0 є Я (А ) 
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в ЯР (А) знайдеться єдиний елемент у, для якого виконується рівність 42 - у. Тому повинна 
існувати взаємно однозначна відповідність елементів "Ф? (4) елементам ЯР (47), тобто по- 
винне існувати з обернене відображення | 4": ЯР (А)-Ф (А"), неперервне |на парі 
і ЯР (А), ЯР (А) ) підпросторів 1, взагалі кажучи, розривне (можливо і необмежене) на усьому 
ЗР(0, 2л). До того ж область значень ЯР (4) оберненого відображення повинна бути за- 
мкнутою (по відповідній метриці, що не збігається з метрикою вихідного простору) множи- 
ною першої категорії, оскільки ця множина породжується в результаті неперервного відо- 
браження множини першої категорії 1» (4) за допомогою обмеженого на цій множині опера- 
тора 4", тобто повинна збігатися з ЯР (4). Одночасно, при побудові оберненого відобра- 


ження щоразу виникають певні труднощі, оскільки в рамках гільбертового простору, приро- 
дним чином зв'язаного з задачею, що розглядається, відсутні ефективні критерії приналеж- 
ності елементів підпросторам "Р (А) і ЯР (47). Характеризацію підпросторів ЯР (А) і (А) 
доводиться здійснювати на підставі додаткової інформації про вихідні дані 1 розв'язок задачі, 
що забезпечує умови однозначного розв'язання відповідних рівнянь. Тому питання про існу- 
вання розв'язку оберненої задачі потенціалу типу (2.1) доцільно вивчати не на всьому вихід- 
ному гільбертовому просторі, а на більш тонкій його підструктурі. На підставі цього цікаво 


відзначити ще одну точку зору на підпростір ЯР (А") як на множину можливих розв'язків за- 

дачі. Розглянемо підпростір 

- ЛВ 

ФА) аа са): Є є ФР0, дл; | 2" «| С9| 
каб 


функцій з 8 (0, 27), що мають (узагальнені) похідні аж до т-го порядку включно, кожна з 


яких належить до (87 (0, 27). Очевидно, що при К « т виконується нерівність | 0 ? «є ? 
і тому 

УР 1 і і оо! а Р. М ра о і 
тобто 


Фа)-Пе ОСА"). 


Отже, підпростір ЯР (А7) (також як і підпростір "Р (А)) є зліченно-нормованим. Виходить, що 


зліченно-нормовані простори є природними в постановках обернених задач потенціалу, оскі- 
льки ці задачі, як правило, зводяться до функціональних рівнянь з компактними (лінійними 
чи нелінійними) операторами, які в тому чи іншому банаховому просторі мають потужні 
нуль-множини. Однак, цей погляд навряд чи доцільно використовувати на практиці, оскільки 


вихідні дані задачі задаються в метриках просторів С чи 637, інколи в С"? чи Ф 0? (коли 
похибка мала разом з першою її похідною). 

Умови існування розв'язку задачі визначення контуру обмеженого зіркового тіла роз- 
криває 

Теорема 2.2 (існування). Якщо збурююча область 9 , що породжує зовнішній потен- 
ціал и(Ф, г), і круг 9,(0, р) радіуса р, р«Ппетіпг, що індукує потенціал и, (Ф, г) , мало 
відрізняються один від одного в сенсі виконання нерівності 

тах (а) З тіп К(а-ф,р,г), Оха,фхадОл, 


де Са) ре р(а) -- однозначна функція, що описує конфігурацію границі 2 З області 9, 
а функція 
/ 
Ва-ф, р,г)« (р' -Зргсова-фяко)?, 


то послідовність наближень (С (а)), яка генерується процесом 
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Га ру Є, (ада - ме узифог)- ач ву ума) 


п 0,1,2,..., ж; С,(а) «0, 
збігається (за метрикою С(0, 2л)) до функції С(а), тобто до розв'язку р(ас) нелінійного 
інтегрального рівняння 








2л р(а) 1 
Г4о | о здз зи(ф, г) 
о 5 5 -дзгсов(а -фУ ЧИ" 


за умови, що кожне з наближень с (а) належить до компактної множини 


ФРС? є ЗРО, лу а калу з Са; |6 |, с о, 1, 


























І? 
(0 - 2 й и 
с""-яі2сч(2п А) а | 2с чи 
і що виконується нерівність 
1 ме 
г 4 ФА де ше р''тах|с(а)| і узир. 
чи вано 








Доведення. Для встановлення збіжності послідовності (2, (а)) складемо різницю 


Гока, фу уУЄ, (а) - б, (аа а ПОЛО» бо УЗ ГУ ба» Уа, 


і спробуємо оцінити знизу ліву 1 зверху праву частини отриманої рівності. Нагадаємо, що тут 
К(а, Ф, У) - - т К(Ф, р, г), 


А СО А Р А - 
2 2 У)- Щ | 
Л(9, с, у) 2, йо р" дріЙ" др" Вир,г) 











Введемо позначення 
с, з тах о (02), б» з тіп 2 (02), бо з тах (а). 


і на підставі теореми 2.1 запишемо систему очевидних нерівностей 


ФІГ) є реа ФД)- ГО), 


К(у,, р, г) 
14) -ФИДУ Є ра, ад є ГОД) ФСУ), 
К(у, р, т) 


з якої визначаємо 


(АФ) є рій (а) - Є, (ан АФ, 
К(у, р, г) 


де АФОФИД)-ФИ), дГеГС)- ГО Ф)- ріс х ГО, тар. С, », 


Ф(2) «Ф(а, С р; Ф, т). Інтегруючи цю нерівність, одержуємо після нескладних, але громізд- 


ких викладок, що з одного боку 
2л 


Гока, фу мЇС, (а) - 6, (аа |2 


Пе РА нео о! пара 











22 





п" ЗФпч Ї- 
пе? 
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п-1 


2 4лРрП- ш) п ДХ ж 





а з другого 


їз Гогов, би ДУГИ Со МО є 


























йо зайоьЇ Це РА (ря со е г у ра « 
«Ж2лр0чш) о) бар Іс 








Спеціально підкреслимо, що як нижня ліворуч, так і верхня праворуч оцінки виразу 


Г ока, фу мс, (а) - Є, «(аа 








отримані за припущення, що не тільки сама гранична функція / (и), але і всі її наближення 
С, (0) мають, принаймні, перші похідні, що інтегруються, тобто, що всі вони є елементами 
простору Соболева Ф 7(0, 27). Якщо ж 2, (а) є Ф (2(0, 27), але є елементом 83 (0, 2л), 
то оцінка знизу виписаного виразу дорівнює нулю 1 обернення рівняння на довільних елеме- 
нтах (8 (0, 27.) не може бути побудоване. 

Порівнюючи тепер між собою отримані вище оцінки для елементів послідовності 
С (с) з Ф Ф(0, 27), приходимо до такого співвідношення 


| аю аб іо - 4 ер ой 








С? 











де 
11 а 
и яр щш СОП5Ї 
"Аз 23 
1-у 





Звідси переконуємось, що послідовність (2, (а/)) збігається в собі при 4 «1, оскільки ряд 


Усла) «4 (а)) 


за цієї умови збігається, тому що 
2 


























ре ПН пан НЕ па АВ ор ро ДИ 
і, виходить 
о Р г а 
У ана)-слау| «ра об 
пеВ С пед -4 














Але 4 «1, якщо виконуються умови теореми. Дійсно, умови теореми виконуються, якщо 
1 
у 
а - НО. мл, (2.6) 
1-у 


тобто, якщо шу«1І-у. А це співвідношення завжди виконується для функцій класу 








аз НУ 
1-у 





ФО )с Ф (0, 2л), тому що воно еквівалентне до вимоги -т«С(а)«т, тепП-р»д, 
У випадку, коли виконується вимога (2.6), варіація параметра ш не перевищує 3, 
тобто амплітуда відхилення області Ф? від круга Ф,(0, р) не перевершує 3 2. 


Множина "Р " (А7) замкнута за нормою 
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неї «бернієтй) 


і тому послідовність (|. (а/)) збігається до деякої функції 2(а/) з Я "(А). 








Переконаємося, що функція р(0/) - р-- С(а/) є шуканим розв'язком нелінійного інтег- 
рального рівняння 








2я Р(а) і 
Гча Ї по зд з и(Ф, г). 
0 0 8" -25гсо5( - Фу 2 


Для цього в обох частинах рівності 
2л 


Г ока» фу у), (аа: з Фу, г)- м, (ру г)- ЦО» бу» Ма 


о 
перейдемо до границі за індексом п-» оо і одержимо вираз 


2л 2л 
1 
Гос зда зизи, - |Про С, ума, 
0 К 0 
який, за попереднім поданням функції / (и, 2, у), можемо переписати у вигляді 


Ж п п-2 
ч 1 1 
а з п-2 Й ою 2 п-2 п да виз Ме 
о паї п! Р ор ор Кк 
Посилаючись тепер на теорему 2.1 і визначення 








2л р 1 
и(фуг)з Го Пп | заз, 
0 0 


«57 - 25гсов(д -фУ ЧІ" 








переконуємося в справедливості теореми. 

В перебігу доведення теореми існування обкреслено контури множин, у яких варто 
шукати нормальні розв'язки рівняння задачі. Це, перш за все, обмежені множини типу 
(47) в просторі Соболєва Ф (7(0, 27) с 82 (0, 27). Незабаром ми переконаємося, що 
на цих множинах задача поставлена коректно. На інших же множинах, як показують обчис- 
лення якісних ознак |105|, задача визначення області за значеннями потенціалу занадто не- 
стійка. Наведені в (105| чисельні оцінки відхилення двох тіл і відповідних їм різниць зовніш- 
ніх потенціалів цікаві самі по собі в тому відношенні, що дозволяють розмірковувати не 
тільки про характер стійкості задач, але і про неперервність оберненого відображення, а та- 
кож 1 про те, що критерії стійкості даної задачі виконуються винятково на відповідних ком- 
пактних класах єдиності її розв'язків. При цьому останні визначаються повністю теоремою 
єдиності Новікова |421. 

Для доведення теореми стійкості розв'язку задачі (2.3) на Я ""(А") знадобиться одне 


допоміжне твердження щодо залежності між полями та «подібними» областями, що їх обу- 
мовлюють. Тут під «подібними» розуміються такі області 9, 1 7,, границі 09,, і- 1,2, 
яких описуються відповідно функціями 

с (0) с (а) 

с,(а) яс(а)чс, Одхахдя, 
де с -- деяка стала. Введемо в розгляд параметри 

ст і са) 6,- піп со), ше ба» узри", 

що характеризують границю 089, області 9, і її місцезнаходження щодо апроксимуючого 
область круга радіуса р 1 круга радіуса г 2 й, на границі якого задані вихідні дані задачі 
(2.3) (з цієї пори 1 надалі, якщо не обумовлене протилежне, будемо вважати поле заданим на 


колі радіуса ", що охоплює збурююче тіло ? ініде з ним не перетинається). В термінах ви- 
значальних параметрів справедлива 
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Лема 2.1. Якщо поля и" (ф, г), породжені двома подібними однорідними класу Ляпу- 


а) 


нова Зар'"'" областями 9., і- 1, 2, незначно ухиляються одне від одного, тобто задовіль- 


і? 
няють умові 

1; 2 
| и (ф, гу-иЗ(ф, г) І. «є, 


за досить малого є» 0, то і самі області мало відрізняються між собою в сенсі виконання 
нерівності 


| ст(а)-с»(а) Ї. зає, 
де аза(и, у) стала, пов'язана з параметрами границі однієї з подібних областей. 
Для доведення леми досить оцінити знизу різницю 


и (р, гу - м? (р, гу з ФС (а) - ФІЗ, (а))зда: 


на Р З (47), Посилаючись на теорему 2.1 і використовуючи очевидне співвідношення для 
різниці Ф(2, ) - Ф(2,) - Ф"(/)с будемо мати відповідно з прийнятими раніше позначеннями 


параметрів і формулами (2.1а) 


и (р, г - и? (ф, г), 2 Злерії - ші) 





ша У. 
1-и 


Якщо тепер покласти 
-1 


ру) з р змін ) 


то твердження леми буде доведено. 
На підставі встановленої леми доводиться наступна 
Теорема 2.3 (стійкості). Якщо однозначні функції С(а), Очах2, що описують 


границі 09, зірчастих однорідних обмежених областей 5 ,, належать до компактної 
1 жо ; ; 5 з 

множини Я "(А") і задовільняють (крім обмеженості в середньому похідних) додатковим 

вимогам 


«то, 
Го 


46 (ог) 
аа 


а породжувані областями потенціали и"? (ф, г), і -1, 2, незначно ухиляються один від одно- 
го в сенсі виконання нерівності 








(р, ручці (Ф, г) , «є 


за досить малого є» 0), то самі області мало розрізняються між собою, так що згадані 
функції підкоряються умові 


| сі(а)-65(п) Ї- «ає, 
де аза(и, у) -- стала, залежна від параметрів класу ЯР " (А). 


доведення, Складемо різницю 


ою ко дур о0- с (аа зи? (ф)- и? (ф)- ле (2, Ма, 


і визначимо на підставі теореми 2.1 ланцюжок нерівностей 


|тис) -ФІСБУма! - П/С- ЛОСУа є р Па Сад -уа)на є 
0 0 5 К(2-Ф) 
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зл 
«(ру - морі, ПЛР - Гена. ол) 
Звідси видно, що доведення теореми зводиться до ет знизу функціонала 
)-ФіСМа 
і оцінки зверху функціонала 
2 
ПІС)- лена. (2.8) 
0 








Введемо позначення 2), з пах 2, (0), 0»; - піп 2 (8), Су, - тах 2, (8), 2.» - пап 25 (0) 


; бо з нари РИ 


ж -. . 
ких, перетворень на Р " (47), аналогічних до тих, які використовувалися для доведення тео- 


) і на підставі теореми 2.1 за допомогою нескладних, але громізд- 








реми 2.2, знайдемо, що 


Пеєо- Фа з Ї» (ба вбита - раз Їс, (С Сива фаз 

















аа со очарбувк 2 "шва - фо . 





"оп! Ср" 
Згадуючи (2.1) і позначаючи ш ей,р ,уУ-рг,одержимо 
2л у 
Песс)-ФИ на з || 2лр- ик Ю жо) с (а) - сад, (29) 
0 











Визначимо тепер оцінку зверху для функціонала (2.7), з огляду на приналежність фун- 
кцій С(а) до компакту ЯР "( А"). Для цього введемо в розгляд 6 - мережу 
па (а) - б (а), 
ов(0) - СФ-(а)-б 


у вигляді границь «подібних» областей. Варіацію функціонала (2.3) можемо подати в цьому 
випадку в такий спосіб 











ї УРЕ ПИ ПОР БАРИ о 
Пуер - Гома)з р - Га о! г о кину су (реф) 


леви Р 








нему нд) КИ фун нд) Су фун) и 


деФф-ФЯІаФ,а ДФІ «2"'є. Звідси, з огляду на умову теореми, яка стосується оцінки похі- 








лоз тера «тен, а очо| 52 др' пКОрац, 











А так як відповідно до леми 2.1 

сь кИ ні) би кру), З а, у)є, 
і за визначенням дб - мережі 

СИ) СИ Фу. «а(и, у)є, 
то, посилаючись на (2.1а), можемо записати 
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«2лр(1- и) (т-За(и, у))є. (2.10) 


Йиео- тера щі- б 


Для завершення доведення теореми визначимо з ланцюжка нерівностей (2.7) і отриманих 
оцінок (2.9) 1 (2.10) відхилення однієї функції 2, (с), що описує границю збурюючої області, 














від іншої такої ж функції 0, (сг) в залежності від міри є » 0 відхилення потенціалів (створю- 
ваних близькими областями). Будемо мати 
са)-єма, є Ве, 
де поріг Б - Б(2, т, и, у), як неважко обчислити, визначається співвідношенням 
-1 








Кроти ро з тив 1-3 б | 5 
І-и І-м І-м 
-1 -1 
рана), м-н і кін) зі д2лр(- р) цін б) 
ПР П-у) І-у І-у 1-у 

















Стала Б - Б(, т, ш, у) характеризує коректність постановки задачі визначення одно- 
- ої - - ж - го ДК 
рідної області в класі ЯР "(А ) зірчастих областей, обмежених неперервно диференційова- 


ними контурами з обмеженими похідними. 
2.4. Стійкий спосіб розв'язування задачі (2.4) -- заміна задачею мінімізації функ- 
ціонала нев'язки (105). Задачу розв'язування інтегрального рівняння Фредгольма першого 


роду 
А ау, С є (А), ує 820, 2л), (2.4) 


з самоспряженим компактним оператором А замінимо задачею мінімізації квадратичного 
функціонала нев'язки 


кеу-|ла змо (44 - у, А4-т) 
або функціонала типу нев'язки 
Ес(2) - (49 - 2У, с) 
на компактній множині ЯР "(4") з 87 (0, 27). В зв'язку з цим позначимо через и- ВУ 
проекцію точки у є 83 (0, 2л) на множину 
ЯР (А) з ш(;: и є 82 (0, 27), 1є|0,2лР5 из 40, С е я Ф(А7)), 


через 


гсчаге, їщо РО нан ага ідо, РССУ)) гебо,2а) 
точку нижньої грані Е,. - Е(0,) й або ЕС.) функціонала К(2) (чи К.(2)) на Ф "(4)),а 
через 


ФО М С є (А) ге |0,2яр Ф(Д) 2 Е, з іп. Ф(С) 


"Й ) 
- - - - - ж 

множину всіх точок мінімуму функціонала Ф(2) з К(2) (чи К-(2)) на області Ф (47) 
с ЗР) (0, 27) визначення нормальних розв'язків лінійного інтегрального рівняння (2.4). 

Таким чином, задача розв'язування інтегрального рівняння (2.4) 1 задача мінімізації 

, ув 

кожного з функціоналів К(2) і К.(2) еквівалентні одна одній тільки на множині Ф "(4)). 
Множини всіх точок мінімуму функціоналів на 1" (4 ) непорожні і складаються кожна з 
єдиної точки, що слугує нормальним розв'язком заданого рівняння. В той же час відношення 


- - - . - - ж 
еквівалентності між задачами зникає при переході від підпростору Р " (47) до всього прос- 


тору 87(0,2л). Це пов'язано з тим, що властивості розглянутих функціоналів над 
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2 (0, 2л) у значній мірі визначаються властивостями оператора задачі, що діє в тому ж 


просторі. Цей оператор прямої відповідності є обмеженим, неперервним і компактним (пункт 
2.2). А відображення множини вихідних даних на множину розв'язків представляються вза- 


галі розривними в 57 (0, 27). Безумовно, що ці особливості оператора задачі повинні відби- 
ватися в характеристичних властивостях розглянутих функціоналів. Частково їх розкриває 
Теорема 2.4 |78|. Функціонал Ф(С) (К чи ЕК.) над гільбертовим простором 


ЗУ2(0, 2д) є опуклим. Множина ЯР? всіх точок мінімуму Ф(2) незамкнута й у будь-якій 
кулі 5(х.,є) радіуса є» 0 з центром у точці х, є В? знайдеться хоча б одна точка 
сах, нує 82(0,2л), ує 2 (4), відмінна від х., у якій функціонал Ф(С) досягає своєї 
нижньої грані Ф,. Гіперповерхня, що описується функціоналом, являє собою циліндр з твір- 
ними типу параболи, а множина 0, усіх точок мінімуму функціонала Ф(2) над 87 (0, 27) 
- це гіперплощина ЯР? Ф З (А). 

Висновок 1. Задача мінімізації функціоналів Е(2) і К.(С) над С, с 82 (0, 27) неко- 


ректна. 
Висновок 2. Критерій близькості «теоретичного» и - АС, і «спостереженого» у 


полів у метриці простору 87 (0, 2) не є ознакою оптимального розв'язку (квазірозв'язку) 
х, є ЯР? задачі. 

Дійсно, будь-яка мінімізуюча послідовність 42, (с/)), пт Ф(г, ) - Ф, при п-»оо, зав- 
жди збігається до множини М, всіх точок 2, мінімуму Ф(/) на 57 (0, 27), у той час як до 
множини "Р (? єдиності розв'язку рівняння 4/2 - у може і не збігатися при Фе ЕК чи Фо Кк, 
. У свою чергу, оскільки множина 0, усіх точок мінімуму К(2) на 87 (0, 27) -- це гіперп- 
лощина потужності множини 2 (А), то мінімальність нев'язки Ре не може гарантувати 
вдалий вибір розв'язку задачі. Обраний за нев'язкою -/Е, розв'язок С, ех. у, х, є Р), 
ує З (А), може відрізнятися за нормою (5/7 (0, 277) від точного нормального х, на довільно 


велику величину. Тому критерій нев'язки для вибору на 87 (0, 27) ефективного розв'язку 
задачі варто вважати недоцільним. 
У зв'язку з пошуком мінімумів функціоналів К(0) 1 Е-(20) поглибимо вивчення їх 


аналітичних властивостей розглядом диференціальних характеристик цих функціоналів. 
Теорема 2.5 |78|. Функціонали К(С) і К.(С) двічі неперервно диференційовані над 
З (0, 27), їхні градієнти Е"(2) і Е(2) належать до класу С" ( У) з константами ЛІі- 


пшиця І - 2 А А і 1.22 А| відповідно. Кожний з функціоналів над будь-якими елементами 


2, п з 8(0, 2л) задовільняє одній із наступних нерівностей 
ФС) 2 Ф(п) -(Ф"(п), 2 - т), (Ф'2)-Ф"п), с -п)20; (Ф"(2)г), п) 2 0; 
деФеКчи Фен, 
Нерівності, які визначають опуклість функціоналів К(2) і Е.(С), не можуть бути по- 
силені до нерівностей 
(Ф2)-Фп), 2 -пу2 щі - п, (Ф"(2Уль пог шт 


при ши» 0 тому, що оператори А А і А мають потужні нуль-множини 2 (4 4) - З (4). І 





тільки в зв'язку з цим функціонали К(2) і Е.(0) не є строго опуклими; константа // в при- 
ведених вище нерівностях неминуче дорівнює нулю, а множини Ї, мінімуму всіх точок 
Е(2) чи К.(2) над 8/7 (0, 27) необмежені. Виявляється цікавою інтерпретація отриманого 
результату. В банахових просторах умови 
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Ф(С,)-0, (Ф"(2,)е,е)» 0, ех, 


на відміну від аналогічних вимог для звичайних функцій не є достатніми для (глобального) 
мінімуму функціонала Ф - К або Фо Б. Дійсно, якщо глобальний мінімум Ф, функціона- 
ла Ф(/) досягається в точці х, є ЯР; Ф с ,, то відповідно до теореми 2.4 він досягається 1 в 
будь-якій точці 2, з х. 5 у простору 830, 27), де у є З (4). Тому хоча будь-яка мінімізу- 
юча послідовність (2,) функціонала Ф(2) завжди збігається до множини Ї,, до множини 


ЯР? коректності може і не збігатися, тому що в будь-якому околі 5(х., є) глобального міні- 
муму присутні точки з нуль-множини 8 (А) оператора А. 

Множина Лебега 

ЗЕ, є) з Са): |Ає-»Ї «єї; бує 92(0, 2л)) 

за досить малого є » 0 занадто широка. Вона містить елементи, близькі до (точного) норма- 
льного розв'язку х, і як завгодно далекі від нього. Вибір елемента з цієї множини як оптима- 
льного розв'язку задачі, за точністю порівнянного з точністю вихідних даних, можна здійс- 
нити за допомогою коректування мінімізуючих послідовностей функціонала Ф(/) на мно- 
жині 2 - У, ПРО), що містить, очевидно, шуканий нормальний розв'язок рівняння (2.4). 

Для коректування мінімізуючих послідовностей Ф(г ) використаємо стабілізатори на- 
ступного типу 





оче, рад | са а? є ча. 


оо-с, роз | гад са) за леї а | 


які визначені на множині 


іт (Ж, с) з ах с") є 890,27) с" (ачдл)-с"(а), Чо) єс, Ю-0,1) 
за допомогою диференціального оператора ДВ «Ечн ас де Е -- тотожний в 8 (0, 27) 
а 


оператор, Ф «ОО чи Ф-0,, 420 -- число, що означає вагу похідної, с - соп5і . Ясно, що 
множина Ф0) (9, с), як перетин множин Ї/, і ЯР/?, містить точний розв'язок рівняння 
(2.4). 

Властивості наведених стабілізаторів описує 

Теорема 2.6. Функціонали О(2) і 0, (2) на множині Ф)(Ф, с) додатні, сильно опу- 
клі, диференційовані функції класу С (ФУ, мають градієнти класу С" (ФО) і задовіль- 
няють нерівності 








(Ф"(с)п, п)2 





де ФебОчи Ф «б.. 


Доведення. Функціонал ьо додатній за визначенням. А так як, очевидно, 


2 і го ар «і з РИ 


то і функціонал 


б "У 


осде і "оЖека) ча чо а-| Сама - (5, 0)»0 


додатній. Сильна опуклість першого з б аа випливає з того, що для довільних чи- 
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сел а і В з відрізка |0, |, підлеглих умові а -/0 - 1, і довільних елементів 2 (с), П(0/) мно- 
жини Ф0)(Ф, с), функціонал можна записати у вигляді 


дає ві) -« Сас --В п) - (а - а) О(2) - (В- фО( з а О(с)ч В От) - 


2л 
сс п 718 2 Сп 
- -т'за'Сьочбу ат дань я2ат-ЧА, 
ав пут ач (ту 2) 
який на підставі рівностей (2.11) 1 (2.12) 
те й 
Ї 25 Петеуар --(рітФії, 20 (22) 
0 
перетвориться до такого 
(ас - Вп) «аО(2)- ВО(1)-аф|2 - пі. (2.13) 
Міркуючи аналогічно, з тими ж числами а й | і елементами з 0) (М, с), одержимо 
О, (аб я Вп) «а0, (2) ВО, (п) - аа пі (214) 


Для доведення властивості диференційованості розглянемо відповідні прирости фун- 
кціоналів. Якщо функція /(а/) та її варіація л(а/) суть елементи з 9) (Р, с), то, з огляду на 


рівності типу (2.11) і (2.12) і визначення самоспряженого оператора Р « Е- а будемо 
77 
мати 
2л 1 2 
ЛО) - ОС пу - ОД - 2 | р атадіач з |2р?надтадаа, 
0 0 


де 
0 9 а 
РФ-р'раЇ|Е-ач-|Еч -Е-д? | 
2 :2) 1 дай 


З отриманого подання видно, що 





, д'є(а) 


О(2)е2р9а еа2р'ре «2 С(а)-4 зай 


р аура (2.15) 





Аналогічно отримуємо 
2л 1 2л 
ло 22|(атадач з. |2ра)тадаа, 
0 0 


що дає 
Ос(С) - 20(а), Осо) 2 2р. (2.16) 
Для доведення приналежності градієнтів О"(2) і О.(2) до класу С" "(Ф0)) звузимо 
клас розглянутих функцій, жадаючи, щоб елементи 2 (а/) мали узагальнені похідні до друго- 
го порядку, а прирости м(а/) - - узагальнені похідні до четвертого порядку включно. Позна- 
0) » (9 
чимо множини цих функцій відповідно через Ф . (0,2л) 1 ? » (0,27). Безпосередньо з 
формул (2.15) 1 (2.16) видно, що градієнти О"(2) і О/.(2) -- це елементи 87 (0, 27) і тому 
г 3) 
для першого з них можемо записати при 7 є Ф » (0,27) рівності 
2; 


| |ОС нт) - ОС) | - Ап, п), 


2 


СД 


|ОС чт) -О(2)|/-4| 7-4 2 


тому що 
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2 2 2 2л 3 3 дл 2л 4 
дай да | да" да а да бо до о 2 да" 
2? б 2? 2 
р -Е-а б -|Е-ад? Е-а? - рОро, 
і до? Є ба! й 25Ї й 2) 


Звідси, застосовуючи нерівність Коші-Буняківського, знаходимо нерівність 
| нр-ОО|«2|09пД лі 
Посилаючись на цю нерівність, переконуємося, що коли число ря | обмежене, то 
| сн) -О(2)|-»0 | при  |7|-»0. 
А це у свою чергу означає, що градієнт О"(2) є елементом класу С" "(Ф0)), Далі градієнт 
О(2) на множині Є (О,., с) також належить до класу С" "(Ф0)), оскільки 
| Ос(с яп) - О,(2)| 2 2|п| -» 0. 


В той же час для операторів 02"(2) 1 ОС(С) з їхніми зображеннями у вигляді (2.15) 1 (2.16) 
відповідно матимемо 
|«ФОМС зр ОСЬ тр | - 2| (Р - р у, ур | 2 0, 
ФО е тр - ОГС, п) | -2| (Ір - РУ, п)|-0 
і тому О(/), 0, (2) є СФОФ), 

Доведення заключної частини теореми проводиться з використанням прийомів дове- 
дення теореми 2.5 на підставі визначень (2.13) 1 (2.14) сильної опуклості функціоналів. 

Висновок. Функціонали О(2) і О,(2) на множинах ФЗ)(Ф, с), Ф «ОО, не ма- 
ють локальних мінімумів, що відрізняються від глобального, який досягається на елементі 
с(а)єе0, а є|б, 27). 

Вивчимо більш докладно області визначення 90) (ЧЬ, с) функціоналів О(2) і 0,2). 
Зосередимо спершу увагу на області Є")(О,, с) і переконаємось, що вона компактна в 
СІО0, 22| за будь-якого значення постійної с2 0. В зв'язку з цим відзначимо, що на відрізку 
Г0, 27) завжди існує хоча б одна точка Ф, така, що для будь-якої функції С(Ф) з нерівності 


О (2) - с стане справедливим ланцюжок нерівностей 


2л| (ро) є Пеад| «а єс, 


1 


с )|? 
який у свою чергу приводить до обмеження і (Фо) | « 2, . З урахуванням цього обме- 
я 


ження для будь-якої функції 2 (ф) з множини Ф (0, с) можемо записати (на підставі не- 
рівності Коші-Буняківського) 


5 9 р. 
|«(фу|- «Г2) МахЇіеноа| з М з сопіі. 


б У . . 1 
Звідси випливає рівномірна обмеженість множини Ф (0, с). А тому що (з урахуванням 





(ро) | адаа 





тієї ж нерівності) 
1 


б. и 
|Кау -ЄСВ)| - «(а-в)? етіче| «(да- рі)? 





Їс (рудф 
В 
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для всіх функцій ((а) є ФО , су, то ФФ(О,, с), очевидно, одностайно неперервна (у 
метриці СТО, 277). На підставі відомої теореми Арцела |29| множина 92 (О,., с) рівномірно 
обмежених і одностайно неперервних функцій є компактною. 

Для того щоб переконатися в компактності множині (2)(О, с) (у тій же метри- 
ці) достатньо функціонал О(2) подати на підставі рівності (2.11) у вигляді 


я 2 
ор-| сауна 289) ча 
5 а 
і відзначити, що з нерівності О0(/) 4 с випливають нерівності 
2 2 2 
Гоа «с, 2Д 229 да «с 
б о ЗО 
Звідси, модифікуючи відповідним чином попередні викладки, неминуче встановимо 
компактність множини 90) (0, с). 

Звернемо тепер увагу на спектральні характеристики операторів, за допомогою 
яких сконструйовані функціонали О(2) і О0,(2). Виділимо в множині 930, с) підм- 
ножину 

ЗО, дз рух се ЄО, с), с"є 820, 27); с"(ф-дл) - с"(ф)Ь 
і подамо на цій підмножині функціонал О(2) у вигляді 
(2) « (2, Б'рДе) є (2, РО). 
Якщо 2 (ф) і п(ф) суть припустимі функції (з 22 (О, с)), то 


(п, Фо) -(2, рт) зд -підує"(ФіЇ, «Єр 20 


і тому диференціальний оператор 1)? симетричний (самоспряженийї) в ФО (О, с). Окрім 


2 
0 





цього, так як 2 (2)- 0 є єдина за висновком теореми 2.1 припустима функція така, що 
С(а)еб, то при 2 (а) є 0, як стверджується в тій же теоремі, буде (2, ДС)» 0. Отже, 


оператор Др"? додатний. Визначимо його власні числа і власні функції. 
Теорема 2.7. Послідовності 





2 2 Фо 21.2 
1, 14-47, 1-47, ...19ап, 1--4"п",..., 
1 1 - . 
-уч-, зуос05Ф, --58іпФ,..., о сО5пФ, --8іппФ,..., 
2 «Я Я АТ Я 


утворюють повну систему власних чисел ії відповідних їм власних функцій неперервного си- 
: 2 й 2. 
метричного диференціального оператора Др?) « рр в ОО, с). 


Доведення. Власні функції і власні значення ПД)? визначаються наступною граничною 
задачею 


б и( ) ' ' 
о ПФ) - Ам) що) - мл), иЧ0)-м"Ол), 
в якій граничні умови задаються природними обмеженнями на контури збурюючих тіл. Пе- 
ретворене рівняння 


д'и(ф) 5 2. А-1 
оо жф"'и(ф)з0), - ; 
дог В" цф) В Я; 
має такий загальний розв'язок 
иФ) -азіпфф-рсо8воф, а,Р- сопбі. 


З першої граничної умови знаходимо, що а- 0 і со5/)2л - 1. Отже, послідовність 





52, власних чисел визначається у вигляді 2, -144"п', а послідовність (х'")7 власних 


56 
функцій, що їм відповідає, -- у вигляді ірсо5 пФ, . З другої граничної умови одержуємо 
Ь«0 і со8/92л -1, тобто ту ж послідовність власних чисел А, -1ч4'п', п0,1,...,осо, 


яким тепер відповідає послідовність 4 (абіп // Фі, власних функцій.  Позначаючи 


х"(ф) Бсозпф, х"""(ф) завіппф, под,1,..., оо, і вводячи нормування власних функ- 
цій, щоб виконувалася рівність 
-(п) (т) і, пат, 
(оз еде 
0, пет, 


завершимо доведення. 
, КР С КА й 
Згадуючи пункт 2.2, переконуємось, що власні функції (їх ), диференціального опе- 


ратора Д"? збігаються з власними функціями інтегрального оператора А задачі відновлен- 
ня фінітної області. А цю властивість уже можна ефективно використати при конструюванні 
алгоритмів чисельного розв'язку задачі відновлення збурюючої області. 


Розглянемо тепер необхідні характеристики оператора Д) над Ф")(0О,., с). Ясно, що 
диференціальний оператор Д) лінійний в Ф(О,, су. Якщо п(ф) і С(ф) суть припустимі 
функції, то 


(п р2)- (Є, рп) є ЧИТФІС(ФіЇ, - 2 Пп'Фіс (рату з -24 | п'(Фіс (рат є 0 


і тому оператор Р не є ні додатним, ні самоспряженим в Ф"(О,., с) і його спектральні ха- 
рактеристики відрізнятимуться від аналогічних характеристик оператора Д"). 

Теорема 2.8. Послідовності 1--йп); й а/ я е "ХХ утворюють у комплексному гі- 
льбертовому просторі |0, 27) повну систему власних чисел і відповідних їм власних фун- 
кцій неперервного в Ф(О,., с) лінійного диференціального оператора , а послідовності 
(з йті й а/ ру е"?)7. утворюють у спряженому (співпадаючому з вихідним) гільберто- 


. . КУ ж 
вому просторі повну систему власних чисел і власних функцій спряженого оператора Р) . 
Для того щоб переконатися в справедливості сформульованого твердження досить ви- 


. . . . "не . -- ж - 
значити власні числа 1 відповідні їм власні функції оператора 0 (чи Р ) як числа 4, такі 
періодичні розв'язки наступної граничної задачі 


д 
цруча - - и(д), цду з ц(дя), и'(0) з и"(27). 


Повертаючись до пошуку нормального розв'язку 2, інтегрального рівняння 44 - у з 
компактним оператором А у вигляді границі мінімізуючої послідовності і//")) функціоналу 
Е(2) (чи Е.(2)) на ЗР)(0, 27), звернемо особливу увагу на наступні обставини. Множини 
Лебега ЗХФ, а) с 87 (0, 27), обумовлені опуклими функціоналами Ф - ЕК при досить ма- 
лому аа є! чи Фо БЕ. при досить великому а » 0), відповідно до теореми 2.4 є банаховими 


просторами, що містять класи припустимих нормальних розв'язків рівняння (2.4). Однак мі- 
німізуючі послідовності функціоналів ГК (0) і К.(С) на цих множинах можуть не збігатися 


до точного розв'язку 4, (а/) є Р ? (47). Вибір серед них таких, котрі збігаються до множини 


єдиності розв'язків рівняння (2.4) Яр с ЯР" (А"), можна здійснити за ознаками існування і 


стійкості однозначних розв'язків, встановлених теоремами 2.2 1 2.3. Реалізація вибору зв'язу- 
ється з побудовою на ЗХФ, а) непорожніх (за висновком теореми 2.6) компактних множин 


З (Р, с) за допомогою функціоналів О(2) і О, (С), що відбивають характеристичні влас- 


тивості нормальних розв'язків. У зв'язку з цим розв'язок 2, (а) рівняння (2.4) може бути 


57 


знайдений в компактній множині 
З (Ф, Че, 4) з (ба): с є Є, (Ж, с) ПЗФ, а); ФС) «4 ЧУ) « с) 

припустимих функцій. Замикання множини випливає з неперервності стабілізаторів 
Ф «0,0, (теорема 2.6). Окрім того, замість пошуку мінімуму стабілізатора Ч(с) на мно- 
жині 90 (Ф, З; с, а) задачу визначення нормального розв'язку рівняння (2.4) можна звести, 
відправляючись від теореми |73|, до задачі мінімізації стабілізуючого функціонала Ч(с) на 
множині ДФ) (9, сус 83 (0, 27) за умови, що гранична функція 4/7 (а) мінімізуючої послі- 
довності задовільняє рівність АС -у| зє, де є», якщо ухис- Р, ду і6єз-0, якщо 
узв ду. Це типова класична задача варіаційного числення на умовний екстремум і 


розв'язується вона простіше за все методом невизначених множників Лагранжа. Складається 
лагранжіан - - згладжуючий функціонал Тихонова 


Т(с5 а) з Ф(2) на Ж(с) 
і відшукується на (0, 2л) такий елемент 0" (ф), Фе0,271, на якому функціонал 
Т(С) а) досягає своєї точної нижньої грані Т(079) а) з параметром регуляризації а, що 
визначається з умови 


| 4292 -м| (2.17) 





820, дл) ю 
На питання про те, існує чи ні визначена параметром а, функція (79 (ф) з множини 
ОЗ Б, З; с, а) така, що кожний з функціоналів 
Тес(в5 а) З Ес (2) - 20, (2), Тсо(б5 а) 3 Ес (5) на ОС) (2.18) 
приймає на ній мінімальне значення, а нев'язка рівняння 402 - у досягає величини (2.17), дає 
відповідь 
Теорема 2.9 (існування регуляризованого розв'язку). Якщо С, (Ф) нормальний 
розв'язок інтегрального рівняння АС - у з лінійним невід'ємним компактним оператором 
420 у 8? (0, 27), то для всякої функції у(ф) з 83 (0, 27) і будь-якого числа а» 0 існує 


(а) 
єдина неперервна функція С? (ф) з регулярною на |0, 21) похідною и для ФР - Оо чи 
Ф 


з похідною, що інтегрується в квадраті на |0, 22), для Ч - О, на якій згладжуючий функці- 
онал Тихонова 
Т(су а) з Ес (0) а (с) 

досягає своєї нижньої грані на оФар, с) за досить великого с «За ЦМ|/, а сімейство фун- 
кцій (С (фу) із ОФ(Р, с) для й «а ча, збігається рівномірно за метрикою, індукованою 
стабілізаторами Ч'- СО, чи Ф - О, до нормального (точного) розв'язку при а-»0. 

Доведення. Згладжуючий функціонал Тихонова відповідно до тверджень теорем 2.5 1 
2.6 є сильно опуклим диференційованим класу С" (Ф0)У функціоналом на опуклій компак- 
тній множині 9) (Ф, с) і тому для будь-яких елементів у є 81 (0, 27) і числа а.» 0 існують 
його точна нижня грань і мінімізуюча послідовність іс чай (р)) функцій з Ф9)(Ф, с) така, 
що 

т ГО. тс,» реа 


пою п 
Дійсно, умова мінімуму для першої з двох модифікацій (2.13) функціонала Тихонова 
має вигляд 


(Азавуг? зу, ГО -0, (219) 


а для другої - - 
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(Ачарус? щу, 29, -2-4 б, (2.20) 
де |, 


Оператори кожного з цих рівнянь Ейлера строго додатні при є2»0, тому що 
ТС ура)» ої Т(С9, у, а) » 0, і у згоді з альтернативою Фредгольма кожне з рівнянь 
має єдиний розв'язок С" (ф) за фіксованого са» 0 і за будь-якої правої частини м(ф) з 
ЗР (0, 27), в тому числі такої, що задовільняє точному розв'язку задачі. 


Існування мінімізуючої послідовності обгрунтовується в такий спосіб. Нехай (у, (Ф)) 


9 








-- послідовність функцій з б87(0,27) така, що послідовність (є,) чисел є,- |У, -и 
из Роду» Збігається до нуля при п -» оо. Тоді із-за неперервності і сильної опуклості фу- 


нкціонала Т(/, у; а) послідовність (0, (ф)) розв'язків 2,7 (ф) послідовності рівнянь Ей- 
лера Т"(С, у,, а) - 0 буде, очевидно, мінімізуючою. До того ж можна вважати, не обмежую- 
чи узагальнення розгляду і посилаючись на опуклість функціонала, що для кожного п21 
буде справедливим наступний ланцюжок нерівностей 
т « «Т(оп) «Топ « 
А на підставі залежностей (2.19) 1 (2.20) матимемо 
(ап) -1 (а, п) - (295-1 
Тесс" азу (А НИЄ) у,) «0, То "з чу, (АчабУ?) у, ) «0, 
тому для будь-якого п 2 1 і довільно фіксованого с/ » 0 виконуватимуться нерівності 
(а) -ГАА52 (а) -4442 
ОС )«2а |, ОС )« а "М. 


п 


«ори 








З цих оцінок випливає, що кожна з послідовностей (4,7 (ф)) функцій 2, «(Ача)"!у, чи 


й. ЩІ «(Азчар'ру"», належить відповідно до множини 27 (О,, с) чи ФС/О, с) рівномір- 
но обмежених в середньому і одностайно неперервних функцій, тобто 
шо 279) -0799| -0 

б є (0, дл) 


пою 


Переконаємось тепер, що ці послідовності функцій фундаментальні за метрикою, ін- 
дукованою стабілізаторами О,.(2) чи О(2), тобто в просторі Соболєва Ф (7 (0, 27). Зосере- 


димося спочатку на мінімізуючій послідовності функціонала ТТ. (С, у; с). Так як для будь- 


яких цілих індексів 1, / -і--К 2 1 з урахуванням рівності (2.11) можемо записати 


2 
она - С) о 2 У ; 














то 

219) 49, ре) СОУ «Аза», -у ої «а"Че,не,) є Дей, 
тому що за визначенням є, 5 є, Звідси при є, -» 0 видна фундаментальність послідовності 
2/9 (ф)) функцій 2,7 (ф) і їх похідних (ДС? (ф)) у сенсі слабкої збіжності в Ф (0,27) і 
існування слабких границь функції 2 (З(ф) і її похідної в силу повноти Ф (7(0, 27) як гіль- 
бертового простору. Так само і для мінімізуючої функціонал То (С, У, с) послідовності 


а 4 (Ф)) визначаємо 
2 


От? -6)0)9| 





Ра 








ж (20,22) 
5 иа 4 
Ла В В Са У СЕ о об 


та переконуємося, що послідовність (2, (ф)) функцій 2," «(Ачар'Ру!у, при є -» 0 


фундаментальна в  Ф (2(0,2л) і тому сильно збігається до граничної функції 


За 


й «є Ф Ф(0,2л). За побудовою кожна з мінімізуючих функціонали Р. (2, у, а) і 
Коо(С, У, Ф) послідовностей (с й (Ф)) крім того рівномірно збігається за метрикою 
ЗРУ(0, 2) до екстремалі 0" (ф). А тому що Ф (7(0, 2л) с 87 (0, 27) і екстремаль кожно- 
го з функціоналів І. (2, у, а) чи (2, у, а) єдина, то С  (ф) е С (ф)є 99 (ФР. с). 
Для завершення доведення теореми залишається показати, що кожен з операторів 
А"(у,а)з(Ачавбу» чи АЧу а)«(Ачар'ру'» 
зі значеннями в Ф (7(0, 27) є регуляризуючим для рівняння 42 - у, або що регуляризовані 


розв'язки 0/7 (ф) за досить малих а » 0 добре апроксимують нормальний розв'язок даного 








рівняння. 
З попередніх міркувань ясно, що 
2 
(71 а а Є 
о -бн) а т ня ре -бн)) « ук й 
2 
(24 а а Єї 
ОС ну (ФСУ зби ВСІ си) Є (2.21) 


Наступні висновки щодо цих нерівностей істотно спираються на таке твердження 
Лема 2.2. Якщо функціонал ача(є), 0«єх«є,, визначається на множині 


ЄСК, єн: Їу-иі «є МАВ ду), елементів уУ(Ф) з ЗР (0, 2) нев'язкою є - Ї» - 409 

лінійного інтегрального рівняння Ас - у з невід'ємним компактним оператором, то його 

нижня і верхня грані на ЗК, Єв) оцінюються відповідно залежностями 

3? 5 Є 
, 


та ав(в) є Б'є?'ї, 
а 


ан(є) з 


(а) 
- )Ї5Ь? -- додатні числа. 
а 





(а) 
де а! « п СТАНИ 
да 


З умови леми випливає, що є «(у - 40), у- АС)зф(а), де С? (ф) -- екстре- 


маль функціонала Т(/, у; а) на Є) (М, с), Звідси визначаємо диференціальне рівняння 














Фа 

2є4є- 90 аа, (222) 

у якому 
(а) 

да да 

Посилаючись на визначення (2.19) 1 (2.20), цю похідну можемо переписати у вигляді 
(а) (а) 
о у Є (2.3) 
да да да 
звідки переконуємось, що існує таке додатне число а", яке задовільняє умові 
(а) (а) 
а" «(А св а | 
да да 


З урахуванням цього зауваження на підставі залежностей (2.22) 1 (2.23) після нескладних об- 
числень знаходимо нижню оцінку для функціонала єа-с(є), 0«єх2є,, на множині 


ЗА , Є) Якщо ж вираз (2.22) подати в такий спосіб 


сне - за анаву оС р идчануб о ) 47 с у анавубо ) 
Г2И 











а да да да 
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і оцінити його зверху, то прийдемо до визначення с,(є). Отже, функціонал а -с(Є) на 
множині ЗК, Є) Є монотонно зростаючою функцією, неперервною на відрізку |0, є, і 
такою, що с(0) -0. 

Повертаючись тепер до визначень (2.21) і приймаючи в них для оцінки параметра с 
границю а, , установлюємо, що 


С нь С Су) є 2) 50 -»0, 


з 
ре ен «(29 є оз0 


9 Ф(0,2я) 3 
при |» -ці- 0. в і ість (2,9 (фу) б В 
ри |у, - ці - Є, -» 0. Виходить, послідовність 16, (Ф); досить добре апроксимує нормаль 





ний розв'язок с, (Ф) при с є|0, с. |, де с з (є,). При цьому кожне сімейство функцій 
С ка (ф)), що відповідає будь-якій послідовності (є,), що прямує до нуля, збігається при 
а-»0 у випадку Ч «СО, слабо, а у випадку ФР «ОО -- сильно до функції 2, (Ф), що має 
узагальнену похідну, за метрикою простору  Ф (0,27). А так як відображення 
ЛОТО, су з ФА) с З(0,2л) неперервне, нормальний розв'язок 2, (ф) рівняння 
АС су за теоремою Новікова |42) єдиний для всякого у є Я (А) і множина Ф0)(Ф, с) 
компактна в 2 (0, 27) для досить великого с» 0, то відповідно до теореми Тихонова |62| 
обернене відображення В " (А) -3 Ф)(Ф, с) також неперервне, Це в свою чергу означає, що 


для кожного 6 » 0 знайдеться таке є - «(0)» 0, що з нерівності |» - ці є (6), у, м е ЯР (4) 


суди 














буде випливати нерівність «б, якщо у - А из АС,. Теорема повністю доведена. 


у? 
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2.5. Регуляризуючі алгоритми. Задача визначення границі тіла 9 , близького до кола 
Ф., радіуса р, зведена до побудови послідовності (0, (са)) наближень 0 (с) границі 
с(а) є р(8)- р у вигляді розв'язків послідовності лінійних інтегральних рівнянь першого 
роду (2.3), де и(ф, г) -- поле на колі ДФ :х'-у"' зх радіуса г» р,а Г(а-ф,б,;р,г) 
- функція, визначена рівністю (2.2). Точний розв'язок 27 (ф) кожного рівняння послідов- 
ності (2.3) при точно заданій правій частині у" (ф) можна шукати у вигляді розгортання за 
власними функціями оператора задачі з ядром К(а., Ф; у), тобто, у вигляді 


с, (фузад я), а собер- ВІ" "зіпкф, 
Каї 


оскільки власні функції оператора задачі утворюють базиси в підпросторах ЯР (4 ) і Р (А) 


гільбертового простору 3/7 (0, 27). Розкладаючи в цьому базисі у" (ф) є ЯР (А) в абсолютно 


збіжний ряд 


У"(фу зач зало сор"? віпф (2.24) 
Каї 
одержимо 
(п, Т УГ ут 
(п, Ту. а? (п. Ту. ка?" сьту. КБ? 
2 2 
"  длОюї) 7 4 4 | 


Ясно, що цей розв'язок в умовах, коли функція у"" (ф) задається з похибками, практично ви- 


користовувати не можна через сильну нестійкість підсумовування рядів Фур'є з приблизно 
заданими коефіцієнтами. 
Приступаючи до реалізації розв'язку задачі (2.3) у виді пошуку послідовності 12, (Ф)) 


границь С, (Ф) кожної з мінімізуючих послідовностей і2,,(Ф)) згладжуючого функціонала 


Т (С) а) для п-го інтегрального рівняння послідовності (2.3), перш за все відзначимо, що в 
зв'язку з завданням замість точного значення потенціалу и(ф, г) на колі 20, наближення 
их (ф) зих(Ф, г) і відповідного йому наближення правої частини УЧ" (ф) рівняння у вигляді 


відрізку ряду (2.24) 
М 
уг (фу за? - У ар" соб -нЬ!"? віпкф, 
Кеї 


кожне з наближень 2, (ф) мінімізуючої послідовності функціонала Т""(2) а, ) будемо від- 


тп 


шукувати у вигляді 


М 
т Ф) з а У а сор ВУ" "віп. 


Кеї 
Це забезпечить абсолютну точність визначення згладжуючого функціонала на даному еле- 
менті, оскільки система ісо5 КФ, 5іп КФ), , Є власними функціями не тільки інтегрального 
оператора задачі, але і диференціального оператора, за допомогою якого з 7. (0, 27) відо- 
кремлюється підмножина ЯР З (А") коректності для кожного інтегрального рівняння типу 
(2.4). 
З умов мінімуму на множинах 20) (9, с) і відповідно до напрямку іс, ) для кожного 
з функціоналів 
Тез о.) З То у?) а (С, АС, 2) ча, (2, Де, 
То (85 и) То ра) б, 46, ду а ФС, ре), (2.25) 


п п 9; 
уро (аб о) -Тро(б я» у па -|АС, ТАЙ таб,» Ре,) 
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одержуємо 
(п) (п) 
(туп) ау (туп) 0 КА (ту КБ, 
0 - 2 Кк б 2 Кк г 
2л(-Іпг)ча, лу зак лу"за, К 
для функціонала Ту (С, У"; а, )Ї 
(п) (п) (п) 
(т, п) | ад (т, пу | Ка (т, пу | Кк 








є з , Рі - 


9 2л(-пг)ча, лу чо, Кз) 


для функціонала Т.о(0, У") а, ) у і, нарешті, 


(п) (п) (п) 
(туп) 2л(-Пп гад, ОТО лукаї нібоо лукр 


, Я 7 бл 
і їгл(- УР ча (лу У ча, КАЧКУ" (лу)? за, КПз 70) 


для функціонала Тоо(С,, У"; а, ). В усіх трьох випадках регуляризовані коефіцієнти розк- 


лу за, Казча у! 











ладання 2,,(Ф) при а, - 0 тотожно дорівнюють абсолютно точним коефіцієнтам (2.24). У 
той же час складність обчислення регуляризованих значень коефіцієнтів зростає в напрямку 
визначень від функціонала Ту(0,, У"; а, ) до функціонала Тоо (2, У) а,). У тому же на- 
прямку через округлення зменшується точність обчислення коефіцієнтів. 

Оптимальні значення коефіцієнтів Фур'є знаходяться методом нев'язки. Задається 
спадаюча послідовність и, ) значень параметра регуляризації, наприклад, у виді спадаючої 
геометричної прогресії 1 за допомогою методу дихотомії відшукується таке значення параме- 
тра а зд,» ЩО задовільняє (наближеній) рівності 


(п) 
(мб (а) -у 
де є - - середня квадратична міра похибок. 
Якщо оптимальні коефіцієнти знайдені, тоді складається наближення 2, (Ф) у вигляді 








-є, 


відрізку ряду Фур'є й обчислюється наближення функції 
2л 
Фр) - | Фа, фСзр га, 
0 


що визначає наближення потенціалу би(фФ, г) «надлишкового» тіла ДФ -919,, і знахо- 


0? 
диться величина відхилення 


АФ) Фр) - ру" (рі. 
Якщо ЛФ"? «є, то п-е наближення 4, (ф) границі С(ф) приймається як остаточний 


розв'язок задачі; якщо ж ЛФ? 5 є, то здійснюється перехід до обчислення п -- 1-го набли- 
ження 2, (Ф) границі тіла ДФ? за наведеною схемою обчислення п -го наближення 2, (Ф) 


тощо. У зв'язку з пошуком наближення 2, (Ф) ще раз підкреслимо, що спроби визначити йо- 


го на відміну від побудови мінімізуючих послідовностей функціоналів типу (2.25) за допомо- 
гою мінімізації функціоналів нев'язки з додатковою вимогою мінімальності норми розв'язку 
практично приречені на невдачу, тому що ці вимоги не обмежують рамки пошуку розв'язків 


на компактах ЯР" (А). 


Висновки. Виведено нелінійне інтегральне рівняння з цілком неперервним (компакт- 
ним) оператором, що визначає за значеннями потенціалу відхилення контуру збурюючого 
тіла від заданого. Для пошуку розв'язку рівняння в припущенні, що відхилення невеликі, за- 
пропоновано ітераційний процес, на кожному кроці якого наближення відхилення обчислю- 
ється у вигляді розв'язку лінійного інтегрального рівняння з симетричним невід'ємним ком- 
пактним оператором. 

Доведено локальну теорему існування і теорему стійкості розв'язку задачі. Вказано 
магістральний напрямок пошуку нормального розв'язку задачі. 
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Розроблено ефективний метод розв'язування нелінійних інтегральних рівнянь, що 
описують відхилення контурів зіркових збурюючих областей від заданих. Глобальний регу- 
ляризуючий оператор для нелінійного рівняння відшукується у вигляді послідовності лока- 
льних для відповідних лінійних інтегральних рівнянь. 


3. ПРО ПРОЦЕСИ ДООРТОГОНАЛІЗАЦІЇ ДЕЯКИХ СІМЕЙ ВЕКТОРІВ 


Вступ. Основні задачі, які визначили в кінцевому рахунку розвиток класичної лінійної 
алгебри, пов'язані з визначенням власних значень і відповідних їм власних векторів деяких 
матриць п - го порядку, а також з розв'язком систем лінійних рівнянь. Нагадаємо, що влас- 
ними значеннями (або числами), матриці п - го порядку називають |З, 37, 68, 69| корені її 
характеристичного полінома р, (ЛА) п - го степеня, тобто корені рівняння 


ри(А) « фе (4-Л ЕФ) «|А- ЛЕ"? 











ай "А ПРОДУ ПРУИРИ ам 
й НО рів аз 
РУ ПАРИ РР Р аз (3.1) 
б 77 ЛО о антЯ 
ззуЧ-ад'-а4-..-адеб, 
де В"-(б, і) - - одинична матриця, діюча в п - мірному просторі щої і; з СИМВОЛ Кро- 


некера. Власними векторами и, і - 1,2,...,п, матриці А називають нетривіальні вектори 


-- розв'язки системи п однорідних лінійних рівнянь з п невідомими, тобто розв'язки систе- 
ми 


(4-4, Е"Уи «0, (3.2) 
і ж ки . 
де ці? - (п з Мої-зМ,;) -- власний вектор матриці, що належить власному значенню Л ,, 


і - 1,2,...п, а зірочкою позначена (і цього позначення дотримуватимемося надалі) операція 


транспонування (спряженості) вектора - стовпця (або матриці). Тут же відзначимо, що не для 
кожної матриці, заданої над полем дійсних чисел, знайдуться п дійсних власних значень Л, 
і відповідних їм власних векторів и, 

Центральною проблемою при визначенні власних чисел є відшукання коефіцієнтів Л, 
, і - 1,2,...,п, характеристичного полінома. Це зумовлено тим, що безпосереднє визначення 
коефіцієнтів Л, із співвідношення (3.1) у вигляді сум всіх мінорів 7 - того порядку визнач- 
ника матриці, що спираються на головну діагональ, не тільки надзвичайно громіздко, але і 
вимагає величезного числа операцій особливо при великих розмірах матриці. Остання обста- 
вина в свою чергу приводить до відомих труднощів навіть при використанні сучасних ком- 
п'ютерів. Тому для знаходження 4, були створені спеціальні обчислювальні методи, в яких 


пошук коефіцієнтів Л, реалізується без стомливих 1 безперспективних обчислень численних 


мінорів визначника матриці. Першим серед них був, певно, метод У. Дж. Леверьє (1840), 
який залишався ще досить трудомістким. Більш витончені і менш трудомісткі методи з'яви- 
лися значно пізніше, майже сто років потому. Це -- методи А. Н. Крилова (1931), А. М. Да- 
нилевського (1937), К. Хессенберга (1941), Д. К. Фаддєєва (1949), К. Ланцоша (1950-1952) та 
інших. 

Обчислювальні аспекти цих і багатьох інших методів глибоко проаналізовані в осно- 
воположних творах К. Ланцоша |37|, Д. К. Фаддєєва і В. Н. Фаддєєвой |б69Ї, 
Дж. Х. Уїлкінсона |63| 1 В. В. Воєводіна |12|. В цих роботах після детального дослідження 
впливу помилок заокруглювання проміжних обчислень на точність остаточних результатів 
сформульовані загальні принципи складання стійких алгоритмів, які реалізовують на ЕОМ 
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той або інший метод. Однієї з важливих, якщо не найважливішою особливістю таких алгори- 
тмів є врахування того, що арифметичні операції над числами в комп'ютерах (із-за обмеже- 
ності діапазону використовуваних чисел і представлення їх кінцевим числом розрядів) не пі- 
дкоряються звичайним законам комутативності, асоціативності і дистрибутивності. В цьому 
плані різні методи визначення коефіцієнтів 4, характеристичного полінома (3.1) можна кла- 


На підставі міркувань, про які мова піде нижче, можна сподіватися, що найточніші оцінки 
2.; можуть бути отримані при реалізації методів, заснованих на ортогоналізації цілком ви- 
значеної послідовності векторів, тісно пов'язаних з досліджуваною матрицею. Тому серед 
всіх методів непрямого визначення коефіцієнтів 4, перевагу потрібно віддати методу міні- 


мальних ітерацій К. Ланцоша 1 йому спорідненим методам. 

А зараз -- декілька слів на виправдання зробленого вибору. При реалізації на комп'ю- 
тері методу мінімальних ітерацій або будь-якого іншого методу, заснованого на ортогоналі- 
зації деякої системи векторів, доводиться обчислювати рекурентно одна з одною зв'язані за- 
лежності, в яких в якості коефіцієнтів використовуються значення скалярних добутків векто- 
рів. При такому способі підрахунку неминуче накопичення помилок заокруглення, які мо- 
жуть в значній мірі спотворити остаточний результат. На щастя, використання властивостей 
ортогоналізації векторів дозволяє виправити отриманий результат. Відзначимо тут же, що 
подібна ситуація, коли в процесі обчислень вдається знайти "жорсткий репер" (в даному ви- 
падку -- ортогональність векторів), зустрічається не часто. Принаймні, в методах У. Дж. Ле- 
верьє, А. Н. Крилова, А. М. Данилевського, К. Хессенберга, Д. К. Фаддєєва і інших ми їх не 
знаходимо. В зв'язку з цим слід очікувати, що оцінки Л, в методах, заснованих на ідеї орто- 


гоналізації, будуть кращими в порівнянні з такими в інших методах. 

Ідею відновлення втраченої при обчисленнях точності результату на випадок ортого- 
налізації векторів запропонував К. Ланцош одночасно з описом методу мінімальних ітерацій 
1116, 117). Відповідно до його пропозиції покращення ортогональних властивостей деякої 
системи векторів, близьких до ортогональних, повинно проводитися в результаті повторної 
їх ортогоналізації. При цьому поправки, що поліпшують ортогональні властивості векторів, 
знаходяться з певної системи лінійних рівнянь. Уточнюючі в деяких аспектах ідею К. Лан- 
цоша пропозиції були висловлені в роботах |9, 12, 115|. В статтях |33, 95| ідея К. Ланцоша 
розвивається з інших позицій, які дозволяють уникнути розв'язку допоміжної системи ліній- 
них рівнянь. При цьому в роботі |95| для випадку ортогоналізації сім'ї векторів методом 
Грама -- Шмідта процес ортогоналізації поданий в замкнутому вигляді. Розвиваючи викорис- 
товуваний в |95| підхід на інші способи ортогоналізації, супутні різним канонічним зобра- 
женням дійсної матриці п - го порядку, ми отримаємо можливість вирішити задачу: визна- 
чити з високою точністю коефіцієнти характеристичного полінома заданої дійсної матри- 
ці п - го порядку. 

В свою чергу розв'язок поставленої задачі, як відомо | 15, 37, 68, 69), забезпечить на- 
дійне обчислення, по-перше, власних чисел матриці, по-друге, її власних векторів 1, нарешті, 
по-третє, (узагальненого, нормального) розв'язку системи лінійних рівнянь з розглядуваною 
матрицею. Насправді, при побудові анулюючого матрицю полінома методом К. Ланцоша або 
йому спорідненими на проміжних етапах конструкції ми одержуємо послідовність поліномів 
Штурму, що дозволяє розділити корені полінома, а для їх наближеного обчислення викорис- 
товувати або ітераційний метод І. Ньютона, або більш швидко збіжний процес наближень П. 
Л. Чебишева |8, 37). Отримавши власні числа і маючи послідовність поліномів Штурма, по 
відомих співвідношеннях |69)| легко відновлюємо власні вектори матриці. Після цього не 
складає труднощів розв'язати систему лінійних рівнянь з цією матрицею при будь-якій пра- 
вій частині системи | 15, 69. 

Перш ніж перейти до детального розгляду сформульованої проблеми відзначимо одну 
важливу характеристику с( А) матриці 4, що відіграє фундаментальну роль в питаннях стій- 


кості розв'язків систем лінійних рівнянь, яку називають мірою (або числом) обумовленості 
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матриці і визначають співвідношенням 
сСА) з | АЇА" 








2 








де | А, 14" | - деякі норми матриці А і її оберненої А, узгоджені з нормою даного векто- 
рного простору. Якщо в цьому просторі збурення Д, і Д, вектора фФ і матриці А системи 
Ахаб 


визначаються відповідно величинами А М | і Їд 4Ї, а їх відносні похибки --- числами 



































Р 1 р 
ІБ |А 
то, як відомо |9), відносна похибка розв'язку системи оцінюється нерівністю 
ДЛ. 
Я і, С бАчОБ) 33) 


Її 1-ссАбА 
Звідси витікає, що відносна похибка розв'язку збуреної системи лінійних рівнянь залежить 
значною мірою від числа с( А) обумовленості матриці А: з одного боку, чим більше це чис- 
ло відрізняється, скажімо, від одиниці, тим нижче відносна точність розв'язку системи. З 
другого боку (і на це вказує знаменник в правій частині оцінки (3.3)), чим більше міра обу- 
мовленості, тим з більшою відносною точністю повинна бути відома матриця даної системи 
для забезпечення існування її розв'язку. Таким чином, з погляду стійкості вельми істотно, 
щоб число обумовленості не було надто великим. Проте з рівності 
А'АЗАА ЗЕ? 

витікає, що незалежно від вибору норми матриці ці числа ніколи не можуть бути занадто ма- 
лими, принаймні, будь-яке з них (залежно від вибору норми матриці) не може бути менше 
одиниці. В роботі |95| виділені класи перетворень початкової системи, що не міняють її обу- 
мовленості. До їх числа відносяться ортогональні перетворення. Надалі буде показано, що 
міра обумовленості матриці відіграє істотну роль не тільки в питаннях стійкості розв'язків 
систем, але 1 в розв'язку сформульованої проблеми. В світлі останнього зауваження і голов- 
ним чином того, що клас ортогональних перетворень не змінює міри обумовленості матриці, 
слід ще раз підкреслити те виняткове положення, яке займають методи, засновані на ортого- 
налізації деякої системи векторів, серед інших методів непрямого визначення коефіцієнтів 
характеристичного полінома. В обчислювальному аспекті ці методи виявляються найстійкі- 
шими. 


3.1. Метод ортогоналізації послідовних ітерацій 169). Будемо розглядати дійсну ма- 
трицю 4, діючу в п - мірному евклідовому просторі 2", котра має п різних дійсних зна- 
чень 4, і відповідних їм п власних векторів и"), і- 1,2,...,п. У відповідності до теореми 
Келі - Гамільтона матриця А являється коренем свого характеристичного полінома р, (Л), 
тобто вона задовольняє рівності 

р.(4) з А" ча ага"? -.. за, Аз, (34) 
Цю чудову рівність А. Н. Крилов використовував при визначенні коефіцієнтів /, полінома 
р, (2) шляхом попереднього перетворення р,(Л) в йому еквівалентний поліном Ф, (Л). 
Перетворення, за допомогою якого з р, (А) виходить Ф, (4), задається в базисі 


2 -1 
и, Ли, А"и,..., А", (33) 
де ц є ХХ") -- такий вектор, що його ітерація А"и матрицею А належить лінійній оболонці 
З (и, Аи,..., А""и) тільки при Кк п. Іншими словами, вектор и повинен підбиратися так, 


щоб система векторів (3.5) була лінійно незалежною. Тут 1 впритул до п. 3.5 вважатимемо, 
що вибір вектора и задовольняє поставленим умовам. 
Плідна ідея А. Н. Крилова про перетворення р (Л) в Ф (Л) за допомогою базису 
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(3.5) по суті використовується в усіх методах непрямого визначення 4 ,, які таким чином ві- 
дрізняються один від одного тільки тим, як визначаються компоненти вектора Ав цьому 
базисі. В розглядуваному нижче методі компоненти | А"у в базисі оболонки 5 
(и, Ли, ..., А""и) визначаються в процесі ортогоналізації сім'ї ГАЮ Ти) Мр з ЙО ЛИЮНЯ 

На відміну від опису |69| викладемо обчислювальну схему методу для ортонормованої 
послідовності векторів и"), і - 1,2,...,п. В якості початкового вектора ортонормованої сис- 


теми приймемо 





и 
чо ша обеійо з Їм а и, м) з 


Нехай перші К ортонормованих векторів ці, и.,...,и, вже побудовані і 6) 
(ць Ай оз 4 щу) з (цу цу зи,). Тоді вектор и,,, будуватимемо так. Обчислимо іте- 
рацію Аи, вектора и, матрицею А ірозкладемо Аи, на дві складові: 
ча Пен) Но Є сь); 
а для відшукання цих проекцій складемо лінійну комбінацію 























К 
аа по | 
Ша З Ша Мра З Чи 2 Вій : (3.6) 
ізі 
Коефіцієнти /3,, - - компоненти проекції и;,, в базисі Є (й, и»,...,и,) - - Знайдемо з умов 
ортогональності: 
Ви (Ашьц), і-1,2,.2зв, 
Вин Їм 2 ке іньу й, (3.7) 
В силу вибору початкового вектора повинно бути 
/ 
Ін 20 
для всіх К - 1,2,..., п-- 1, що дозволяє визначити вектор и,,, у вигляді 
' 
М 
ща, ї» Кк аї1,2,...,п-і1. (3.8) 
ші 
К-1 - - - 
При цьому за побудовою (и, Аий,... 4" и) з Є (ці и»,...,и,) і, крім того, внаслідок 








рівності и - 0, що витікає з (3.4), одержуємо 


0 РР ЛОВ (3.9) 
ізі 
Для кращого огляду випишемо формули (3.6) 1 (3.9) більш детально 
шцо-й 
Визьзяфиця й, 
Вода Вищої НА, 


за аьий пів і і ос блін (3.10) 
Ваня - РЕ, й -финль т АВ ранній їз Ай 
Ос еф лій -Вь М. Р 2 ся ди, 
Звідси слідує 
1 
рабу ФАСАДИ, к-012.гап, (3.10 А) 


Кк 
де поліноми Ф,(7) обчислюються, як це витікає з співвідношень (3.10), по наступних реку- 
рентним формулах: 
Ф1)-1 
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р) «сову то 
1 Йо 
ФьЬІ) Чо) 3 Оз та 
ПОРТИ РО НЕИРЕОНРИОНОРРИ РІВ ПРИРРРНИРРНЯ (3.11) 
1 - віє 
о(5-а-Ви)з (ФьО-2 от) 


2К0-(-8,99,К0-У вить) 


ке дна П--Ї; 

Таким чином, вирахувавши при вдалому виборі початкового вектора и, всі коефіцієн- 
ти В, і -1,2,...,К--Ї, Кк е1,2,...,п, по формулам (3.7), ми можемо послідовно знайти всі 
поліноми Ф (1), Ф (1), ..,Ф, (1). У випадку, якщо процес обчислення /3),, закінчиться завча- 
сно, тобто коли при К « п виявиться /3,, ,- 0, ми отримаємо мінімальний анулюючий век- 
тор и, поліном Ф, (г). 

Легко показати, що Ф (7) з точністю до постійного множника співпадає з характерис- 
тичним поліномом р, (7) матриці А. Замітимо на закінчення, що співвідношення (3.10) мож- 
на записати у вигляді 


АС -йв, (3.12) 
де (яні, Мо. ---, М, Ї, 
Ви Вь 0 0 рнабоннмні 0 о 1 
Ва Во В 0 наново 0 0 
В Ва Ва бі нні 0 0 
ВІ акрбеннскй занільвовиоко зеборфибієюнюх сЙнербвбеново (вбфтркбвйнон» «Нрюрраннийна обенажов Фк» 
Вл В 2 Й аз Віа раазаланадеаналецькося Па Ро ДАН 
Ва Ва Віз Від е'жзвсанв вовивовявато ЙО і Вип 








Звідси витікає, що матриці А і В подібні, бо 
АЗОВУ"зОВО, 
а це в свою чергу означає, що анулюючі поліноми матриць А 1 В однакові. 


3.2. Метод мінімальних ітерацій К. Ланцоша |З, 69). Метод мінімальних ітерацій 
для вирішення проблеми власних значень є не що інше як описаний в попередньому пункті 
метод ортогоналізації послідовних ітерацій деякого початкового вектора симетричною мат- 
рицею. 

Отже, якщо А: 52") -» СХ") -- симетрична дійсна матриця, діюча в п - мірному евклі- 


. ж го 
довому просторі 2" ,а це (ці Цо»г-оц,) Є 8) -- вдало вибраний вектор, то процес ор- 


тонормалізації послідовних ітерацій и, Ац,..-, Аи, Кк -1,2,...,п, вектора 
и 
ил иЇ чАКи,и 
1 Їм дою 7 | | кт) ( , ) 


матрицею А здійснюється по тричленним рекурентним формулах 
Ванна Вірна Вк КМ, (3.13) 
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Ка2,3,...,п- 1, 
апри Кк «11 Кк - п відповідно по формулах 
В ильззВиш Ай 


0-3 анна ПВ пл МИ, (3.14) 
де 
Віт (4и, зи,) я 0 для і- «ЇХ 
В. я (Ми, зи) «0 для - 21. (3.15) 
Насправді, з (3.10) маємо 
Вічна З ВушоВуль оо Ву мк А, 
Вина С "ВищоВоль го В Ми, (3.16) 
Оскільки одиничні вектори попарно ортогональні, а матриця -- симетрична, то з виразів 


(3.16) отримаємо 
В цу» и,) - у (Ап, и,) ач Ву (цу Аи,) 5 

з Вау (ци В ни ЯР ул» жом, ЯВ зай) (3.17) 

А оскільки и,,Ї. (цу ггьи,), / 2 1,2,..., К, то з рівності (3.17) негайно одержуємо 
Вуч0, /-152, аз -2, 
"Вон Вект 0, і2-Кк-і, (3.19) 
Ват(4ии), /-к. 

Далі з рівності 














К 
Вика - 3 Вил - Ли, во? 0, к-1,2,...,п- 1, 
ізі 
з врахуванням виразів (3.18) знаходимо 
Ви Вина (Ми м, |) 2 0. (3.19) 


Тим самим вирази (3.13) - (3.15) з врахуванням симетричності матриці повністю доведені, 
так як для будь-яких ї, ) - 1,2,..., п мають місце рівності 
(Ли, и, ) с (и, Ми). 
Отже, ми прийшли до співвідношень 
Віщьзояфиця4ц, 
Вод зоРФощофоть Ай, 


орел вбфевОальсиь Ко а сова (3.20) 
Ва фу Віа і РН 
0-й рай "В лай из 
які по аналогії з рівністю (3.12) можна записати більш коротко: 
АС Ов, (3.21) 
де (зи, Мо, ---, М, І, 
У» об Ю.О зада 0 0 0 | 
Ва Во Ваз ооое 0 0 0 
Фо фо оба кіна 0 0 0 
Вер зефавно! забою сабіьіюо Узабді» 0 з ВИНЙ 0 ВІВ 
000... дно орні а 
0 0 9 ДО 0 а со 
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По побудові очевидно 09 (и, Ац,,..., А" цу) в (ци, и,). 

Таким чином, ми переконалися в тому, що якщо симетрична матриця 4 має попарно 
різні власні значення, то матриця В буде також симетричною, причому трьохдіагональною з 
позитивними елементами /,, , 1/3, ,, К - 2,3,..., п, розташованими відповідно під і над 


головною діагоналлю /3,,, Кк - 1,2,..., п. Далі, з подібності матриць А і В витікає рівність 


п п 
Уа - У, Зі, 
ра 1 


що можна використовувати як контроль при обчисленнях. 


З виразів (3.20) з посиланням на рівність (ці, Аці,..., А" ц) з Є (цу и,...,й, 
робимо висновок, що ортонормована послідовність векторів и, и,,...,и, отримується у ви- 
гляді 

иа оч ФОАці; Кк -0,1,...,п, (3.20 А) 
В Кі 
де 
ФІ) -1 
1 
Ф (2) 27? - ФС) 
12 
1 
фі) 35ЧЧ11-Во)Ф()-Во ФО, 
Віз 
зіазаньнаніні боби баіанінібаньнимо б о борні іо (3.22 
1 
ФА оНЧО В Ф А) В о Фо) 
К-1К 
Ф (2) те (1 Вид Ф РН 1 9 зак є ДН Ф лобі)» 
каз ТЗ, 
причому при нормальному перебігу процесу обчислень, коли вдасться обчислити всі вектори 
лінійної оболонки 67 (и, из,...,и,), поліном Ф (4) співпадає з точністю до деякої постій- 
ної з характеристичним, Відмітимо, що сім'я поліномів Ф.(Л), Ф/ (2), ..., Ф,(А), обчислю- 


вана, як і будь-яка інша сім'я ортогональних (по деякій мірі) поліномів, по тричленним реку- 
рентним співвідношеннях (3.22), утворює в силу додатності коефіцієнтів Й, ,,, в, 


Кк -2,3,..., п, послідовність поліномів Штурма. 


3.3. Метод А - мінімальних ітерацій 169). Метод А - мінімальних ітерацій полягає в 


побудові сукупності 4 - ортогональних векторів на лінійній оболонці 65 (4,, А4,,...ь 4"4,), 
Кк -1,2,...,п, послідовних ітерацій деякого початкового вектора додатньо визначеною мат- 
рицею. 


З визначення додатньо визначеної матриці слідує, що матриця А невироджена 1, отже, 
для неї існує обернена 4", яка, очевидно, теж додатньо визначена. З кожною додатньо ви- 
значеною матрицею, діючою в евклідовому просторі 57"), можна зв'язати скалярний добуток 
двох векторів х, у є 52" 


(АХ, у), 
який на відміну від старого позначатимемо квадратними дужками, вважаючи 
Го У1 5 (АХ, у). (3.23) 
Легко переконатися, що скалярний добуток (3.23), який часто називатимемо А - скалярним 
добутком двох векторів х, у є 2" , задовольняє всім аксіомам звичайного скалярного добут- 


70 


ку, заданого в 52". Скалярний добуток |х, у|, також як і (х, у), індукує нову метрику на 
с, перетворюючи його на деякий новий евклідів простір, який позначатимемо через ї2,"? 
. Норму елемента х в цьому просторі позначатимемо через |х| і задамо її наступним співвід- 
ношенням: 


ЇХ "жа ША пар Сх, хе (Сх, ), хо є КУ, 


в 


(3.24) 
Будемо називати її А - нормою. 
Спряженими відносно матриці А векторами (або А - ортогональними) назвемо від- 
мінні від нульового вектори х, у є ї2" , що задовольняють рівності 
(АХ, у) «(х, Ау) «0, (3.25) 
А - ортонормованими називатимемо вектори х,, у, з 82, для яких має місце (3.25) за умо- 


ви, що 
Х 





хо- Й 


З 
5 77 





ІІ 


х, у - - відмінні від нульового довільні вектори ст), 

Після цих зауважень перейдемо до розгляду обчислювальної схеми методу А - міні- 
мальних ітерацій. Як і раніше припускатимемо, що початковий вектор а вибрано вдало для 
його послідовних ітерацій. Нехай, починаючи з А - одиничного вектора 


й м « (4, 4), 


Л ж 
, ч- (4154125 з Чи) . 
побудуємо по тричленним рекурентним співвідношеннях систему А - ортонормованих век- 


аг 
4 
торів таким чином: 





Б4» 2 -РиЧ Кк 44, 
Роза З У Фо В А), 
Раача З Фо» З Роз 3 43» 


паля Мед ФО ок МЕ ал (3.26) 
і личи З Фулнобиоо були 441 
09 бу У 4, 
що, очевидно, можна переписати у вигляді матричної рівності 
АО-ОВ, (3.27) 
якщо ввести позначення по аналогії з вище наведеними. А саме, покласти 
Ї415 45» --:»4,1- 0, 
ГА4, 44»... 44,1- 40 
КУ ОУН ДИ 5 ЗИ РА 0 0 | 
зі, з еф 0 0 0 
0 В 32 В 33 В 34 0 0 
вда небо» Моб забава о» зно 0 Зьнодох "пе 
0 0 00... Вол. бле ах 
0 0 ФО дао 0 оо пабі 








Анулюючий матрицю 8 поліном Ф (Л) будується точно так, як і відповідний полі- 
ном для симетричної матриці (див. (3.22)), що витікає з співставлення схем (3.20) 1 (3.26). 
Тому маємо 


Ф.(І)-1 
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1 
Ф()- ь ЩОо-ВОФ 1, 


12 


1 
Фь()- р. їЕ-вьа)Ф (7-6, ФІ), 


23 


АРЕНА РИ РРО НРА РАНИ ЛЕО РЕ ТЕНМ п (3.28) 
1 
ФІЗ 5 Фу) Ф 203 Було ФиЗ()Ь 


к-ІК 
Ф,(2) -(1-5,,)Ф,1)-5,а Фо); 
Коефіцієнти ф,, поліномів ф(І), Ф (І), ...Ф,(1) В даному випадку обчислюються по- 
іншому. Втім їх визначення також супроводить побудові сім'ї А - ортонормованих векторів 
4» і-1,2,..., п, які обчислюються індукцією по Кк. Припустимо, що перші К векторів лі- 
нійної оболонки а У о РИН у вже побудовані так, що 
За» Чо 4) 9 За 44. 4 4,)- Тоді вектор 4,,, будуватимемо так. Складемо лі- 
нійну комбінацію 
Чка З "Фусабоа З бук Чу (3.29) 
коефіцієнти Ф,, ,, Б, якої визначимо з умови А - ортогональності вектора 4;,, З кожним з 
векторів оболонки 9 (4, 455.-:54,)- Для цього помножимо 4;,, на матрицю і запишемо ви- 
раз, що витікає з (3.29): 
4» 44) 9 09 ву» 4) бу (4 44.) А"4,), 
звідки 
Бека З 44 44) 22,3, п; 
в. «4. 14.) ю-Ї2 й. (3.30) 
Визначивши коефіцієнти Ф,, ,, ф,, за формулою (3.29) вираховуємо 4.,,, котрий при вдало- 
му виборі початкового вектора 4, може виявитися відмінним від нульового. Якщо 4, 70, 
то в якості (К -- 1) - го вектора сім'ї, що будується, беремо 


! 
Чина 














Чеа ЗТ, |» (3.31) 
КА 
вн о » Каі,2, ...ьП-Ї. 
В результаті таких обчислень можна побудувати всі вектори сім'ї 4,), К - 1,2,...,п, 


визначити елементи ф,, додатньо визначеної тридіагональної матриці В і отримати зобра- 
ження (3.27). 


3.4. Метод біортогоналізації послідовних ітерацій початкових векторів матриця- 
ми А і 4 |37, 69). Метод ортогоналізації послідовних ітерацій, запропонований К. Ланцо- 
шем (116, 117|, можна розглядати як узагальнення його методу мінімальних ітерацій на ви- 


падок несиметричної матриці А: їЖ") -» хЖ"?, Суть методу полягає в тому, що при вдалому 
виборі початкових векторів и, 1 у, будується не одна, як в розглянутих вище методах, а дві 
послідовності (и,) і у.) спряжених між собою векторів и, 1 у,, і - 1,2,...,п, кожна з яких 

- . - . ж ж - 
знаходиться в біортогональному відношенні з однією з оболонок б (у, 4 у... (4 ) м) 1 
З (и А и), К -1,2,...,п-1. Розгляд двох замість однієї спряжених послідовнос- 
тей векторів обумовлений наступним відомим з аналітичної геометрії фактом. Всі прості 
співвідношення між координатами векторів (у вигляді скалярних добутків), заданих в прямо- 


кутній системі, можна зберегти при переході до косокутної системи, якщо кожний вектор 
представляти в двоїстому базисі двома системами спряжених одна одній координат (які в 
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прямокутній системі між собою співпадають). Задаючи вектори - - послідовні ітерації векто- 

ра несиметричною матрицею - - в двоїстому базисі (и,) і (у,) , ми можемо отримати таку ж 

просту по своїй структурі матрицю В, як і у випадку ітерацій симетричною матрицею. 
Дійсно, вибудовуючи дві послідовності векторів за тричленними формулами вигляду 


ца ЗО рей ЩЕ Ми, 
ж 
Уа Вікна ПВ ОО Уа (3.32) 
Кк а2,3,...,п, 


ць в-Є ци, 
Фф 
у зав тАУ, 
і визначаючи коефіцієнти г... 1 8. з умов біортогональності векторів и. 1 у, і -1,2....,п 
І/ 1, і і? ев; 2 Ро 


одразу одержуємо 
- (ЯигзУ,) (цу » 4"»,) 





ат з еВ (3.33) 
(цу У,) (пу Ук) 
Що ж до коефіцієнтів 
а а Лота? о (Ауді) 
Кк » Кк 2 
(мзУраі? (мзУраі? 


то вони теж рівні між собою, що витікає з наступного ланцюжка рівностей 
(Ац Уа) З (цу узі) З (ЗО уро З лій З ЧУ В лока В ельлУл) З 
За о о (ог Ура) ен ель Ук (А УА В о Уро) Р 
яВ ль о Ук) 2 (Ми Ук) (А У му). 
Далі, оскільки 
(А'у, ца) З (Ай Урі) З (цу, Аа) з (цьУТфульоУка ЗВ клік) У (Муз Ущ), 


то, вибираючи орієнтування векторів біортогонального базису так, щоб було (и,, У, )» 0 для 


будь-яких значень ХК - 1,2,..., п, отримуємо остаточно 
(цу У, ) 
Ват ана о но оо20, (3.34 
(пу) 
К2,3,...,п. 


Таким чином, для визначення по формулах (3.33) 1 (3.34) координат Ф,, ,, Б,, векторів - іте- 
рацій в біортогональному базисі (и,), (у,) вимагається вирахувати послідовно скалярні до- 
бутки (и, у» Урі)» (щьу) 1 (ц ду, -(Ац,,у,) для К -1,2,...,п (при К «1 перший ска- 
лярний добуток не обчислюється). Самі ж вектори и,., і У, , К - 1,2,...,п, потім знахо- 
дяться по рекурентним співвідношеннях (3.32). Процес обчислень обривається, якщо при де- 
якому значенні К - Кк, одержуємо (и, , У, ) - 0. Вважатимемо, що якщо К, - п--1, то процес 
обчислень протікав нормально. Якщо ж (, «п--1, то процес вироджується. Можливі декіль- 
ка випадків виродження: 

2 и, ЗУ, 70 (двосторонній обрив) 

а и, 7 0, у, 2 0 (односторонній обрив) 

31 и, 7 0, у, 7 0 (односторонній обрив) 

А и, 7 0, у, 7 0, але (и, , У, ) - 0 (тупиковий обрив). 
Односторонні і тупиковий обриви процесу обчислень являються винятковими і їх можна 
уникнути за допомогою належного вибору початкових векторів и, і у,. Двосторонній обрив 


аналогічний достроковому припиненню обчислень в інших методах, причини яких обгово- 


рюються нижче. 
Виходячи з виразів (3.32) - - (3.34), можемо записати 


ек. 








ацафуцти, 
РРО РЛЬ РАСЛЬ ТИ 
РЕ ДУ р Р (3.35) 
Айр їз В ін оїа у РАВ 7 чи, 
Ми, о В ів аа Ав 
для векторів з оболонки 69 (ц, Ай... А" ц)ії 
А уеРбуіу бу» 
А у "Вал Я зга КУ» 
АР УР не (3.36) 
А Упоі - В поТп-2 Уп-2 а пої яд У, 
А М я І Йой В Ул 
для векторів з оболонки 69 (у, А у. (А )" м). 
Введемо зрозумілі вже нам позначення 
Гц з Мо, М, ТЕО, 
ГАц, Ац.,..., Ли |з АС, 
Ги з Узз-» У, ТЕ У, 
ГА" у 4 у,,..., 4 у, | - А"Й, 
Ви 1 0 0 1'анам 0 0 0 | 
каз рен 1 0 ..... 0 0 0 
0 12 ЗОНИ 2 СРО ПУ 0 0 0 
фарб фбраж зачой- яівюофо о 04 Фаьфі 04звіем За , 
0 0 0 0 ко у; ПЕ Й, о І 
0 0 0 0 ..... 0 Ві «и 
тоді співвідношення (3.35) 1 (3.36) можна записати так: 
АСеуВо АРзУВ (3.37) 


Анулюючий матрицю В поліном Ф (7) будується по тричленним рекурентним фор- 
мулам відповідно до схем (3.35) 1 (3.36). 


3.5. Визначення власних векторів |69). Розглянемо спочатку побудову власних век- 
торів за даними методів ортогоналізації послідовних ітерацій і мінімальних ітерацій 
К. Ланцоша. Вважатимемо, що корені характеристичного полінома Ф (Л)-0 знайдені 
(яким-небудь з методів |8, 37) 1 складають сукупність (Л,), і- 1,2,...,п. Для кожного з 
власних векторів и"? , що відповідають власним значенням 4, :з мОЖна записати зображення 
його в ортонормованому базисі (и,) у вигляді 


и - У сан,» і-1,2,...,п. (3.38) 
Каї 


Компоненти с; 2 (и? и, ), вектора и? в базисі (и. можна знайти таким чином. Запишемо 
розклад початкового вектора и, ортонормованого базису по власних векторах и"), 


і-1,2,..., п.Матимемо 


п 
і (2) 
ця У аи. 
ізі 
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Пригадаємо тепер формули (3.10 А) ї (3.20 А) ї для всіх подальших векторів и,, 
Кк 2,3,..., п, базису знайдемо 


1 п. а, , 
ро оФі( 4) ц з Уфа ди"). 


к-ІК ізі К-їК 
Звідки зразу ж отримуємо 
Сто 
В зу 


Попутно замітимо, що значення поліномів Ф, (4) можна обчислювати безпосередньо по 


і, з (и), и,) ві 


"КО (339) 


відповідним рекурентним формулах (3.11) 1 (3.22). 

Таким чином, підставляючи значення (3.39) коефіцієнтів с;, в зображення (3.38), оде- 
ржуємо з точністю до деякого постійного множника с, (тобто з тією ж точністю, з якою 
можна отримати їх, розв'язуючи систему (3.2)) 


п 
ці? о УР, ендо о (3.40) 

Ка 12814 
У випадку нормального перебігу процесу обчислень по методу А - мінімальних ітера- 
цій власні вектори и? додатньо визначеної матриці визначаються так само у відповідності 
до приведеної схеми з тією тільки відмінністю, що в якості базису 42" вибирається сукуп- 


ність (4,) 4 - ортонормованих векторів, що задовольняють умовам 
(4 44,)2-бі;» 5)212,...,п. 
Звідси з точністю до деякої постійної одержуємо 


М А, 
МО у СОР зар) 


Езі к-ІК 





По аналогії з приведеними обчисленнями так само просто визначаються власні векто- 
риб іо 
ємо 


несиметричних матриць А і А" за даними методу біортогоналізації. А саме, ма- 


Й п А. 

ЕЗі (ну) 

| М А. 
кеї (им) 


Ко 


о ізі,2,...,п. (3.42) 


3.6. Розв'язок систем лінійних рівнянь. Наявність даних, які отримуються за допо- 
могою методів, описаних в попередніх пунктах, дозволяє знаходити також і розв'язки систем 
лінійних рівнянь 


Ахаб 
А: АЮ ж рес, (3.43) 
Дійсно, здійснимо перетворення 
рер (3.44) 
де С -|цні,ц,...,и, | -- ортогональна матриця, складена з векторів-стовпців м, 
і - 1,2,..., п, отриманих реалізацією або методу ортогональних ітерацій (п. 3.1), або методу 


мінімальних ітерацій (п. 3.2). Помножимо рівняння (3.43) зліва на матрицю ЇЛ і з врахуван- 
ням рівностей (3.12) або (3.21), а також підстановки (3.44) отримаємо 


ВЕЗПГЬ. (3.45) 

Ясно, що система (3.45) простіша початкової (3.43), особливо ця простота проявляєть- 

ся у випадку симетричної матриці А. Тоді, як ми пам'ятаємо, матриця В буде тридіагональ- 
ною. В останньому випадку систему (3.43) можна розв'язати за допомогою одного витонче- 


75 


ного прийому |69|, котрий можна розповсюдити на системи з несиметричною матрицею. 
Опис такого узагальнення не зустрічається в роботах |8, 12, 15, 37, 68, 69). Проте, він може 
виявитися важливим для застосування. Тому приведемо тут цей результат. 


Теорема 3.1. Якщо початкові вектори и, і у, є Ау,, ЦІ, У; У, Є 2, і структура не- 
виродженої несиметричної матриці А: ХХ"? -» їж? такі, що забезпечують нормальний 
перебіг процесу обчислень по рекурентним формулам 

иц з Аца фун СВ ло» 
ж 
у "4УнВучнана ПВ уко ко» 
Кк а2,3,...,п, 
ж 
и, з Ац-Виц, у, 24 у,-Фиу! 


векторів двоїстого базису Чи,), (у, в 427, де 


| (цу) ке12 
КИ. , те з--П, 
(шу) 
(муз У) 
В а- у вно 
(пул Ук) 
то розв'язок системи лінійних рівнянь 
Ахаб, 


який будемо шукати у вигляді 


п 
Ха 26 ої; з 
рак 


може бути отриманий за п послідовних кроків наближень по схемі 


() (В, ув) 
х''ясочо М, 
(щу) 
Р 
хін щ 9 ре Ук) , 
(цу У.) 
кі 
варена" хУ -Убін, 
іч 
Ка2,3,...,п. 


Доведення теореми витікає з наступних несуперечливих міркувань. Оскільки справед- 
ливі вирази 


(у,)з 2 боб ЙО -б ие КУБОьноо й, 
зі 


К-1 
(9 т ра рен 
аб удреУ Єицц,уд)єй, Е22,3,...,п, 
ізі 
) задовольняє умовам 
юЮ 
х-х? є цу Моно 4 Ху У мнь Мі) 
Тому для К - ої компоненти вектора-розв'язку х в базисі (и,) можна записати 
щ- (9 щі (ю й що 
5 «и УдД-(хоХ У дУ(хох УВУ СВ лька) З 
- (юЮ й реа о 
«(кох ,4уг) з (Ура) Ка2,3,...,п, 
де п з Ах- Ах «Ь- Ах"? Звідси одержуємо, що 


- (пу Му) 
і (піз У) і 


то вектор х-х" 


Ка2,3,...,п. 


Отже, 
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Ура) М 
о й 
(ць У) 
що з врахуванням визначення першого наближення х"? у відповідності з вибором початко- 


вого вектора у; - Ау, повністю доводить теорему. 
1 0 
(юЮ 


зно а з (3.46) 


Помітимо, що нев'язку вектора х"" можна вираховувати рекурентно за формулою 


(Пе Мк) 

роснорунео со Й 

Ка Кк Я 
(щу У.) 

яка виводиться елементарно з посиланням на (3.46) 
(поз Укл) 
да зрейх ЗВ'авеЛІ Зоо ож ; 22,3,...,пчі. 
ць Уцг 


На закінчення цього пункту відзначимо, що за наявності даних методу А - мінімаль- 
них ітерацій (п. 3.3) система лінійних рівнянь 
Ахаб (3.47) 
(з додатньо визначеною матрицею) розв'язується простіше, ніж в попередніх випадках. На- 
справді, при реалізації методу А - мінімальних ітерацій ми одержуємо матрицю 
0-44 ноЧВ 
складену з А - ортонормованих векторів - стовпців 4, 2 (4115 доро...» 4,,) і тому має місце 


рівність 


ОАО-Е"?, 


Звідси, якщо покласти 


ході, де? 


і помножити систему (3.47) зліва на матрицю О", то негайно матимемо 
ж ж 
О 4Ос-сзОБ 
Отже, розв'язок системи (3.47) визначається у вигляді 
ж 
хаоОО0 Б, 


що, між іншим, дає для 4 -О0). 


3.7. Вибір початкових векторів і побудова базису в 52. Розглядаючи методи не- 
прямого визначення коефіцієнтів характеристичного багаточлена, анулюючого матрицю 
(див. пункти 3.1 - 3.4), ми кожного разу підкреслювали, що нормальний перебіг процесу об- 
числень залежить не тільки від структури початкової матриці, але і від вдалого вибору поча- 
ткового вектора (або пари векторів, якщо йдеться про метод біортогоналізації). Для того, 
щоб успіх завжди супроводив обчисленням по схемі того або іншого методу з'ясуємо можли- 
ві причини припинення обчислень і намітимо шляхи виходу з скрутних ситуацій. Втім, осно- 
вна причина припинення обчислень ясна. Це -- лінійна залежність системи векторів, отри- 
маної ітераціями матрицею початкового вектора. В свою чергу для виявлення зв'язку харак- 
теристик цієї системи векторів із структурою матриці 1 з вибором початкового вектора вико- 
ристовуємо спеціальні базиси в 57"). 

Сконцентруємо нашу увагу перш за все на алгоритмах знаходження характеристично- 
го полінома для симетричної матриці. Ці алгоритми простіші інших в тому відношенні, що, 
як відомо з лінійної алгебри, характеристичні багаточлени для симетричних матриць завжди 
розв'язувані в полі дійсних чисел. А це означає, що для всякої симетричної матриці існує 
(над полем дійсних чисел) хоча б один власний вектор. Втім кількість цих векторів не менше 
ніж кількість простих власних чисел. Це витікає з визначення (3.2) власних векторів. 

Нехай и? є 22?) -- власні вектори симетричної матриці А, що відповідають простим 


власним значенням Л,, і - 1,2, 3,...,К « п. Розкладемо початковий вектор р? (в методі мі- 


німальних ітерацій р? зи, в методі А - мінімальних ітерацій р? «4,) по базису и?) з 


ЛІ 


власних векторів матриці А 


К 
р«Ура'зро 348 
ї«1 


де Р - (п-Ю - мірний вектор, що належить доповненню 52" 4 6и),...,иО) лінійної 


(5 


оболонки Єи?,...,и?) до всього простору 42". Відмітимо, що якщо серед власних чисел 


2, є кратні, то А « п. Для ітерацій вектора р? матрицею А одержуємо 


Ко 
Ар З роати коси; (3.49) 
ізі 
де К, 4 К, якщо серед власних значень є рівні нулю. Зображення (3.48) 1 (3.49) виписані для 
найсприятливішого вибору вектора, коли р,; 70, і - 1,2, 3,...,к - п. Проте може трапитися, 
що навмання вибраний вектор р"? належить Є(и"?),...,и?), І « К, і тоді отримаємо 


Кк, з пайо(, Ко). 
Таким чином, ітерації 

УЛ Гор РУ па о (3.50) 
вектора р(? симетричною матрицею А будуть лінійно незалежними тоді і тільки тоді, коли 
серед чисел 2, і-1,2,..., п немає кратних (К є п) і рівних нулю (К є Кк еп), авектор р"? 
вибраний так, що жодна з компонент р,, із-1,2,...,п не дорівнює нулю. Якщо ж 
р би), ...и?У, (І « п), то в послідовності (3.50) незалежними виявляться тільки век- 
тори А"р?, т-0,1,...,1-1, тоді як вектори А"р'?, те 1,1--1,...,п-1, будуть лінійно від 


Ф або для К-.2, якщо 


них залежати. Це трапиться вже для І-І, якщо р" зи 
о ДР (1 0) 
ро з-ри Ярий 


вибору початкового вектора р"), але принаймні з'ясували, від чого залежить припинення об- 


і т. п. Цими міркуваннями ми не знайшли ще критерію для сприятливого 


числень. Це може бути зв'язано або з наявністю в спектрі" матриці кратних власних чисел, 
або чисел, рівних нулю, або з тим, що початковий вектор належить лінійній оболонці влас- 
них векторів матриці, не співпадаючої зі всім простором. 

Повернемося до рівностей (3.49) і відзначимо, що за наявності в спектрі матриці одно- 
го або декількох власних чисел, які перевершують по модулю інші, ці великі по модулю чис- 
(2 


ла, складові р, Аг и"", що містять ці більші за модулем числа, виявляться переважаючими. 


Звідси виводимо наступне правило вибору початкового вектора р"). Беремо довільний век- 
тор, наприклад, вектор - праву частину системи рівнянь, якщо така задана, або яку-небудь 


сторону ау З (45 425 г» ду) Заданої матриці і послідовно обчислюємо 
Ра З 44» Роз З Раз» ит З Фан): 


Якщо т «5-6, то найбільші по модулю власні значення вже гратимуть основну роль і век- 
тор Рі) буде близьким до лінійної оболонки, натягнутої на власні вектори, що відповідають 


великим за модулем власним значенням. Нехай серед компонент р,;, і - 1,2,..., п, вектора 


Фипт) КОМПОНЄентиИ Р,» Ї З 1,2, 3,...,1 5 п, виявилися "малими" за модулем. Змінимо їх дові- 


7 Спектром матриці як правило називають дискретну множину ІЛ.) коренів її характе- 
ристичного полінома. 
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льним чином, додавши до них деякі числа у ік В результаті отримаємо 
2 ри о ж 
Рию З (Фоі» з Фан) 
-, Риї 2 І- ї, 
Риі ру ЧУї» іт» К1,2,. 
Рі У І-Ц, нь ль КАЛ 


Далі, обчислимо лінійну комбінацію 


т-1 


Р п Уу іа) тр 
ї-1 


при заданих у,, і-1,2,..., т-1, і в якості початкового вектора р"? беремо відповідним 
чином пронормований вектор р"). 

Проте такий вибір початкового вектора зовсім не гарантує того, що процес обчислень 
в будь-якому з методів ітерацій (див. пункти 3.2 - 3.3) завершиться нормально. Тепер його 
протікання визначатиметься структурою матриці. Нижче для стислості поведемо мову тільки 
про процедури, регулюючі процес обчислення в методі мінімальних ітерацій. Якщо обчис- 
лення за цим методом завершуються достроково на Кк - ому кроці при К «и, то це означає, 
що ф,(2) буде мінімальним анулюючим вектор р" - и, поліномом. Слідуючи пропозиці- 
ям 169), покажемо, як можна добудувати систему векторів 

йо пвойво 

до ортонормованого базису всього простору 47", якщо вектор и, виявився лінійно залеж- 


ним від векторів оболонки б(и,,...,и, ,), тобто, якщо на К - ому кроці ми отримали 


ко 
Ос -У В раз б Ж А 
рик 


Вчені, кп. 
Продовжувати обчислення будемо за наступною схемою. 
Виберемо вектор 57, що не належить Є(и,,...,й, |), У вигляді 
к-1 


ВО ех и У пе йо -| 


із! 
де у, і - 1,2,..., К, задані числа, а е, - - одиничний вектор з компонентами рівними нулю за 


0) ТМ 59), 


виключенням Х - ої, рівній одиниці. Розкладемо вектор ро на дві складові 
(2 : 2 
ре «тн 1 ро Є Оу йо 


ів якості и, послідовності, що будується, беремо вектор 


0) 
о Рі 


" Чреї 
0) 


Проекцію р" відшукуємо тим же способом ортогоналізації, який застосовувався для побу- 





дови всіх попередніх векторів и, и,,...,и, | Наступний (К -- 1) - ий вектор обчислюється 


точно так, як і визначався вектор и., тобто знаходиться спочатку 


' 
ца З7Риц тб 4й, 


і 2д, 17, 

Ітерації вектора и, матрицею 4 можуть виявитися знову лінійно залежними на де- 
якому / - ому кроці обчислень. Якщо / « п, то доведеться повторити описану вище проце- 
дуру, починаючи з вибору нового вектора Р" і т.д. Для того, щоб завершити обчислення 
можливо доведеться побудувати декілька лінійних оболонок 


В о зорбівосо п ЕРЖИНИМО яру пря Й т (цим) 
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ортогональних одна одній. При цьому весь простір 42" 


пряму суму цих оболонок 
зна пе (уусарібр а Я 8 ко (щу ПАР У ДБ РА Й ут (ооо їВ.) (3.51) 


А оскільки кожна з цих оболонок отримана на підставі ітерацій деякого початкового вектора 


Б, із1,2,..., т« п, анульованого відповідним поліномом ф "(г), то характеристичний 


розвивається природним чином в 


поліном матриці А з точністю до деякої постійної співпадає з добутком мінімальних полі- 


номів Ф З) , що анулюють послідовні "початкові вектори" р , і-1,2,..., т, тобто 


ск. (7) 
Ф «0) т По (г) , 
ізі 
де кожний з поліномів ф "(г), відповідний вектору р, знаходиться за схемою (3.11). Ви- 


пишемо для огляду матрицю В переходу від системи векторів ші) до системи ГА» и , що 
відповідає випадку 
(по - 
ОЗ У УПА Уч 3 С У 7 УНН У АТ 


Матимемо 
і би б о о о о о те о о 
Взр Во» Ваз 0 0 0 0 
0 0 0 Вузюз Вуоюо Врач 0 0 0 0 0 
В"- 0 0 0 0 Вулко Веокчі 0 0 0 0 0 
"Ро 0 0 0 0 0 Віко Вик 0 0 0 
0 0 0 0 0 0 Верн Вени Врзткчо 0 0 
УР С з ЛК дай і Ж б, "УЗВАР. РАН ВИ 
ДО РР Ми а М МУ 0 0 ооо «оо 





Як бачимо, в описаному випадку матриця В складається з двох тридіагональних кліток, ко- 
жна з яких співпадає по своїй структурі з матрицею перетворення (3.20). Очевидно, в загаль- 
ному випадку, відповідаючому розбиттю (3.51), матриця В складатиметься з т тридіагона- 
ЛЬНИХ КЛІТОК. 

Побудова базису а) в ЖЖ") і відшукання коефіцієнтів характеристичного багаточле- 
на Ф, (А) - 0 за допомогою методу А - мінімальних ітерацій здійснюється майже у відпові- 
дності з описаною вище схемою за винятком однієї обставини, коли замість ортонормованої 
будується А - ортонормована сукупність векторів. 

Описаний спосіб вибору початкового вектора і схема добудови базису в 47" легко 
узагальнюються на випадок несиметричної матриці. Дещо складніше обгрунтувати це уза- 
гальнення, оскільки не кожний характеристичний поліном (з дійсними коефіцієнтами) неси- 
метричної матриці вирішуваний в полі дійсних чисел і в зв'язку з цим не кожна несиметрич- 
на матриця має хоча б один власний вектор над цим полем. За недоліком місця і часу ми не 
торкатимемося тут цього питання, а тільки відмітимо наступне. Процес обчислень в методі 
ортогоналізації послідовних ітерацій (п. 3.1), якщо він обривається достроково на К - ому 
кроці, Кк « п, можна продовжити, вибираючи в якості вектора и, знову підібраний вектор 


р), із1,2,..., тс п, за описаною вище схемою без яких би то не було змін. 

Що ж до методу біортогоналізації (п. 3.4), то при односторонніх або тупикових обри- 
вах обчислення доводиться починати спочатку із знову вибраними початковими векторами 
и, 1 у, Якщо ж має місце двосторонній обрив, то отримана послідовність векторів взаємного 
базису 

і сіаейної 
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У Пд (3.52) 
відповідним чином добудовується. При виборі пари початкових векторів и,, У,, як і при ви- 
борі початкового вектора в методах мінімальних ітерацій, заздалегідь можна обчислити пос- 
лідовності 

2 - Ацау» из - Ай, б НГ Арлу» 
ж 
У З 4ЯУ 


У "Я Ущі о зо з Ус "ЯМ 


п» (те-1)2 


й з б й б ж 5 
де в якості и, 1 Ук, МОЖНА, наприклад, взяти відповідно який-небудь рядок ду і стОВПеЦЬ 


а) 
аю матриці А і на цій основі визначити и, і у,. Добудова системи (3.52) починається з ви- 


бору таких векторів р" 147), і-1,2,...,т, т« п, які задовольняють умовам 
р? є и цоцуа) і РОДУ, уз Ура) 
49 є (уруьуці) і 4916цу му) (3.53) 


Якщо на К - ому кроці має місце двосторонній обрив, а вектори р"? і 40? підібрані відпові- 


(5 


дно до вказаних вимог, то вважаємо и, - р), у, «4"). Процес побудови послідовності 


(3.52) продовжуємо до наступного / - того обриву. Він може виявитися останнім, якщо 


Із п. При І « п обчислення продовжуються з новими векторами и, « р", у, «47, які за- 


довольняють умовам (3.53) при К - / іт.д.і т.п. 


3.8. Процеси доортогоналізації при чисельній реалізації канонічних зображень 
(3.12), (3.21), (3.27) і (3.37). В попередніх пунктах описані деякі методи непрямого визначен- 
ня коефіцієнтів характеристичних поліномів симетричних і несиметричних дійсних матриць, 
а також засновані на цих методах способи обчислення власних векторів (при заданих влас- 
них числах) і розв'язків систем лінійних рівнянь з такими матрицями. Ми бачили, що для ви- 
значення коефіцієнтів анулюючого матрицю полінома вимагалося виконати певний об'єм 
обчислень. Так, в методі ортогоналізації послідовних ітерацій доводилося вираховувати лі- 
нійні комбінації (3.6) - (3.7) вигляду 


ї 
14 ра 
ца - 4й, - 5,8 кій; 
ра 


Ви (Аци), Вані | ша (3.59 
Чо ав! і- 1,2,....К; Ка1,2,...,п- і. 

В методі мінімальних ітерацій (п. 3.2) проводилися аналогічні обчислення (див. (3.13), (3.13), 
(3.19)), тільки в меншому об'ємі. Складніші обчислення (див. (3.30) - (3.31) виконувалися в 
методі А - мінімальних ітерацій за формулами 

Чіа З "Фу З Фу Чу 

Бір З (44, 44) Кк -2,3,...,п, 
Ба з (44, 44), К-1,2,...,п, (3.55) 


вка З АС А» Чен)» 


'! 


Че зба» Кк - І, да зано 








І, нарешті, в методі біортогоналізації (п. 3.4) - - за формулами 


ца З Фрі З Фу Ми, 
Ура 2 Фрі 7 Фк /к З А"у, (3.56) 
А 2 З 
о | м), Кк -1,2,...,п; ПЛР С м, , Кон 
(цу з Ук) (цуол» Урі) 


Далі, для відшукання власних векторів матриць у вигляді (3.40), (3.41) 1 (3.42) перш за 
все повинні бути визначені коефіцієнти /8,,, б,, по відповідних формулах (3.54) - (3.56), а 
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для отримання розв'язку систем лінійних рівнянь необхідно мати крім цих коефіцієнтів ще 1 
значення векторів и,, У, Ац,, Ду, , ч.» і- 1,2,..., п, по можливості більш точні. Проте в 
зв'язку з тим, що при реалізації на ЕОМ майже кожної арифметичної операції неминуче ви- 
никають, а при великому числі операцій нагромаджуються помилки заокруглення, немає ні- 
яких підстав сподіватися на отримання добрих результатів при реалізації схем (3.54) - (3.56). 
Це підтверджується дослідженнями |12, 15, 33, 37, 68, 69, 95|, де також сказано, що для 
отримання з деякою гарантованою точністю результатів типу скалярних добутків (при мно- 
женні матриці на вектор), необхідно спеціальним чином організувати весь обчислювальний 
процес. Зокрема, щоб отримати мінімальну похибку заокруглення суми деяких чисел, необ- 
хідно підсумовувати ці числа, ранжирувані за величиною їх модуля, починаючи з найменшо- 
го (за абсолютною величиною). Надалі, не зауважуючи кожного разу особливо, широко ко- 
ристуватимемося "впорядкованим підсумовуванням". Втім, використання одного тільки 


"впорядкованого підсумовування" не завжди приводить до бажаного результату. Так, напри- 
(ю 


клад, вже при нормуванні вектора и"") за схемою 
и ої? ї ї - 
фер Їмо а (и? бу а У ші, (3.57) 
із 








часто не вдається досягти рівності |12, 95|. Для того, щоб обчислити орт и у напрямі вектора 
чі г пі п ж . - . 
п з (й, йо,..., й,) , потрібно, як це рекомендується в |95|, провести циклічні обчислення за 


схемою (3.57) для К - 1,2,..., т, починаючи обчислення з вектора и"? - її і закінчуючи їх 








отриманням вектора и, - и, для якого має місце нерівність І, и» -1| «є, де є» -- ма- 


шинна точність зображення чисел. На приведеному прикладі проілюстрована основна ідея 
обчислень того або іншого результату: для відновлення втраченої при заокругленнях точнос- 
ті початкових даних потрібно будувати деякий (швидкозбіжний) ітераційний процес з враху- 
ванням непокори машинної арифметики звичайним законам комутативності, асоціативності 1 
дистрибутивності. Ідея використання послідовних наближень (з вказівкою необхідного числа 
ітерацій), направлених на відновлення значущих цифр шуканого результату, являється осно- 
вним, якщо не єдино можливим, засобом розв'язку поставленої в даному розділі задачі. Без її 
реалізації не можна визначити жодним з методів, розглянутих в пунктах 3.1 - 3.4, коефіцієн- 
ти характеристичного полінома з високою точністю. 

Першим кроком в реалізації даної ідеї є побудова за допомогою відповідних ітерацій 
(донормування по |95|) одиничних або А - одиничних векторів, що дозволяють більш надій- 
но визначити як самі вектори и, і 4,, так і коефіцієнти /),, і Б,, рекурентних зв'язків між 


ними у вигляді (3.34) - (3.56). Ми з самого початку передбачили застосування цієї процедури 
і тому подали опис методів в пунктах 3.1 - 3.3 на відміну від традиційного їх викладу в робо- 
тах |8, 37, 69| з врахуванням нормування системи векторів, що будується. 

Другий крок полягає в побудові ітераційного процесу для обчислення ортогональної ( 
А - ортогональної або біортогональної) системи векторів. Описані в пунктах 3.1 - 3.4 проце- 
си ортогоналізації являються дуже нестійкими по відношенню до помилок заокруглень про- 
цедурами, особливо в умовах поганої обумовленості матриць. В зв'язку з цим нагадаємо, що 


спектр погано обумовлених? матриць належить досить широкому інтервалу, який включає 
як дуже великі, так і дуже малі за модулем значення власних чисел |95|. А це означає, що при 
поганій обумовленості матриці А в послідовності 


що Йо Й Що 2-9. с. тї-Ї 
(або в послідовностях (А"д. ), А" )"»,)) обов'язково присутні близькі по нормі ї2" (або 


52292 У вектори. Останнє в свою чергу може привести до того, що при визначенні з невисокою 


2 Погано обумовленими називають |68, 69, 95| матриці А, для яких міра с( А) - | А: |А7 | приймає значення, 


що значно перевищують одиницю. 
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точністю Кк - го ортогонального (4 - ортогонального або біортогонального) вектора и, лі- 
нійної оболонки би, и.,...,и,) (або їй відповідних в інших методах) початкові вектори 
К-1 


раї 22 В ги; (або Їм відповідні) для К «2,3,...,п--1, можуть виявитися майже лінійно 
ізі 


залежними 1, як наслідок цього, вектор и, не буде ортогональним до лінійної оболонки 
Ви йо... й, |) (або їй відповідних). Тому особливого значення набувають наступні про- 
позиції. 

Теорема 3.2. Нехай початковий орт и, є 47") і структура матриці А: Ж" -» 4Ж"? 
такі, що послідовність 

п-1 
ГРА 17 АРЕНА ца 7 
складається з лінійно незалежних векторів п - мірного евклідового простору 42. Система 
векторів Чи,) буде практично ортогональною еквівалентній системі векторів лінійної обо- 
лонки 
К 
Ви, Ай, д"и,), К-б,1,...,п- і, 

тоді і тільки тоді, коли кожний вектор и,, К -1,2,...,п, буде визначатися за допомогою 
наступного процесу доортогоналізації 
шо що 





у - йти, а ІС 


і-»оо 


о -1 1 1 1 ї 
що УВГ В) ек 338) 


кт? 


К-1 


(і-1) (і) сю (і-1) 
к- ій |- и -У, калу» КОП 
уч 


(0) - що 
шо Но КН Зоо їй 
який для отримання доброго наближення и? вектора и,, що задовольняє нерівності 
(5 
| ри | «є, 


достатньо припинити на І - ому кроці, коли буде досягнута нерівність 


(0) 


ра . 0 . У 
шр|12 "іі и)|, а 5 -- кількість 
І 











Є 
де в якості "порогу" є з г береться число є 4-2 


значущих цифр мантиси чисел (в двійковому коді). 
Теорема 3.3. Нехай початковий вектор 4, є 12" вибраний так, що сукупність 
Упр 1 ПЕНЮ, о 
отримана ітераціями вектора 4, додатньо визначеною матрицею А, складається з лінійно 
незалежних векторів п - мірного евклідового простору 42". Система векторів і4,) буде 
А - ортогональною і еквівалентною системі векторів оболонки 
4» А 4"4,), Кед,1,...,п- 1, 
тоді і тільки тоді, коли кожний вектор 4,, К «2,3,..., п, визначатиметься за допомогою 
наступного процесу А - шо 


з тар, 4074 |4 








19 | 


жо , 


2 (і) 1... (1-1) (і-1) (1-1) 
р зу Фу Чак» 


Б б сб 2- Сн б 44,2» 
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р 0 944); 
БА о 4802, 40), 
49 54. 0 -2,3,..а,п 


який для отримання доброго наближення 40? вектора 4,, що задовольняє нерівності 


Їчк як | 














в г «є, 
досить припинити на І - ому кроці при досягненні нерівності 
5 
ь Чер зи РВРНХ Р ар 35 бо 
РОБКИ Є (0) (0) 
де як "поріг" є ,- йо береться число є 227 са таїп ар 
п 





Я Й ж (п)  . ; 

Теорема 3.4. Нехай початкові вектори и, є ї2", у, є і? і структура матриці 
АС 542? такі, що кожна з послідовностей 

п-1 . ж жу ної 
иь Айьоз А" й і У,АУоо(А) У, 
РО ; ; ; 5 (пу з цк (п) 
складаються з лінійно незалежних векторів п - мірних евклідових просторів 57" і «7  ,де 
ж (п) о . Р . з . 
5" --спряжений з 42" простір. Сукупності векторів Ки.) і у.) будуть еквівалентними 
відповідно лінійним оболонкам 
Кк . ж жі 
Ви Ди А'цу) і Ку, Ау(А у), Кк -б,1,...,п- 1, 


5 ; ра є . - : ж (п) 
1 біортогональними тоді 1 тілЬьКИ тоді, коли кожний з векторів и вза У єї 


7 


к -2,3,..., п, будуть визначатися за допомогою наступного ітераційного процесу 
- йти), - ту, 
і-»ою і-»ю 


ца (1-1) (1-1) (1-1) 
-Вупобо Вул М» 


ре (і-1) (1-1) (1-1) 
"Ву ьо пРульаУнка ЧУк о» 


(і-1) (1-1) 
(і-1) (ур 75 Щоз) собу (щу) 
Коко , катеої т 
(ці 2з Ур-2) (цу КО 
0 0 ж 
ца Ану, Мов Ура» 


який для отримання добрих наближень 


(7) 
Я Р 








(1 
| или Ї в? «є, 


векторів и, і у, досить обірвати на І - ому кроці при досягненні кожної з нерівностей 


(г І 
| коло Му Вл р З 





(г І 
| ко1когУрор в СР клУкеть 





кб 


(0) 
УР) 











М 


Доведення кожної з цих трьох конструктивних теорем відрізняються одне від іншого 
незначними нюансами. Тому достатньо привести хоча б одне з них. Обмежимося доведенням 
теореми 3.2. 

Покажемо, перш за все, що процес доортогоналізації збігається, тобто, що при 
і-1,2,3,..., вектор 


Є 
2 ге 22» 0 
де в якості "порогу" є ,- кі береться число є 2 т, 2 ті | тм 


(2 
ці -зи, 1 (ци, й, ,). 


Позначимо через В,,; матрицю переходу від системи векторів щ,,и,,...5 м |, М") до системи 


цьо Що, яка внаслідок рекурентних співвідношень (3.54) має вигляд 


В,; задовольняють нерівностям 


(р 
м З! 





ої) і 


(2 








ка) 


(2 
к-12 


««1 при / «5-1, і-1,2,3,.... 


(1) 


к-ік | 
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(3.59) 


Оскільки вектори й, іи,, / - 1,2,..., К - 1, майже ортогональні, то елементи /3 ми ; матриці 


Далі, з очевидної рівності 


5 9 ; 
і Їм 








К-1 Е-1 
2 У 
«о пмідува ве СО 
ГА ГА 











витікає, що 
(і) (і) 
о«|я і) | «ВуиЧі (3.60) 
Тому кожна з матриць 8,,, і - 0, І, 2,..., невироджена і має обернену 
Г 1 0 0 | 
0 1 0 
-1 
Ві- 3.61) 
(і) (і) (бу 
ПРАВ С ВЕ Й как | 
у к-1) р. і 1 Р ву 1 001 З 5 б 
де, як неважко отримати У)27- ст при /«к-1і у0,2- . З співвідношень 
К-ІК Ви 
(3.56) виходить, що елементи у (),; оберненої матриці В.; підкоряються нерівностям 
і (і) 
оку кі; (3.62) 
Позначимо тепер через (/""". майже напівортогональну матрицю, складену з векторів - сто- 
впців и, из, й 0), а через С матрицю -- з векторів - стовпців щ,1,...Мо 5 Ай: 


Тоді процес доортогоналізації (3.54) можна записати у вигляді наступної матричної рівності 
Р ж 
п,К -1 
99 «сФф| вт Ро (3.63) 
із0 


і доведення збіжності процесу (3.54) зводиться, таким чином, до визначення границі 


шило ОО 
р-ою 


по нормі А" ,К «п, 








Оцінимо різницю | др: -уво , вважаючи для визначеності, що 42» р. З рівності 


(3.03) маємо 
б. р 
усю -утю а сеї |в" -4 | ві" Ве 
і-о 9 


р р Р 
і б "ег'-ЕФ (п,Ю "а'-ЕЮ9 
- СФЇ | ві," 9 | ві," -ЕОре о а | Ві т о ЕОР. 
із0 ізр ізр 
Звідси, посилаючись на нерівності (3.60), одержуємо 
(пЮ уу(ак) (.ЮЦ | ( р-1уч-ро рю 
|осмо «ро Бе ормо Б сво - ко |. 


Остання нерівність з врахуванням співвідношень (3.62) при 4, р-» со показує, що послідов- 


85 
ність матриць ЇЇ/ ні збігається в собі по нормі простору матриць. Отже, при 4, р-» со 
пого |Аоме ьо, 


тобто числова послідовність Цен || збігається в собі і, значить, обмежена. А це в свою 


чергу означає, що 
БРОУ а Р 


чений 
занормою А? К «п. 

Для завершення доведення теореми помітимо, що в умовах обмеженої точності пред- 
ставлення чисел немає необхідності в отриманні наближень и? з великими значеннями ін- 
дексу і (при цьому обчислення на комп'ютері "зациклюються"). Достатньо обмежитися та- 


ким наближенням и/?, для якого поправка 
К-1 
(0 
2 ку Мі 
ЧІ 


(5 (5 


в будь-яку з компонент ці, р 1,2,..., п, вектора и, не перевищує за абсолютною вели- 


чиною точності її зображення в ЕОМ. Теорема доведена повністю. 
В ході доведення теореми ми знаходимо (зокрема з співвідношення (3.63)), що елеме- 
нти ф,; матриці В в зображенні (3.12) потрібно вираховувати за формулами 


р і-1 
Вуз шву РОЙ ВК єю, 
їзі те 


р 
з Ша Ф? ), Ка12,..,п-1 3.64 
Вин іт ЛП кн п ( ) 


причому раціонально при цьому обмежитися значеннями ро /. Аналогічно визначаються 


елементи матриць В ів канонічних розкладаннях (3.21), (3.27) 1 (3.37). 
Для того, щоб оцінити точність одержуваних таким чином елементів матриць В увсіх 
зображеннях (3.12), (3.21), (3.27) 1 (3.37) введемо до розгляду матриці похибок Д,, Ду, Ду, 


Ло відповідно матриць В, У, У, О і визначимо відносні похибки останніх співвідношен- 


ЛЕ Р рр Р р С | 
|8! 01 ІЙ 10| 
Якщо, виходячи з теорем 3.2 - 3.4, можна сказати, що точність обчислення ортогональної ( А 


- ортогональної або біортогональної) системи векторів практично обмежена лише точністю 
зображення чисел в даній ЕОМ, то вже відносно точності елементів /),, або Б,, матриць В 


нями 














бВеа , ба 











справедливі слідуючи висновки з згаданих теорем. 
Висновок І. Відносні похибки матриць В із зображень (3.12) і (3.21) оцінюються не- 
рівністю 
5В«2с( АС, 
де с(А)- | А І А" | - міра обумовленості початкової матриці А. 
Доведення витікає з очевидної рівності (Г/Л, - АД, - Д,,В. 
Висновок 2. Відносні похибки матриці В з розкладу (3.27) оцінюються нерівністю 


56 «с (О1-«0))0, 


2 п 
де при обчисленні числа с(О) слід враховувати, що | о | з Уа , 4)» 
ізі 





107-Їо"4ї 


Насправді, з виразу ОД, - АД, - Д.В з врахуванням того, що мають місце рівності 


86 
О'А0- Е"?, А М 00, 


відразу ж знаходимо 


А, | «с СОЇ А «(Об0. 


А з матриць А 1 В витікає, що 


14Ї «с(ОЇ| 8| 
звідки знаходимо необхідний результат. 

Висновок 3. Відносні похибки матриці В в зображеннях (3.37) оцінюються нерівніс- 
тю 8 В2 С/)б О, де С(/)-| ЦИ - «(. 

Доведення аналогічно попереднім. 

Таким чином, висновками теорем 3.2 - 3.4 підтверджується відомий в теорії стійкості 
чисельних розв'язків збурюючих систем лінійних рівнянь факт: в питаннях стійкості чисель- 
них методів непрямого визначення коефіцієнтів характеристичного багаточлена матриці ви- 
рішальну роль відіграє міра обумовленості матриці. Чим гірше обумовлена матриця тим з 
меншою точністю визначаються коефіцієнти характеристичного полінома. 


Висновки. В розділі викладені за літературними даними деякі методи непрямого ви- 
значення коефіцієнтів характеристичного полінома матриці, засновані на ортогоналізації де- 
яких послідовностей векторів. Охарактеризовані області застосування даних методів при 
розв'язуванні інших проблем обчислювальної алгебри, зокрема при визначенні власних век- 
торів 1 при відшуканні розв'язків систем лінійних рівнянь. 

В розвиток положень роботи |95| встановлено, що методи визначення коефіцієнтів 
полінома в класичному їх виконанні мало придатні для практичних обчислень на комп'юте- 
рах, особливо в умовах поганої обумовленості матриць. В зв'язку з цим сформульована зада- 
ча визначення з високою точністю шуканих коефіцієнтів 1 запропоновані способи стійкої ре- 
алізації класичних методів на комп'ютері, що враховують непокору машинної арифметики 
звичайним законам комутативності, асоціативності і дистрибутивності. 

Під цим кутом зору дано на відміну від традиційного |З, 12, 15, 37, 68, 69| 1 первинне 
трактування методів, що включає побудову ортонормованих або А - ортонормованих послі- 
довностей векторів. Вказані способи вибору початкових векторів сукупностей, що будують- 
ся, котрі сприяють нормальному перебігу того або іншого процесу обчислень. Попутно дове- 
дена теорема 3.1, яка узагальнює на випадок несиметричної матриці відомий спосіб 
розв'язку систем рівнянь |69|. За допомогою прийомів, описаних в роботі |95|, знайдені стій- 
кі процедури доортогоналізації ортонормованої, 4 - ортонормованої або біортонормованої 
систем векторів, які дозволяють визначити з високою точністю, як самі вектори конструйо- 
ваних множин, так і коефіцієнти характеристичного полінома. Доведені теореми 3.2 - 3.4 про 
збіжність процесів доортогоналізації ортонормованих, А - ортонормованих або біортонор- 
мованих послідовностей векторів. Вказані гарантовані оцінки точності шуканих коефіцієнтів 
на ЕОМ, звичайно, завищені, так як при їх отриманні розглядалися найбільш несприятливі 
можливі погрішності машинних операцій. 


ВИСНОВКИ 

За звітний період досліджень отримано нові результати щодо вирішення проблем ко- 
ректності та розв'язності основних обернених задач теорії потенціалу, в числі яких методи 
розв'язання цих задач для конкретних математичних моделей геологічного середовища, що 
мають як суто теоретичне, так 1 прикладне значення. Основні результати досліджень відо- 
бражені в фахових наукових виданнях та двох кандидатських дисертаціях, захищених за цей 
час авторами звіту. 

Новим в теорії обернених задач для визначення контактної поверхні є введення в нау- 
ковий обіг нового класу єдиності розв'язків лінійних обернених задач потенціалу - класу 
Нумерова, введеного на основі аналізу 1 узагальнення існуючих модельних уявлень і поста- 
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новок контактних задач. Крім того, новим є створення комплексу алгоритмів (які реалізують 
узагальнені способи відновлення контактної поверхні за вхідними даними, заданими як на 
довгих, так і коротких профілях), позбавлених прив'язки до величини середньої глибини Й, 
характерної для відомих методів вирішення структурних задач граві- і магнітометрії. Дета- 
льне опрацювання питань чисельної збіжності та точності скінченновимірної апроксимації 
операторів лінійних та нелінійних задач для контакту з можливістю їх позитивного вирішен- 
ня за вхідними даними на обмежених інтервалах вносить завершальні штрихи в теорію стій- 
кого вирішення вищезгаданих умовно коректних задач визначення форми та щільності збу- 
рюючих тіл. Підсумки порівняння результатів чисельного моделювання розроблених алгори- 
тмів та відомих класичних способів дозволили чітко окреслити місце перших в процедурі 
інтерпретації даних гравіметрії методом підбору та рекомендувати їх у складі автоматизова- 
них систем для експрес-аналізу гравіметричних даних на територіях шаруватої будови. 

Головний підсумок першої частини звіту полягає в обгрунтуванні залежності точності 
обчислень контакту не лише від похибок поля, але і від довжини інтервалу його задання та 
розробці нового методу відновлення контакту на основі інтегральних конструкцій зі швидко 
спадними ядрами, для якого потрібен значно менший обсяг вхідної інформації, ніж в існую- 
чих методах. 

В другому розділі розроблено ефективний метод розв'язування нелінійних інтеграль- 
них рівнянь, що описують відхилення контурів зіркових збурюючих областей від заданих, 
тобто вказано напрямок пошуку нормального розв'язку задачі відновлення границь зіркових 
тіл, близьких до круга. Глобальний регуляризуючий оператор для певного нелінійного рів- 
няння відшукується у вигляді послідовності локальних операторів для відповідних лінійних 
інтегральних рівнянь. В рамках цієї задачі досліджено проблему єдиності розв'язку, а саме 
доведено теореми існування та стійкості. 

Елементи новизни містить метод розв'язання систем лінійних рівнянь, який наведений 
у третій частині звіту. В ньому використовуються процеси доортогоналізації деяких сімей 
векторів, що виникають при побудові характеристичних поліномів матриць. Встановлено, 
що методи визначення коефіцієнтів полінома в класичному їх виконанні мало придатні для 
практичних обчислень на комп'ютерах, особливо в умовах поганої обумовленості матриць. В 
зв'язку з цим сформульована задача визначення з високою точністю шуканих коефіцієнтів і 
запропоновані способи стійкої реалізації класичних методів на комп'ютері. 

На закінчення необхідно відмітити, що 27 травня 2003 року раптово обірвалося життя 
Арнольда Володимировича Чорного - керівника науково- дослідної роботи, доктора фізико- 
математичних наук, професора, головного наукового співробітника відділу глибинних про- 
цесів Землі і гравіметрії. Вчена Рада Інституту геофізики затвердила співкерівниками теми 
кандидатів фізико-математичних наук О. А. Чорну та А. І. Якимчика (протокол Хо 6 від 2 
липня 2003 року). На жаль, авторам звіту не вдалося завершити деякі етапи дослідження. 30- 
крема, передбачалося розробити уніфікований метод доведення глобальних теорем єдиності 
розв'язування лінійних та нелінійних обернених задач теорії потенціалу, створити ряд мето- 
дів розділення сумарного потенціального поля на складові, за допомогою якого розроблені 
способи пошуку збурювальних тіл за відповідними аномаліями можна було б поширити на 
складні реальні ситуації. Заплановані та розпочаті Арнольдом Володимировичем теоретичні 
дослідження автори цього звіту безперечно будуть в силу своїх можливостей продовжувати. 
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